Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,149,912 G→A 5.8% intergenic (‑153/+104) PP_1008 ← / ← gapA RNA polymerase subunit sigma‑24/glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,149,9120GA 5.8% 298.2 / 6.9 103intergenic (‑153/+104)PP_1008/gapARNA polymerase subunit sigma‑24/glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (48/49);  new base A (4/2);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

Only 100 of 105 total aligned reads displayed.

GGGACCTGGCCGACGAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTT  >  NC_002947/1149763‑1150060
                                                                                                                                                     |                                                                                                                                                     
gggACCTGGCCGACGAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTgg                                                                                                                                                       <  2:450786/150‑1 (MQ=255)
 ggACCTGGCCGACGAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGgc                                                                                                                                                      <  1:1805026/150‑1 (MQ=255)
        gCCGACGAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATcc                                                                                                                                               <  1:394319/150‑1 (MQ=255)
           gacgaAAGGTTAATGATAATGCTTATCATCCATGCGGCAGGTCAACTGTGGAATTATCATGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTg                                                                                                                                           <  2:527822/150‑1 (MQ=255)
            acgaAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGTGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTgg                                                                                                                                          <  1:1163440/150‑1 (MQ=255)
            acgaAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTgg                                                                                                                                          <  1:2144108/150‑1 (MQ=255)
              gAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGG                                                                                                                                        <  1:1868893‑M1/150‑1 (MQ=255)
                   gTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCG                                                                                                                                   >  2:2062940‑M1/1‑150 (MQ=255)
                         gATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGGAGCGGGCtt                                                                                                                              >  2:1775033/1‑150 (MQ=255)
                          aTAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGc                                                                                                                            >  2:1814108/1‑150 (MQ=255)
                           tAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCCACCCCAGCTACAACCACGGCACAGGGCAGTCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGcc                                                                                                                           <  1:285580/150‑1 (MQ=12)
                                cTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGC                                                                                                                                    >  1:1540146‑M1/1‑136 (MQ=255)
                                cTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGC                                                                                                                                    <  2:1540146‑M1/136‑1 (MQ=255)
                                 ttATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGaa                                                                                                                     >  1:1620026/1‑150 (MQ=255)
                                 ttATCATCCGTGCGGCAAGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTTCCCGCGaa                                                                                                                     >  1:1179952/1‑150 (MQ=255)
                                      agccGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACgg                                                                                                                >  2:1610940/3‑150 (MQ=255)
                                       tcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCcgc                                                                                                                        <  2:810800/141‑1 (MQ=255)
                                       tcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCcgc                                                                                                                        >  1:810800/1‑141 (MQ=255)
                                       tcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACggg                                                                                                               <  1:661661/150‑1 (MQ=255)
                                        ccGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGc                                                                                                              <  1:1814108/150‑1 (MQ=255)
                                        ccGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGc                                                                                                              <  2:1439520/150‑1 (MQ=255)
                                          gTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATc                                                                                                                                                <  2:575975/115‑1 (MQ=255)
                                          gTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATc                                                                                                                                                >  1:575975/1‑115 (MQ=255)
                                           tGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAACAGATGGCGCGATCCCTGGGGGCGCCGGCTTGCCCGCGAACACGGGGGaa                                                                                                           <  2:94941/150‑1 (MQ=38)
                                            ccgggccaacgggcggatgctcaaggtgctggcctggtatgacaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGACGCGACCCGTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAg                                                                                                          <  1:958302‑M1/49‑1 (MQ=255)
                                                 aGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGgcttgcccgcgaagccaacaccgtggtgg                                                                                                     >  1:1993211‑M1/1‑121 (MQ=255)
                                                      aaGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGACCCGTGTGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggc                                                                                                >  2:2323386‑M1/1‑116 (MQ=255)
                                                        gtgtGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATcc                                                                                              <  1:2205402/150‑1 (MQ=255)
                                                              aaTTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTGTTGCTCAATTGCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCCAACACTCAGGACGCCGGTGCCATCCACCACg                                                                                        >  1:125910/1‑150 (MQ=25)
                                                               aTTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACgg                                                                                       <  1:2062940/150‑1 (MQ=255)
                                                                ttATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGgt                                                                                      >  2:1824518/1‑150 (MQ=255)
                                                                   tCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCtt                                                                                                                              >  1:1443143/1‑107 (MQ=255)
                                                                   tCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCtt                                                                                                                              <  2:1443143/107‑1 (MQ=255)
                                                                    cAGGAAAAGGTAAAGATTGTAGTCACTGGATGCAGTTCGCCTCATCAGCTAGCTACTGTTGTGCTCAATTTCTAACAGCGGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTg                                                                                  <  1:542076/150‑1 (MQ=17)
                                                                            ggTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGgcttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcg                                                                          <  2:1993211‑M1/150‑57 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                         >  1:351104‑M1/1‑93 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGtt                                                                                   >  1:2189223/1‑140 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGtt                                                                                   <  2:2189223/140‑1 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         >  2:2049961/1‑150 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         >  1:1094394/1‑150 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         <  1:1610940/150‑1 (MQ=255)
                                                                             gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGATCGCAACTTTGTTGATCTATTTTTAAAAGCTGGCGCGATCACTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATACACCACGGTATTGGCttggtgg                                                                         >  2:83977/1‑145 (MQ=18)
                                                                             ctggtatgacaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGACGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         >  1:2362992‑M1/69‑150 (MQ=255)
                                                                             ctggtatgacaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGACGCGACCCGTGTGGGAGCGGCCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         >  1:95482‑M1/69‑150 (MQ=255)
                                                                               aaaGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGc                                                                       >  2:1249938/1‑150 (MQ=255)
                                                                                aaGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCa                                                                      >  1:311329/1‑150 (MQ=255)
                                                                                   aTTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGACCCGTGTGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagc                                                                   <  1:2323386‑M1/150‑64 (MQ=255)
                                                                                    ttGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGc                                                                  >  2:1446565‑M2/1‑149 (MQ=255)
                                                                                    gacaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGc                                                                  >  1:135328‑M1/62‑150 (MQ=255)
                                                                                     tGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagccc                                                                 <  2:351104‑M1/150‑66 (MQ=255)
                                                                                     acaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGACGCGACCCGTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         <  1:1472301‑M1/82‑1 (MQ=255)
                                                                                     acaacgagtgggggttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGACGCGACCCGTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCgg                                                                         >  2:1472301‑M1/61‑142 (MQ=255)
                                                                                          tCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGccgctcc                                                            >  1:1824518‑M2/1‑143 (MQ=255)
                                                                                              tGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCcaca                                                        <  1:2049961/150‑1 (MQ=255)
                                                                                                 gtGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTTCTCAATTTCTGCGACCTGGCGCGATCCCTGGGGGAGCGGGGTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACgg                                                     <  1:264647/149‑1 (MQ=21)
                                                                                                  gttcagcaaccggatgctggacaactgcctggcgctgcaccgggctgGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACggg                                                    <  1:1190760‑M1/103‑1 (MQ=255)
                                                                                                     cGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcg                                                 >  1:809105/1‑150 (MQ=255)
                                                                                                      gTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgc                                                <  1:1991063/150‑1 (MQ=255)
                                                                                                        ttGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGt                                              <  1:1249938/150‑1 (MQ=255)
                                                                                                        ttGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGt                                              <  2:1007034/150‑1 (MQ=255)
                                                                                                        ttGCTTTTTTTGGGAGCAACTGTCTCGTGGCATGTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgcgg                                              >  1:1642902/1‑149 (MQ=21)
                                                                                                         tGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTc                                             >  2:1779702/1‑150 (MQ=255)
                                                                                                          gCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTca                                            >  1:735638/1‑150 (MQ=255)
                                                                                                          accggatgctggacaactgcctggcgctgcaccgggctgGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTca                                            >  2:1190760‑M1/40‑150 (MQ=255)
                                                                                                          accggatgctggacaactgcctggcgctgcaccgggctgGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTca                                            >  2:1535524‑M1/40‑150 (MQ=255)
                                                                                                            tttttttGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                          <  2:2262272/134‑1 (MQ=255)
                                                                                                            tttttttGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                          >  1:2262272/1‑134 (MQ=255)
                                                                                                            tttttttGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCAGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagc                                          >  2:597211/1‑150 (MQ=255)
                                                                                                              gatgctggacaactgcctggcgctgcaccgggctgGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagcca                                        >  2:1107008‑M1/36‑150 (MQ=255)
                                                                                                                     gAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGc                                                                                                              >  1:4466/1‑73 (MQ=39)
                                                                                                                     gAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGc                                                                                                              <  2:4466/73‑1 (MQ=39)
                                                                                                                       gCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGGCCGGTGc                               >  1:962297/1‑150 (MQ=255)
                                                                                                                        cAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCa                              <  2:311329/150‑1 (MQ=255)
                                                                                                                          aCTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                         <  2:1657525‑M1/105‑58 (MQ=255)
                                                                                                                          aCTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGccttgcccgcgaagccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                         >  1:1657525‑M1/1‑48 (MQ=255)
                                                                                                                           cTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGc                                                                  <  2:1149113‑M2/111‑2 (MQ=255)
                                                                                                                           cTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGc                                                                  >  1:1149113‑M2/1‑110 (MQ=255)
                                                                                                                            tGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGGCCGCGAACACGGGCGAAGCCCGTGCCATCTACCACGGGCGTGTCTTCGCGGGCAAGGCCGCTCCCCCACGGGTCGCGTCAGCCAGCCCGGTGCAgcgc                          >  2:1882815/1‑150 (MQ=21)
                                                                                                                             gTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGcc                         >  1:1631025/1‑150 (MQ=255)
                                                                                                                              tCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                          >  1:376894/1‑116 (MQ=255)
                                                                                                                              tCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                          <  2:376894/116‑1 (MQ=255)
                                                                                                                               cTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAg                       <  2:1094394/150‑1 (MQ=255)
                                                                                                                                  gCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCACCCCGGTGCAGCGCCAGGCa                    <  1:1399027/150‑1 (MQ=255)
                                                                                                                                    tCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGgtgt                                                                                    >  2:477949/1‑84 (MQ=255)
                                                                                                                                    tCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGgtgt                                                                                    <  1:477949/84‑1 (MQ=255)
                                                                                                                                    tCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAgc                            <  2:2273020/140‑1 (MQ=255)
                                                                                                                                    tCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAgc                            >  1:2273020/1‑140 (MQ=255)
                                                                                                                                     cAATTTCTAAAAGCTGGCGCGCTCCCTGGGGGAGCGGGCTCGACCGCAAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGtt                 <  1:519256/150‑1 (MQ=255)
                                                                                                                                     cAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGtt                 <  2:382782/150‑1 (MQ=255)
                                                                                                                                        tttCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagcca                                        <  1:897793/124‑1 (MQ=255)
                                                                                                                                        tttCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagcca                                        >  2:897793/1‑124 (MQ=255)
                                                                                                                                         ttCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTc                                                  <  1:772967/113‑1 (MQ=255)
                                                                                                                                         ttCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTc                                                  >  2:772967/1‑113 (MQ=255)
                                                                                                                                          tCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGTCAGTTGTCCa            <  2:532031/150‑1 (MQ=255)
                                                                                                                                          tCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCa            <  2:1620026/150‑1 (MQ=255)
                                                                                                                                          tCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCa            <  1:1417975/150‑1 (MQ=255)
                                                                                                                                          tCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCa            <  2:125910/150‑1 (MQ=255)
                                                                                                                                            tAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGGCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGc          >  1:740881/1‑150 (MQ=255)
                                                                                                                                                    GGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCGGTGCCATCCACCAGGGTGTTGTCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGtt  <  1:1113876‑M1/150‑1 (MQ=255)
                                                                                                                                                    GGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGtt  >  1:1645169‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGACCTGGCCGACGAAAGGTTAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTT  >  NC_002947/1149763‑1150060

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