Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,378,540 C→T 7.5% G38S (GGT→AGT)  rplM ← 50S ribosomal subunit protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,378,5400CT 7.5% 162.4 / 4.3 53G38S (GGT→AGT) rplM50S ribosomal subunit protein L13
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (2/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-02

TTTGATACCACCGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCAATAAATAGTTACACGTTGGTGAA  >  NC_000913/3378397‑3378688
                                                                                                                                               |                                                                                                                                                    
tttGATACCACCGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCa                                                                                                                                                 <  2:16705/149‑1 (MQ=255)
tttGATACCACCGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCa                                                                                                                                                 >  1:245455/1‑149 (MQ=255)
      accaccGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACTGCGCAGGcgac                                                                                                                                           >  2:166749/1‑149 (MQ=255)
      accaccGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACTGCGCAGGcgac                                                                                                                                           >  2:347957/1‑149 (MQ=255)
        caccGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGcgacga                                                                                                                                         >  3:472225/1‑149 (MQ=255)
         accGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAg                                                                                                                                        >  2:347512/1‑149 (MQ=255)
           cGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGcc                                                                                                                                      <  1:184903/149‑1 (MQ=255)
                   ccGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGt                                                                                                                              <  1:347512/149‑1 (MQ=255)
                    cGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACTGCGCAGGCGACGAGCCAGTTCAGTa                                                                                                                             <  1:166749/149‑1 (MQ=255)
                    cGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACTGCGCAGGCGACGAGCCAGTTCAGTa                                                                                                                             <  1:347957/149‑1 (MQ=255)
                           gtgATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACGGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCTCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGAc                                                                                                                      <  1:420097/149‑1 (MQ=255)
                           gtgATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGAc                                                                                                                      >  2:92355/1‑149 (MQ=255)
                                  tACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCag                                                                                                               >  1:284562/1‑149 (MQ=255)
                                  tACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCag                                                                                                               >  2:320384/1‑149 (MQ=255)
                                   acacTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCaga                                                                                                              <  3:407021/149‑1 (MQ=255)
                                              gTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGt                                                                                                   >  2:221259/1‑149 (MQ=255)
                                                    ttgttgCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTc                                                                                             >  1:150766/1‑149 (MQ=255)
                                                       ttgCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTcaac                                                                                          <  4:283057/149‑1 (MQ=255)
                                                       ttgCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTcaac                                                                                          <  1:92355/149‑1 (MQ=255)
                                                                  aGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGt                                                                               >  4:129204/1‑149 (MQ=255)
                                                                     aaCTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTcgc                                                                            >  4:124593/1‑149 (MQ=255)
                                                                       cTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGt                                                                          >  1:107132/1‑149 (MQ=255)
                                                                              aGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACgg                                                                   <  2:150766/149‑1 (MQ=255)
                                                                              aGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACgg                                                                   >  2:409533/1‑149 (MQ=255)
                                                                               gCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTc                                                                                             <  1:462963/122‑1 (MQ=255)
                                                                               gCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTc                                                                                             >  2:462963/1‑122 (MQ=255)
                                                                                    cAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTg                                                             <  1:389738/149‑1 (MQ=255)
                                                                                     aGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTgg                                                            >  1:10880/1‑149 (MQ=255)
                                                                                       aaCGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGtt                                                          >  4:305308/1‑149 (MQ=255)
                                                                                          gatgatGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAg                                                       <  1:409533/149‑1 (MQ=255)
                                                                                          gatgatGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAg                                                       >  2:328714/1‑149 (MQ=255)
                                                                                           atgatgTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGc                                                      <  1:277444/149‑1 (MQ=255)
                                                                                            tgatgTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCt                                                     >  4:209560/1‑149 (MQ=255)
                                                                                                  aaTCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTaaaa                                               >  3:268321/1‑149 (MQ=255)
                                                                                                  aaTCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTaaaa                                               <  3:124593/149‑1 (MQ=255)
                                                                                                     cACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGtt                                            <  2:245455/149‑1 (MQ=255)
                                                                                                         ggTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCa                                        >  3:225337/1‑149 (MQ=255)
                                                                                                               tACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTaaaa                                  >  4:76617/1‑149 (MQ=255)
                                                                                                               tACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTaaaa                                  <  2:284562/149‑1 (MQ=255)
                                                                                                                aCGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTCACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAg                                 <  2:362913/149‑1 (MQ=255)
                                                                                                                    gCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTa                             <  1:328714/149‑1 (MQ=255)
                                                                                                                        aGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCa                         >  4:32374/1‑149 (MQ=255)
                                                                                                                            tATTCCGTTTTGCGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAAAATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCaataa                     >  4:402164/1‑149 (MQ=255)
                                                                                                                              ttCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTaaaa                                  >  2:240842/1‑134 (MQ=255)
                                                                                                                              ttCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTaaaa                                  <  1:240842/134‑1 (MQ=255)
                                                                                                                              ttCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCaataaat                   >  3:353143/1‑149 (MQ=255)
                                                                                                                              ttCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCaataaat                   >  2:320311/1‑149 (MQ=255)
                                                                                                                              ttCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCaataaat                   >  2:271733/1‑149 (MQ=255)
                                                                                                                               tCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCaataaata                  <  4:472225/149‑1 (MQ=255)
                                                                                                                                 cGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTAcc                           >  2:353735/1‑138 (MQ=255)
                                                                                                                                 cGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTAcc                           <  1:353735/138‑1 (MQ=255)
                                                                                                                                           ttACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCAATAAATAGTTACACGTTgg      >  1:78614/1‑149 (MQ=255)
                                                                                                                                             aCCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCAATAAATAGTTACACGTTGGTg    <  2:376037/149‑1 (MQ=255)
                                                                                                                                               cgcgCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCAATAAATAGTTACACGTTGGTGaa  <  2:78614/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TTTGATACCACCGATGTGGCCGGTGTGGTGATAGTACACTTTGTCAGTACGCTTGTTGCCGGTTACAGCAACTTTGTCAGCGTTCAGAACGATGATGTAATCACCGGTATCTACGTGCGGAGTGTATTCCGCTTTGTGCTTACCGCGCAGGCGACGAGCCAGTTCAGTAGCCAGACGGCCCAGAGTTTTACCGGTCGCGTCAACAACATACCAGTCGCGTTTTACGGTTTCTGGTTTAGCTGTAAAAGTTTTCATTAAAAGCTTACCCAATAAATAGTTACACGTTGGTGAA  >  NC_000913/3378397‑3378688

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