Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,662,436 G→A 9.6% P657L (CCT→CTT)  pbp3 ← penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,662,4360GA 9.6% 142.3 / 2.3 51P657L (CCT→CTT) pbp3penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (2/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACCATGCGCTACACAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGTACTACT  >  CP000730/1662299‑1662573
                                                                                                                                         |                                                                                                                                         
aCCATGCGCTACACAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACTTCTCTACCTAATTCTCCCCCTATTAACCATGTTGGTGGAACAGGc                                                                                                                                         >  2:412898/1‑140 (MQ=255)
aCCATGCGCTACACAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCCCCTGTTAAACATGGTGGTGGTCCAAGc                                                                                                                                         >  2:255863/1‑140 (MQ=255)
aCCATGCGCTACACAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAAGc                                                                                                                                         >  2:702602/1‑140 (MQ=255)
         tACACTATCGTCAGGTTAAATGTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGta                                                                                                                                <  2:470566/140‑1 (MQ=255)
          acacAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTGCCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAGGTCTCCACCTATTAAACACGTTGGAGGAACCGGGTGAttaagat                                                                                                                               >  2:558227/1‑135 (MQ=255)
          acacAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGtat                                                                                                                               >  1:425992/1‑140 (MQ=255)
              aaTCGTCATGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTCCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCTCCTGTTAACCATGGTGGTTGTACAGTCTGATTTGTATATAc                                                                                                                           <  2:203063/140‑1 (MQ=255)
                tCGTCAGGTTAAATTTTATTGGTCTTTGACTTTATTTTTATCATTTTTACCTAATTGCTTAAAGTAGTAGTTAATTACAGCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACaa                                                                                                                         >  1:491498/1‑140 (MQ=255)
                           aaTTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCt                                                                                                              >  1:581212/1‑140 (MQ=255)
                                tATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAAttt                                                                                                         <  2:402645/140‑1 (MQ=255)
                                 aTTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTATACCTAAGTCTCCACCTGTTAACCATGGTGATGGTACAGACTGCTTTGGAACTACAATTgcctgcgaaaatttt                                                                                                        >  2:549478/1‑126 (MQ=255)
                                       ctttgtctttATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGatca                                                                                                  <  1:509401/140‑1 (MQ=255)
                                        tttgtctttATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGatcat                                                                                                 <  1:478862/140‑1 (MQ=255)
                                            tctttATTTTTATCATTTTTACCTAACTGCATAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAAt                                                                                             <  1:214808/140‑1 (MQ=255)
                                              tttATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCCTCAATTg                                                                                           >  1:392123/1‑140 (MQ=255)
                                                tATTTTTATCATTTTTACCTAACTGCTTAAAGGAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCAACGGTTGgc                                                                                         <  2:277691/140‑1 (MQ=255)
                                                  tttttATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGgcgc                                                                                       >  2:759388/1‑140 (MQ=255)
                                                       aTCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATATGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATc                                                                                  >  1:248741/1‑140 (MQ=255)
                                                         cATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCt                                                                                <  1:719099/140‑1 (MQ=255)
                                                           tttttACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCtat                                                                              >  1:70745/1‑140 (MQ=255)
                                                           tttttACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTAACTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATACtat                                                                              >  2:169550/1‑140 (MQ=255)
                                                             tttACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCCCCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGTCGCGTATCCtatat                                                                            <  1:396982/140‑1 (MQ=255)
                                                             tttACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCtatat                                                                            <  1:655393/140‑1 (MQ=255)
                                                                       tGCTTAAAGTAGTAGTTAATTACCTCTCTACCTAAGTCTCCACCTTTTAACCATGGTGGTGGTACAGGCTGATTTGTATATAAACATGAAAACCCTAATCTTGGATCATACATTGCTGCGTATCCTTTATAAGTAGAGtt                                                                  <  2:444204/140‑1 (MQ=255)
                                                                              agtagtagtTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCtt                                                           >  1:196904/1‑140 (MQ=255)
                                                                              agtagtagtTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATAAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCtt                                                           >  1:503960/1‑140 (MQ=255)
                                                                                       tAATTACATCTCTACCTAAGTCTCCATCTTTTAACCATGGTGGGGGTACAAGCTGATTTGTATATACAATTGAAACCGCTAATTTTGGATCATCAATTGGCGCGTATCCTATCTAAGTAGAGTTAACTCTTGGCTCTCCg                                                  <  1:668651/140‑1 (MQ=255)
                                                                                       tAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCAGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCATATATAAGTAGAGTTACCTCTTGGCTCTCCg                                                  >  1:509653/1‑140 (MQ=255)
                                                                                       tAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAAGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCg                                                  <  1:702602/140‑1 (MQ=255)
                                                                                       tAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAAGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCg                                                  <  1:255863/140‑1 (MQ=255)
                                                                                        aaTTCCACCTCCACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGt                                                 <  2:173348/140‑1 (MQ=255)
                                                                                         aTTACAGCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTAAAGGCTGATTTGTAAATACAATTGAAAGCGCTAATTTTGGATCATCAATTGGCGCGGATCCTAACTAAGTAGAGTTAACTCTTAGCTCTCCggt                                                >  2:419375/1‑138 (MQ=255)
                                                                                               tctcTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTTAAAACGCTGATTTTGGATCCACAATTGGCGAGTATGCTATATAGGTAGAGTGAACTCTTGGCTCTCCGATTTggga                                          >  2:199712/1‑138 (MQ=255)
                                                                                                     ccTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAAAACTTc                                    >  1:446048/1‑140 (MQ=255)
                                                                                                        aaGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGc                                 <  1:412898/140‑1 (MQ=255)
                                                                                                          gTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCgg                               >  2:187572/1‑140 (MQ=255)
                                                                                                                 ccTGTTAACCAGGGTGGTGGTACAGGCTGATTATTTTATTGAATTTAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGt                        <  1:586200/140‑1 (MQ=255)
                                                                                                                    gTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCATATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGtttt                     >  1:208623/1‑140 (MQ=255)
                                                                                                                          cATGGTGGGGGTACAGGCTGATTTGTTTATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACcagc               <  2:36058/140‑1 (MQ=255)
                                                                                                                          aaTGCTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATTATCAATTGGCGCGTATCCTATATAAGTAGGGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACcagc               <  1:355175/139‑1 (MQ=255)
                                                                                                                           aTGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCATATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGACCACTTCTGCGGTACCCGTTTTACcagca              >  2:562950/1‑140 (MQ=255)
                                                                                                                             ggtggtggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTTTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGt            <  1:401932/140‑1 (MQ=255)
                                                                                                                             ggtggtggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGt            <  1:83575/140‑1 (MQ=255)
                                                                                                                             ggtggtggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGt            <  1:647604/140‑1 (MQ=255)
                                                                                                                             ggtggtggtACAGGCTGATTTGTAGATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGt            <  1:731239/140‑1 (MQ=255)
                                                                                                                                 gtggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGt        >  1:682434/1‑140 (MQ=255)
                                                                                                                                  tggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGta       <  2:70745/140‑1 (MQ=255)
                                                                                                                                  tggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATAATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACAAGCAGTAGGta       >  2:48210/1‑140 (MQ=255)
                                                                                                                                   ggtACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGtac      >  1:446982/1‑140 (MQ=255)
                                                                                                                                      aCAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGtactac   >  2:564227/1‑140 (MQ=255)
                                                                                                                                       cAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGTactact  >  1:689728/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
ACCATGCGCTACACAATCGTCAGGTTAAATTTTATTTGTCTTTGTCTTTATTTTTATCATTTTTACCTAACTGCTTAAAGTAGTAGTTAATTACATCTCTACCTAAGTCTCCACCTGTTAACCATGGTGGTGGTACAGGCTGATTTGTATATACAATTGAAAACGCTAATTTTGGATCATCAATTGGCGCGTATCCTATATAAGTAGAGTTAACTCTTGGCTCTCCGTTTTGGAACACTTCTGCGGTACCCGTTTTACCAGCAGTAGGTACTACT  >  CP000730/1662299‑1662573

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑