Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 476,666 G→C 5.2% A319P (GCC→CCC)  tsaD → tRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947476,6660GC 5.2% 360.4 / 9.3 115A319P (GCC→CCC) tsaDtRNA N6‑adenosine(37)‑threonylcarbamoyltransferase complex transferase subunit TsaD
Reads supporting (aligned to +/- strand):  ref base G (51/58);  new base C (2/4);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 117 total aligned reads displayed.

TCGCGGGTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCC  >  NC_002947/476517‑476810
                                                                                                                                                      |                                                                                                                                                
tatcgGGTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTgg                                                                                                                                                  >  1:833676/4‑150 (MQ=255)
tCGCGGGTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTgg                                                                                                                                                  >  2:341124/1‑150 (MQ=255)
   cgGGTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCg                                                                                                                                               <  2:393802/150‑1 (MQ=255)
    gggTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCgg                                                                                                                                              <  1:1010648/150‑1 (MQ=255)
    gggTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCgg                                                                                                                                              <  2:463344/150‑1 (MQ=255)
         gCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGgcagc                                                                                                                                         >  1:668922/1‑150 (MQ=255)
         gCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGgcagc                                                                                                                                         >  1:668923/1‑150 (MQ=255)
         gCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGgcagc                                                                                                                                         <  1:1300147/150‑1 (MQ=255)
         gCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGgcagc                                                                                                                                         <  1:1580568/150‑1 (MQ=255)
          cGTGAGTGACAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCTTCTGATGCTACGCGAGTGTGACTCGTACGCAATTCTGTACAGCTACTGGTGCGATGCTTTCCTTCGCCGGTTGCCAGCGCCTGCTGCCCGGccagca                                                                                                                                        <  1:896546/150‑1 (MQ=255)
           gTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGCCCGGccagcag                                                                                                                                       >  2:896546/1‑150 (MQ=255)
              aGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCaa                                                                                                                                    <  2:833676/150‑1 (MQ=255)
                tGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGCCCGGccagcaggcgct                                                                                                                                  <  1:1337549/150‑6 (MQ=255)
                tGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGGACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATAGCCTACGCCGGGTGCCAGCGCAGGCTGGCCGGGCAGGAACaaaa                                                                                                                                  >  2:228565/1‑148 (MQ=255)
                ggCCAACAAGGCCCTGCGCGCGTTCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGa                                                                                                                                  <  1:21141/149‑1 (MQ=255)
                  ccAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGAcc                                                                                                                                <  1:773053/150‑1 (MQ=255)
                  ccAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGCCCGGccagcaggcgctgg                                                                                                                                <  1:1444440/150‑8 (MQ=255)
                    aacaaGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTg                                                                                                                              >  2:723578/1‑150 (MQ=255)
                     acaaGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTgg                                                                                                                             <  2:1643470/150‑1 (MQ=255)
                           cccTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCa                                                                                                                       <  1:288315/150‑1 (MQ=255)
                             cTGCGCGCGTACAAAAAGGACGTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGc                                                                                                                     >  2:1912945/1‑150 (MQ=255)
                                 gcgcgTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGc                                                                                                                 <  2:56628/150‑1 (MQ=255)
                                     gTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGc                                                                                                             <  2:560692/150‑1 (MQ=255)
                                         cTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGc                                                                                                                            >  1:1051590/1‑131 (MQ=255)
                                         cTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGc                                                                                                                            <  2:1051590/131‑1 (MQ=255)
                                          tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGCTGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCt                                                                                                                               >  2:449738/1‑127 (MQ=255)
                                          tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCt                                                                                                                               <  1:449738/127‑1 (MQ=255)
                                          tGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCt                                                                                                        <  2:1127787/150‑1 (MQ=255)
                                           ggaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTg                                                                                                       <  1:1082285/150‑1 (MQ=255)
                                            gaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGCCCGGccagcaggcgctg                                                                                                                                 >  1:1515706/1‑117 (MQ=255)
                                            gaggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGCCCGGccagcaggcgctg                                                                                                                                 <  2:1515706/123‑7 (MQ=255)
                                             aggaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGc                                                                                                     <  2:1501207/150‑1 (MQ=255)
                                               gaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTg                                                                                                                              <  1:1282184/123‑1 (MQ=255)
                                               gaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTg                                                                                                                              >  2:1282184/1‑123 (MQ=255)
                                               gaCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGACAGCGGCAGGGCCTGGCGATCAGCGTGCAGGCACGCTGGCCg                                                                                                   >  2:152434/1‑150 (MQ=255)
                                                aCATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGa                                                                                                  <  1:1360315/150‑1 (MQ=255)
                                                  aTGCTCTGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                <  2:1047421/150‑1 (MQ=255)
                                                  aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                >  2:1600048/1‑150 (MQ=255)
                                                  aTGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATg                                                                                                <  1:1069113/150‑1 (MQ=255)
                                                    gCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGa                                                                                              >  1:1900447/1‑150 (MQ=255)
                                                     cTCGGTAGCATCAAGTGCAACGTGCACTCCGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACCCCGGTTGCCAGCGCCTGCTGGCCGCGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTCGCCGAttgca                                                                                             >  2:811439/1‑146 (MQ=255)
                                                     cTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGaa                                                                                             >  1:769840/1‑150 (MQ=255)
                                                        ggTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAg                                                                                          <  1:1912945/150‑1 (MQ=255)
                                                             cATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgcc                                                                                     <  1:1582669/150‑1 (MQ=255)
                                                               tCAAGGGCACCGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccgc                                                                                   >  1:259125/1‑150 (MQ=255)
                                                               tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCATGCACGCTGGCCGATGGAACAGTTgccgc                                                                                   >  2:897081/1‑150 (MQ=255)
                                                               tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccgc                                                                                   >  1:946595/1‑150 (MQ=255)
                                                               tCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTgccgc                                                                                   >  2:1448224/1‑150 (MQ=255)
                                                                 aaGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGccgccg                                                                                 >  2:1675371/1‑150 (MQ=255)
                                                                  aGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCcgccgc                                                                                >  2:165213/1‑150 (MQ=255)
                                                                  aGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGACAGGTACGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCcgccgc                                                                                >  1:1572853/1‑150 (MQ=255)
                                                                   gggCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCt                                                                               >  1:1118254/1‑150 (MQ=255)
                                                                    ggCAACGTGTACTACGCGCGTCCGCAATTCTGTACTGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCg                                                                                                                                               <  2:674461/85‑1 (MQ=255)
                                                                    ggCAACGTGTACTACGCGCGTCCGCAATTCTGTACTGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCg                                                                                                                                               >  1:674461/1‑85 (MQ=255)
                                                                    ggCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACa                                                                                           <  1:1328081/137‑1 (MQ=255)
                                                                    ggCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACa                                                                                           >  2:1328081/1‑137 (MQ=255)
                                                                     gCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGc                                                                                                         <  2:230161/122‑1 (MQ=255)
                                                                     gCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGc                                                                                                         >  1:230161/1‑122 (MQ=255)
                                                                      cAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTa                                                                            >  1:745436/1‑150 (MQ=255)
                                                                       aaCGTGTCCTAAGCGCGTCCGCAAATCTGTACAAACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTaa                                                                           >  2:696190/1‑150 (MQ=255)
                                                                           tgtACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTa                                                                       <  2:1570948/150‑1 (MQ=255)
                                                                             tactacGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAgc                                                                     <  1:1675371/150‑1 (MQ=255)
                                                                                tacNCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGcgcg                                                                  >  2:1189007/1‑150 (MQ=255)
                                                                                   gcgcgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCg                                                               >  1:856731/1‑150 (MQ=255)
                                                                                       cgTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAg                                                            <  1:165213/150‑1 (MQ=255)
                                                                                        gTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGc                                                                                                                     <  2:1877861/92‑1 (MQ=255)
                                                                                        gTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGc                                                                                                                     >  1:1877861/1‑92 (MQ=255)
                                                                                          ccGCAATTCTGTACCGACAATGGTGCGCTGATAGCCACCGGCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGtt                                                         <  2:1590440/150‑1 (MQ=255)
                                                                                           cGCCATTCTGAACCGCGAAAGGTGGGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTc                                                        >  2:593953/1‑150 (MQ=255)
                                                                                           cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAg                                                                                          <  1:1924858/116‑1 (MQ=255)
                                                                                           cGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAg                                                                                          >  2:1924858/1‑116 (MQ=255)
                                                                                            gCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGaa                                                             >  2:1541494/1‑144 (MQ=255)
                                                                                               aTTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGaa                                                    <  1:341124/150‑1 (MQ=255)
                                                                                                ttCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTgg                                                                                                                             <  1:1387005/76‑1 (MQ=255)
                                                                                                ttCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTgg                                                                                                                             >  2:1387005/1‑76 (MQ=255)
                                                                                                 tCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGt                                                  <  2:1900447/150‑1 (MQ=255)
                                                                                                     tACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGAt                                                                                                 <  1:1647848/99‑1 (MQ=255)
                                                                                                     tACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGAt                                                                                                 >  2:1647848/1‑99 (MQ=255)
                                                                                                      aCCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCCGAAGTGCCGt                                             >  2:69617/1‑150 (MQ=255)
                                                                                                        cGACAATGGTGCGATGATCGGTGACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTc                                           <  1:1531858/150‑1 (MQ=255)
                                                                                                        cGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTc                                           >  1:650765/1‑150 (MQ=255)
                                                                                                          aCAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGa                                                                                                                          >  1:1543785/1‑69 (MQ=255)
                                                                                                                    gatgatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGccgccg                                                                                 <  1:108437/100‑1 (MQ=255)
                                                                                                                    gatgatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGccgccg                                                                                 >  2:108437/1‑100 (MQ=255)
                                                                                                                    gatgatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCaaa                               >  1:1398089/1‑150 (MQ=255)
                                                                                                                       gatCGCCTACGCCGGTTGCCAGCGCCTGTTGGCCGGCCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACaa                            <  2:448457/150‑1 (MQ=255)
                                                                                                                       gatCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTAGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACaa                            <  1:152434/150‑1 (MQ=255)
                                                                                                                             ctaaggcGGTTGAAAGCGACAGATGGACGGGGAGCAGCAAGACCTGGCGATCAGCGTGCAGGCAACATGGAAGATGGAAGAGTTGCAGCCACGGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACACGGcgcg                      >  2:363028/7‑150 (MQ=255)
                                                                                                                             cTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCTATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTcgcg                      <  1:63045/150‑1 (MQ=255)
                                                                                                                             cTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTcgcg                      >  2:343266/1‑150 (MQ=255)
                                                                                                                                     gTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGgttgt              >  1:821495/1‑150 (MQ=255)
                                                                                                                                      ttGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGgttgtt             <  2:473319/150‑1 (MQ=255)
                                                                                                                                       tGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGttgttg            <  1:696190/150‑1 (MQ=255)
                                                                                                                                        gCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGc           <  1:1448224/150‑1 (MQ=255)
                                                                                                                                         ccAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCg          >  2:1245312/1‑150 (MQ=255)
                                                                                                                                          cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCa                    <  2:1895138/139‑1 (MQ=255)
                                                                                                                                          cAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCa                    >  1:1895138/1‑139 (MQ=255)
                                                                                                                                              gcCTTCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTCTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACATGTCGCGCAGGTTGTTGcccttga     <  2:946595/150‑7 (MQ=255)
                                                                                                                                                ctgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGGAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGgcgc   <  2:856731/150‑1 (MQ=255)
                                                                                                                                                 tgctgGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTTCAAACAAGTCGCGCAGGTTGTTGCCGCCGCGcc  >  2:1349923/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TCGCGGGTGGCGTGAGTGCCAACAAGGCCCTGCGCGCGTCCCTGGAGGACATGCTCGGTAGCATCAAGGGCAACGTGTACTACGCGCGTCCGCAATTCTGTACCGACAATGGTGCGATGATCGCCTACGCCGGTTGCCAGCGCCTGCTGGCCGGGCAGCAGCAAGACCTGGCGATCAGCGTGCAGGCACGCTGGCCGATGGAACAGTTGCCGCCGCTGTAAGGTAGCGCGTCGAAGGTTCAGAAGTGCCGTTCGCGCCCTGCAAACAAGTCGCGCAGGTTGTTGCGGTGGCGCC  >  NC_002947/476517‑476810

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