Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,162,830 A→G 5.5% N770D (AAT→GAT)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,162,8300AG 5.5% 204.7 / 3.5 73N770D (AAT→GAT) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (2/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

AAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGT  >  AP009048/2162761‑2162900
                                                                     |                                                                      
aaaGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCaa                                                                       >  1:1714709/1‑71 (MQ=255)
aaaGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCaa                                                                       >  1:1720863/1‑71 (MQ=255)
  aGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATa                                                                     <  1:5225105/71‑1 (MQ=255)
    tggtgATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAAt                                                                      <  1:1978889/68‑1 (MQ=255)
    tggtgATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAAc                                                                   <  1:1218762/71‑1 (MQ=255)
      gtgATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACg                                                                  >  1:4313002/1‑70 (MQ=255)
      gtgATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                 <  1:620887/71‑1 (MQ=255)
        gATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATaa                                                                    >  1:3382983/1‑66 (MQ=255)
        gATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAg                                                               >  1:436807/1‑71 (MQ=255)
         aTGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAg                                                               >  1:3747713/1‑70 (MQ=255)
         aTGGAAGTGGATCCGCGCTATACCCAGCACATCAGTGCGCTGGAAAAAATGTTCGTTATTGATAACGaac                                                               <  1:3898653/70‑2 (MQ=255)
          tGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAAt                                                                      <  1:4708468/62‑1 (MQ=255)
          tGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGc                                                             <  1:3084773/71‑1 (MQ=255)
            gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCaa                                                           <  1:3462851/71‑1 (MQ=255)
            gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCa                                                            <  1:3350998/70‑1 (MQ=255)
              aGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAg                                                         >  1:4924729/1‑71 (MQ=255)
               gTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAg                                                         <  1:2373207/70‑1 (MQ=255)
               gTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGc                                                        <  1:4605936/71‑1 (MQ=255)
                  gATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGa                                                      <  1:2288222/70‑1 (MQ=255)
                  gATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGa                                                      <  1:4296879/70‑1 (MQ=255)
                  gATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGAt                                                     >  1:1251498/1‑71 (MQ=255)
                   aTCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATc                                                    >  1:2899428/1‑71 (MQ=255)
                    tCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATcc                                                   <  1:1224351/71‑1 (MQ=255)
                    tCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATcc                                                   >  1:36431/1‑71 (MQ=255)
                     ccGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                >  1:472006/1‑57 (MQ=255)
                     ccGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAACGTTCGTTATTGATAACGa                                                                 <  1:1105268/56‑1 (MQ=38)
                      cgcgCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATccc                                                  >  1:5317251/1‑70 (MQ=255)
                       gcgcTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGa                                                      <  1:4661339/65‑1 (MQ=255)
                           tataCCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                 >  1:3940139/1‑50 (MQ=255)
                           tataCCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTc                                            >  1:2321946/1‑71 (MQ=255)
                               cccAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTa                                         <  1:3376774/70‑1 (MQ=255)
                                ccAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTAt                                        >  1:3421701/1‑70 (MQ=255)
                                  aGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGa                                                      >  1:3337355/1‑54 (MQ=255)
                                   ggACATCAGTGCGCTGGAAAAAATGTTCGTTATCAAt                                                                      >  1:4811976/1‑37 (MQ=255)
                                    gACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                 <  1:914384/41‑1 (MQ=255)
                                    gACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                   >  1:17121/1‑71 (MQ=255)
                                    gACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                   >  1:2294852/1‑71 (MQ=255)
                                     aCATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                   >  1:1548358/1‑70 (MQ=255)
                                     aCATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCt                                  <  1:1450677/71‑1 (MQ=255)
                                      cATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                   >  1:3408861/1‑69 (MQ=255)
                                      cATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTa                                 <  1:2229944/71‑1 (MQ=255)
                                       aTCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTa                                 >  1:1006748/1‑70 (MQ=255)
                                        tCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCg                                           <  1:2788383/59‑1 (MQ=255)
                                         cAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTa                                         >  1:1803594/1‑60 (MQ=255)
                                          aGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATg                             <  1:3088014/71‑1 (MQ=255)
                                          aGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATg                             <  1:3034621/71‑1 (MQ=255)
                                           gTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGa                                                      >  1:4507051/1‑45 (MQ=255)
                                            tGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGc                           <  1:4209332/71‑1 (MQ=255)
                                              tccTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                         <  1:2949687/69‑1 (MQ=255)
                                              cgcTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTa                                 >  1:980291/1‑63 (MQ=255)
                                              cgcTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCa                          <  1:2351592/70‑1 (MQ=255)
                                               gcTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATgg                            >  1:4202595/1‑67 (MQ=255)
                                                cTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTaa                                <  1:4312687/62‑1 (MQ=255)
                                                cTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAAcc                       <  1:4652555/71‑1 (MQ=255)
                                                   gAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGc                    >  1:2474849/1‑71 (MQ=255)
                                                   gAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCGCTGTCGTATTTCGCTAAATGGCAACCgg                     <  1:1695828/69‑1 (MQ=255)
                                                    aaaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCg                   <  1:402436/71‑1 (MQ=255)
                                                    aaaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCCCTGTCGTATTTCGCTAAATGGCAACCGGCg                   <  1:3978463/71‑1 (MQ=255)
                                                     aaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTa                                 >  1:319561/1‑56 (MQ=255)
                                                     aaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCg                   >  1:2545589/1‑70 (MQ=255)
                                                     aaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGa                  <  1:1029931/71‑1 (MQ=255)
                                                       aaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCg                                           >  1:3168971/1‑44 (MQ=255)
                                                       aaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCg                                           >  1:3175013/1‑44 (MQ=255)
                                                        aaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATg               >  1:4781202/1‑71 (MQ=255)
                                                         aTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGc              >  1:1513641/1‑71 (MQ=255)
                                                          tGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCg                                    >  1:382184/1‑48 (MQ=255)
                                                          tGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                         >  1:437828/1‑59 (MQ=255)
                                                            ttCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGccc            >  1:2883438/1‑70 (MQ=255)
                                                             tCGTTATTGATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCa          >  1:4871153/1‑71 (MQ=255)
                                                             tCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGc              >  1:1830558/1‑67 (MQ=255)
                                                              cGTTATTGATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCTCTaaa                               >  1:424677/1‑49 (MQ=25)
                                                                 tATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTAttt                                      >  1:3078757/1‑39 (MQ=255)
                                                                 tATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTAt      >  1:1584351/1‑71 (MQ=255)
                                                                  aTCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATc     >  1:3527223/1‑71 (MQ=255)
                                                                  aTCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATc     >  1:4237422/1‑71 (MQ=255)
                                                                     aaTAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGt  >  1:4890648/1‑71 (MQ=255)
                                                                     |                                                                      
AAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGT  >  AP009048/2162761‑2162900

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: