Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,949,413 A→G 5.4% L186P (CTA→CCA)  yifB ← magnesium chelatase family protein and putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,949,4130AG 5.4% 226.0 / 3.2 74L186P (CTA→CCA) yifBmagnesium chelatase family protein and putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base G (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACTTCATCTTCGTTATCTTTCGCGACGA  >  NC_000913/3949267‑3949560
                                                                                                                                                  |                                                                                                                                                   
gCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTa                                                                                                                                                     <  4:341670/147‑1 (MQ=255)
      ggCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCg                                                                                                                                              >  8:362681/1‑148 (MQ=255)
        ccATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCgg                                                                                                                                             >  4:63486/1‑147 (MQ=255)
        ccATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCgg                                                                                                                                             >  2:7236/1‑147 (MQ=255)
            tAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTGGCGCCCGGGATAc                                                                                                                                        >  7:38667/1‑148 (MQ=255)
            tAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATAc                                                                                                                                        >  1:399533/1‑148 (MQ=255)
             aaTGCGGCGGGCGCGCATTGTTTTCCCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATCAAATCACTGAGATCATGTTGTCGCGCCCGGGAttc                                                                                                                                        <  3:128069/147‑3 (MQ=255)
             aaTACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATAc                                                                                                                                        <  1:7236/147‑1 (MQ=255)
             aaTACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACt                                                                                                                                       >  2:301278/1‑148 (MQ=255)
              aTACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACt                                                                                                                                       <  1:161453/147‑1 (MQ=255)
                   gCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCAt                                                                                                                                   >  7:417787/1‑146 (MQ=255)
                      ggCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGtt                                                                                                                              >  4:262128/1‑148 (MQ=255)
                       gCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGt                                                                                                                               <  1:375192/146‑1 (MQ=255)
                          aGCATTGTTTTACCTGTTACCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCTGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGtt                                                                                                                          <  2:120946/148‑1 (MQ=255)
                           gCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGttt                                                                                                                         <  1:301278/148‑1 (MQ=255)
                            cATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTTAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCGCGGGATACTGCATCAGTTGGTTTc                                                                                                                        <  1:291421/148‑1 (MQ=255)
                             aTTGTTTTACCTGTTCCCGGCGGCCCATTCAGTAAAAGGTTGTGGGCGCCAGCAGCGGTAATTTCCAGTCGTCGCTTTCCTTGTTCCTGACCGATAACATCACTGTGATCATGTTGTAGCGCCCGGGATACTGCTTCAGTTGGTTTCg                                                                                                                       <  2:325239/148‑1 (MQ=255)
                             aTTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCg                                                                                                                       <  2:330778/148‑1 (MQ=255)
                               tGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATACCATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCggg                                                                                                                     >  6:160043/1‑148 (MQ=255)
                                 ttttACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGCCCGATAACATCACTGAGATCATGGTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCg                                                                                                                   >  2:69461/1‑148 (MQ=255)
                                  tttACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGGCCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGt                                                                                                                  >  3:211586/1‑148 (MQ=255)
                                    tACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTGGCGCCCGGGATACTGCATCAATTGGTTTCGGGCGtt                                                                                                                 <  8:38667/147‑1 (MQ=255)
                                        tGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCga                                                                                                              >  5:74240/1‑146 (MQ=255)
                                         gTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCCGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGagag                                                                                                           >  1:18454/1‑148 (MQ=255)
                                          ttCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGc                                                                                                          <  7:356566/148‑1 (MQ=255)
                                             ccGGCGGCCCAATCAGTAAAAGGGTGTGCCCGCCAGCAGCGGTAATTTCCAGGCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTg                                                                                                       <  5:300078/148‑1 (MQ=255)
                                               ggcggcCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCt                                                                                                     >  5:80692/1‑148 (MQ=255)
                                                      cAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTc                                                                                              >  2:186671/1‑148 (MQ=255)
                                                        aTCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCggg                                                                                                                                            >  6:273640/1‑100 (MQ=255)
                                                        aTCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCggg                                                                                                                                            <  5:273640/100‑1 (MQ=255)
                                                           aGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGaa                                                                                         <  6:36876/148‑1 (MQ=255)
                                                           aGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGaa                                                                                         <  1:186671/148‑1 (MQ=255)
                                                            gTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTGCGGGCGTTCGAGAGCGTGCTTACCTTCCAGaaa                                                                                        <  4:277399/148‑1 (MQ=255)
                                                               aaaGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACg                                                                                      >  7:23245/1‑147 (MQ=255)
                                                               aaaGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGc                                                                                     <  2:91998/148‑1 (MQ=255)
                                                                   gTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCacag                                                                                 >  3:366851/1‑148 (MQ=255)
                                                                    ttGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCg                                                                                                               <  5:8272/117‑1 (MQ=255)
                                                                    ttGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCg                                                                                                               >  6:8272/1‑117 (MQ=255)
                                                                       tgCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCacagaca                                                                              <  3:63486/147‑1 (MQ=255)
                                                                         cccGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCt                                                                           >  7:51956/1‑148 (MQ=255)
                                                                              cagcagCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCCGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAg                                                                      <  2:18454/148‑1 (MQ=255)
                                                                                   gcGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTTCAgatgat                                                                 >  1:187823/1‑148 (MQ=255)
                                                                                     ggTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCg                                                               <  3:149122/148‑1 (MQ=255)
                                                                                         aTTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTa                                                           <  2:399533/148‑1 (MQ=255)
                                                                                           ttCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTAt                                                          >  3:269490/1‑147 (MQ=255)
                                                                                           ttCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTAt                                                          >  8:25371/1‑147 (MQ=255)
                                                                                           ttCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACACACAGACAGCTTGCAGATGATCGGCTAt                                                          >  1:117541/1‑147 (MQ=255)
                                                                                            tCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATc                                                         <  7:258841/147‑1 (MQ=255)
                                                                                             ccAGTCCTCGCTTTCCTTGTTCCTGCCCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCa                                                        <  4:211586/147‑1 (MQ=255)
                                                                                             ccAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCa                                                        <  7:60990/147‑1 (MQ=255)
                                                                                                   cTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATcc                                                 <  4:269490/148‑1 (MQ=255)
                                                                                                   cTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATc                                                  <  1:167213/147‑1 (MQ=255)
                                                                                                    tCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTGGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCt                                                >  7:336651/1‑148 (MQ=255)
                                                                                                    tCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCt                                                >  2:337957/1‑148 (MQ=255)
                                                                                                     cGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCtt                                               >  1:275166/1‑148 (MQ=255)
                                                                                                       cTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCa                                             <  7:362681/148‑1 (MQ=255)
                                                                                                        tttCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCAc                                            >  4:123087/1‑148 (MQ=255)
                                                                                                         ttCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCAcc                                           >  8:47604/1‑148 (MQ=255)
                                                                                                         ttCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTAGCAGATGATCGGCTATCAGGCATCCTACAcc                                           >  3:101730/1‑148 (MQ=255)
                                                                                                              tGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGCCAGCTTGCAGATGATCGGCTATCAGGCATCCTTCCCCGTTaa                                      >  2:121106/1‑148 (MQ=255)
                                                                                                               gTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       <  2:93215/146‑1 (MQ=255)
                                                                                                               gTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       <  5:3339/146‑1 (MQ=255)
                                                                                                                ttCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAAt                                     <  7:25371/147‑1 (MQ=255)
                                                                                                                   cTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAg                                  <  5:160043/147‑1 (MQ=255)
                                                                                                                     gACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGAATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGCCATCCTTCACAGTTAATTAGccc                               >  3:45448/1‑148 (MQ=255)
                                                                                                                      aCCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGCAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCCTCCGTCACCGTTAATTAGccc                               >  6:77419/1‑147 (MQ=255)
                                                                                                                         gATAACATCACTGAGATCATGTTGTGGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACt                           <  8:336651/148‑1 (MQ=255)
                                                                                                                                 cACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAgatga                                                                  >  2:334550/1‑101 (MQ=255)
                                                                                                                                 cACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAgatga                                                                  <  1:334550/101‑1 (MQ=255)
                                                                                                                                 cACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACTTCATCTTc                   >  3:142585/1‑148 (MQ=255)
                                                                                                                                     gagaTCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       >  4:61869/1‑124 (MQ=255)
                                                                                                                                     gagaTCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       <  3:61869/124‑1 (MQ=255)
                                                                                                                                        aTCATGTTGTAGCGCCCGGGATACTTCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       <  4:117041/121‑1 (MQ=255)
                                                                                                                                        aTCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTa                                       >  3:117041/1‑121 (MQ=255)
                                                                                                                                                 tAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACTTCATCTTCGTTATCTTTCGcgacg   <  4:366851/148‑1 (MQ=255)
                                                                                                                                                  aGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACTTCATCTTCGTTATCTTTCGcgacga  >  1:171665/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCAACAGGCCATTAATACGGCTGGCGAGCATTGTTTTACCTGTTCCCGGCGGCCCAATCAGTAAAAGGTTGTGCCCGCCAGCAGCGGTAATTTCCAGTCCTCGCTTTCCTTGTTCCTGACCGATAACATCACTGAGATCATGTTGTAGCGCCCGGGATACTGCATCAGTTGGTTTCGGGCGTTCGAGAGCGTGCTTACCTTCCAGAAACGCACAGACAGCTTGCAGATGATCGGCTATCAGGCATCCTTCACCGTTAATTAGCCCCACTTCATCTTCGTTATCTTTCGCGACGA  >  NC_000913/3949267‑3949560

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