Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,075,473 T→C 5.9% intergenic (+49/‑305) aldB‑II → / → pqqD‑II aldehyde dehydrogenase/coenzyme PQQ synthesis protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,075,4730TC 5.9% 869.7 / 26.6 269intergenic (+49/‑305)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (135/118);  new base C (6/10);  total (141/128)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 273 total aligned reads displayed.

CCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCG  >  NC_002947/3075324‑3075622
                                                                                                                                                       |                                                                                                                                                     
 cGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTg                                                                                                                                                      <  2:5214147/150‑1 (MQ=255)
  ggTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTgg                                                                                                                                                     >  1:1687975/1‑150 (MQ=255)
     gtTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTc                                                                                                                                                  <  1:1648036/150‑1 (MQ=255)
      ttGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTcc                                                                                                                                                 >  2:4392082/1‑150 (MQ=255)
        ggCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAg                                                                                                                                               >  1:4244488/1‑150 (MQ=255)
        ggCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAg                                                                                                                                               >  2:5013399/1‑150 (MQ=255)
         gCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGa                                                                                                                                              >  2:4205310/1‑150 (MQ=255)
          ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGaa                                                                                                                                             >  2:1896214/1‑150 (MQ=255)
            gTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAcc                                                                                                                                           >  1:3460466/1‑150 (MQ=255)
            gTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAcc                                                                                                                                           >  1:647846/1‑150 (MQ=255)
             tGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCg                                                                                                                                          <  1:1554249/150‑1 (MQ=255)
              gagaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCgg                                                                                                                                         >  1:4076933/1‑150 (MQ=255)
                 aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                      <  1:4249152/150‑1 (MQ=255)
                 aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                      >  1:2695667/1‑150 (MQ=255)
                  cccACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaa                                                                                                                                     <  1:1007616/150‑1 (MQ=255)
                       aaGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAAcc                                                                                                                                >  2:1123573/1‑150 (MQ=255)
                         gatgatGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtg                                                                                                                              >  1:3143136/1‑150 (MQ=255)
                           tgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtgtg                                                                                                                            <  1:4603033/150‑1 (MQ=255)
                           tgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtgtg                                                                                                                            <  1:3533311/150‑1 (MQ=255)
                            gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGtt                                                                                                                                                   <  1:931384/126‑1 (MQ=255)
                             atgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggg                                                                                                                          <  2:5152504/150‑1 (MQ=255)
                              tgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCCTTTTGCGGGCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTgg                                                                                                                           <  2:3137064/150‑1 (MQ=17)
                               gcgtGACCATTCTCAGCAGCCCACGACCCTGCTGGTGAGCTACGACAGCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGGTCCAGAACCGGTGAAGAACCTGTggggg                                                                                                                        <  1:3930410/147‑1 (MQ=21)
                                      cATTATCAGCAGACCAAGACCCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGgagcgggcaa                                                                                                                 <  2:3942533/150‑3 (MQ=17)
                                      cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                 <  1:5571469/150‑1 (MQ=255)
                                       aTTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAg                                                                                                                <  1:4485174/150‑1 (MQ=255)
                                        ttATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGc                                                                                                               >  2:582124/1‑150 (MQ=255)
                                          aTCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGccc                                                                                                             <  1:3166554/150‑1 (MQ=255)
                                          aTCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGccc                                                                                                             <  2:1843001/150‑1 (MQ=255)
                                           tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCg                                                                                                                                          >  2:1231363/1‑120 (MQ=255)
                                           tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCg                                                                                                                                          <  1:1231363/120‑1 (MQ=255)
                                             agcagACCCAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCa                                                                                                          <  1:4890023/150‑1 (MQ=255)
                                             agcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                      >  1:1397030/1‑122 (MQ=255)
                                             agcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                      <  2:1397030/122‑1 (MQ=255)
                                              gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGACCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                         <  1:211752/150‑1 (MQ=255)
                                              gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                         >  1:1330276/1‑150 (MQ=255)
                                              gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                         >  1:2821488/1‑150 (MQ=255)
                                                agACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaaaa                                                                                                       >  1:4672608/1‑150 (MQ=255)
                                                  aCCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAgg                                                                                                     <  2:5300358/150‑1 (MQ=255)
                                                    cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                                   >  2:1833918/1‑150 (MQ=255)
                                                     aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                                                  >  1:477190/1‑150 (MQ=255)
                                                     aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                                                  <  2:2028742/150‑1 (MQ=255)
                                                       gAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCgg                                                                                                >  1:1797434/1‑150 (MQ=255)
                                                       gAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCgg                                                                                                >  2:4145978/1‑150 (MQ=255)
                                                       gAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCgg                                                                                                >  2:1240328/1‑150 (MQ=255)
                                                            tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                           >  1:2733088/1‑150 (MQ=255)
                                                            tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                           >  2:5458172/1‑150 (MQ=255)
                                                                   gAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACCGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGc                                                                                    <  1:262122/150‑1 (MQ=255)
                                                                       tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                >  1:927882/1‑150 (MQ=255)
                                                                       tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                <  2:4672541/150‑1 (MQ=255)
                                                                       tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                <  1:5253548/150‑1 (MQ=255)
                                                                       tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                <  2:730471/150‑1 (MQ=255)
                                                                        aCGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAcc                                                                               >  1:611909/1‑150 (MQ=255)
                                                                                aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCaa                                                                       >  2:3020011/1‑150 (MQ=255)
                                                                                aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCaa                                                                       >  2:121708/1‑150 (MQ=255)
                                                                                 aTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                >  1:3011231/1‑140 (MQ=255)
                                                                                 aTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAAc                                                                      <  2:3052675/150‑1 (MQ=255)
                                                                                  tCCGCTGGGCTTATTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAACGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACg                                                                     >  2:2694760/1‑150 (MQ=255)
                                                                                      cTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAcc                                                                 <  1:3817027/150‑1 (MQ=255)
                                                                                      cTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAcc                                                                 <  1:432434/150‑1 (MQ=255)
                                                                                           cttcttCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCa                                                            <  2:3622871/150‑1 (MQ=255)
                                                                                                 ctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCg                                                      >  1:757942/1‑150 (MQ=255)
                                                                                                   aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                    <  2:2054885/150‑1 (MQ=255)
                                                                                                   aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                    <  2:284907/150‑1 (MQ=255)
                                                                                                   aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                    <  2:391308/150‑1 (MQ=255)
                                                                                                   aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCGACCCGTACCAACGCACCCCGCAAGCTCGtt                                                    <  2:2838893/150‑1 (MQ=255)
                                                                                                       cgcgTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaa                                                <  2:3086857/150‑1 (MQ=255)
                                                                                                        gcgTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa                                               <  2:4280394/150‑1 (MQ=255)
                                                                                                         cgTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg                                              >  1:1117596/1‑150 (MQ=255)
                                                                                                           tCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                            <  2:4998061/150‑1 (MQ=255)
                                                                                                                gCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCAtt                                       <  2:1735056/150‑1 (MQ=255)
                                                                                                                gCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCAtt                                       <  1:484962/150‑1 (MQ=255)
                                                                                                                gCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCAtt                                       <  2:3303353/150‑1 (MQ=255)
                                                                                                                     ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGTACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAATGTAGCATTCgg                                    >  1:2853376/1‑150 (MQ=255)
                                                                                                                     ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                    >  1:2352875/1‑150 (MQ=255)
                                                                                                                     ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                    >  1:2352859/1‑150 (MQ=255)
                                                                                                                         cgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCa                                >  1:1023707/1‑150 (MQ=255)
                                                                                                                             ctggCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAACGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTcc                            <  2:3822158/147‑1 (MQ=255)
                                                                                                                              tcgcccgcccCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                           >  2:941025/3‑150 (MQ=18)
                                                                                                                                 cccGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccc                        >  2:4516593/1‑150 (MQ=255)
                                                                                                                                  ccGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGcccc                       >  1:3945502/1‑150 (MQ=255)
                                                                                                                                   cGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccccc                      >  1:2563706/1‑150 (MQ=255)
                                                                                                                                   cGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCGAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCa                                <  2:2982652/140‑1 (MQ=17)
                                                                                                                                    gCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCg                     >  1:3408356/1‑150 (MQ=255)
                                                                                                                                        ccACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                 >  2:3226183/1‑150 (MQ=25)
                                                                                                                                        ccACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                 >  2:4749327/1‑150 (MQ=25)
                                                                                                                                        ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                 >  1:5290159/1‑150 (MQ=255)
                                                                                                                                        ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                 >  1:2968349/1‑150 (MQ=255)
                                                                                                                                         cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                <  2:885388/150‑1 (MQ=255)
                                                                                                                                         cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                >  1:4988519/1‑150 (MQ=255)
                                                                                                                                             ggTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAAt            <  2:2254166/150‑1 (MQ=255)
                                                                                                                                             ggTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAAt            <  1:5379636/150‑1 (MQ=255)
                                                                                                                                               ttCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCACTAATGc          >  2:3507455/1‑150 (MQ=255)
                                                                                                                                               ttCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCATACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGc          >  1:3722347/1‑150 (MQ=255)
                                                                                                                                                  aCCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTCATGCTTc       >  1:5478267/1‑150 (MQ=255)
                                                                                                                                                  aCCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTc       >  1:2436034/1‑150 (MQ=255)
                                                                                                                                                    cGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTa     <  1:4162814/150‑1 (MQ=255)
                                                                                                                                                     gCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCACAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTAc    >  2:3357408/1‑150 (MQ=255)
                                                                                                                                                      cTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTAcc   <  2:477190/150‑1 (MQ=255)
                                                                                                                                                       tGGTTCCAGAACCGGTGAAGAAGCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCg  >  2:4060414/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCG  >  NC_002947/3075324‑3075622

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