Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,162,829 C→T 5.3% I769I (ATC→ATT mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,162,8290CT 5.3% 249.1 / 2.8 75I769I (ATC→ATTmdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (41/30);  new base T (2/2);  total (43/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 = 8.85e-01

TAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATC  >  AP009048/2162760‑2162897
                                                                     |                                                                    
taAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATc                                                                      >  1:2569614/1‑70 (MQ=255)
 aaaGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCaa                                                                    >  1:1720863/1‑71 (MQ=255)
 aaaGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCaa                                                                    >  1:1714709/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:2294852/1‑71 (MQ=255)
                                     gACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                >  1:17121/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:3175013/1‑44 (MQ=255)
                                                        aaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCg                                        >  1:3168971/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)
                                                                     |                                                                    
TAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATC  >  AP009048/2162760‑2162897

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: