Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,075,488 T→C 6.3% intergenic (+64/‑290) 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,4880TC 6.3% 1053.3 / 36.8 330intergenic (+64/‑290)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (174/134);  new base C (12/9);  total (187/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 337 total aligned reads displayed.

AGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACA  >  NC_002947/3075339‑3075637
                                                                                                                                                     |                                                                                                                                                     
agaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                                       <  1:5259039/150‑1 (MQ=255)
agaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                                       <  2:6209476/150‑1 (MQ=255)
  aCCCACAAGATGATGCTTGACCATGATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                                     <  1:8359992/150‑1 (MQ=255)
   cccACCAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCGTCTTCTAAACCGCGTCGCGGCCTACGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTGCAGAACCGGTgaa                                                                                                                                                    <  1:536759/150‑1 (MQ=255)
     cacaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaaga                                                                                                                                                  >  2:7972762/1‑150 (MQ=255)
     cacaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaaga                                                                                                                                                  >  2:6755312/1‑150 (MQ=255)
      acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaagaa                                                                                                                                                 <  1:7861164/150‑1 (MQ=255)
         aGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCt                                                                                                                                              <  2:7666503/150‑1 (MQ=18)
          gatgatGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtg                                                                                                                                             >  1:5505321/1‑150 (MQ=255)
            tgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtgtg                                                                                                                                           <  2:8686413/150‑1 (MQ=255)
             gatgCTTGACCATTATCAGCAGACCACGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTgg                                                                                                                                          <  1:4572137/150‑1 (MQ=18)
             gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                                          <  1:5591259/150‑1 (MQ=255)
             gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                                          >  1:4809604/1‑150 (MQ=255)
                gCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggggg                                                                                                                                       >  2:3174713/1‑150 (MQ=255)
                  ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                                     <  2:2229130/150‑1 (MQ=255)
                  ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                                     >  2:794063/1‑150 (MQ=255)
                       cttgatcatcATACCTAGAACCTGCTTGTGATCTACGTCATCTATCCGCTGGGCTTCTTCTAACCCGCGTCGCGGCCTTCGCGGGCTTACCCGCTCCCCCAGGTTCACCGCTTGTTCCAGAACCGGTGAAGAACCTTTGGGGGCGGGCGa                                                                                                                                <  1:943061/146‑1 (MQ=18)
                       cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTACCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                                <  1:1839841/150‑1 (MQ=255)
                        aTTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAg                                                                                                                               <  2:493235/150‑1 (MQ=255)
                                 agACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaaaa                                                                                                                      <  2:230850/150‑1 (MQ=255)
                                  gACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtgtg                                                                                                                                           <  1:1106748/128‑1 (MQ=255)
                                     cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                                                  >  2:1356891/1‑150 (MQ=255)
                                     cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                                                  >  2:5125071/1‑150 (MQ=255)
                                      aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCACAAGGGcc                                                                                                                 >  1:6049382/1‑150 (MQ=255)
                                      aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                                                                 >  1:4154295/1‑150 (MQ=255)
                                      aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                                                                 >  2:4363028/1‑150 (MQ=255)
                                      aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGcc                                                                                                                 >  1:3262396/1‑150 (MQ=255)
                                         acccTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTGGCCTGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGc                                                                                                              <  1:5156742/148‑1 (MQ=255)
                                           ccTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGcac                                                                                                            >  1:2356217/1‑150 (MQ=255)
                                                tgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAACGCTGGTTCCAGAAACGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgca                                                                                                       >  1:1031950/1‑150 (MQ=255)
                                                       cTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTa                                                                                                <  1:4695850/150‑1 (MQ=255)
                                                       cTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATACATGGCCTa                                                                                                >  1:7326537/1‑150 (MQ=255)
                                                        tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGTGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                               >  1:2100385/1‑150 (MQ=255)
                                                             cATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTa                                                                                                <  2:6284320/144‑1 (MQ=255)
                                                             cATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTa                                                                                                >  1:6284320/1‑144 (MQ=255)
                                                              aTCAATCCTCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGTTTCCAGAACCGGCGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                                                                         >  2:8585363/1‑150 (MQ=255)
                                                              aTCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                                                                         >  2:463285/1‑150 (MQ=255)
                                                                  aTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAAc                                                                                     >  2:5845861/1‑150 (MQ=255)
                                                                  aTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAAc                                                                                     <  1:6949130/150‑1 (MQ=255)
                                                                   tCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACg                                                                                    >  2:7520606/1‑150 (MQ=255)
                                                                     cGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCa                                                                                  <  2:4142489/150‑1 (MQ=255)
                                                                        tGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCTGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAccc                                                                               <  2:7173783/150‑1 (MQ=255)
                                                                                ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                                       >  1:3716853/1‑150 (MQ=255)
                                                                                ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                                       >  1:2150352/1‑150 (MQ=255)
                                                                                ttctAAACCGCGCCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                                       >  1:1694389/1‑150 (MQ=255)
                                                                                 tctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTc                                                                      <  2:1812799/150‑1 (MQ=255)
                                                                                 tctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTc                                                                      >  2:1064904/1‑150 (MQ=255)
                                                                                   tAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGt                                                                    >  1:2762531/1‑150 (MQ=255)
                                                                                   tAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGt                                                                    >  2:5404973/1‑150 (MQ=255)
                                                                                       ccGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCa                                                                >  1:2553202/1‑150 (MQ=255)
                                                                                        cgcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaa                                                               >  2:5930388/1‑150 (MQ=255)
                                                                                         gcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa                                                              <  1:3094135/150‑1 (MQ=255)
                                                                                         gcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa                                                              >  1:444604/1‑150 (MQ=255)
                                                                                         gcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGGGACGAACCTGTGGGGGCGGGCGAGCCCGCACGAGGGCCGGCACATGCATGGCCTACCCGTACCACCGCACCCCGCAAGCTCGTGCCaaa                                                              >  2:6950891/1‑150 (MQ=255)
                                                                                             cgcccccTGCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCACACCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCACAAGGGCCGGCACATGCATGGCCTACCCGTACCACCGCACCCCGCAAGCGCGTTCCAAAGTAg                                                          <  2:5656143/145‑1 (MQ=255)
                                                                                              gcgGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGc                                                         >  2:2130476/1‑150 (MQ=255)
                                                                                                   cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCg                                                    >  1:3021841/1‑150 (MQ=255)
                                                                                                   cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCg                                                    <  1:907617/150‑1 (MQ=255)
                                                                                                   cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCg                                                    <  1:8123415/150‑1 (MQ=255)
                                                                                                    ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                                   <  1:1950277/150‑1 (MQ=255)
                                                                                                        cgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCa                                               <  2:3631589/150‑1 (MQ=255)
                                                                                                         gggCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAg                                              <  1:1248036/150‑1 (MQ=255)
                                                                                                           gCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTc                                            >  2:502749/1‑150 (MQ=18)
                                                                                                            ctcgcccgcccCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTcc                                           >  2:2820419/4‑150 (MQ=18)
                                                                                                            cTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTcc                                           <  1:2397426/150‑1 (MQ=255)
                                                                                                               gCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCTCCCCGCAAGCTCGGTCCAAAGTAGCATACGGTCCAGTCCGcc                                        >  2:7806152/1‑150 (MQ=255)
                                                                                                               gCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCCAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGcc                                        >  2:6869333/1‑150 (MQ=255)
                                                                                                                cccGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccc                                       <  1:3886663/150‑1 (MQ=255)
                                                                                                                 ccGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGcccc                                      >  2:2466173/1‑150 (MQ=18)
                                                                                                                 ccGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGcccc                                      >  2:5517988/1‑150 (MQ=18)
                                                                                                                  cGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccccc                                     <  2:2617872/150‑1 (MQ=255)
                                                                                                                    cTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGt                                   <  2:6559513/150‑1 (MQ=255)
                                                                                                                    cTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGt                                   <  1:2864884/150‑1 (MQ=255)
                                                                                                                       ccACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                                >  1:4298076/1‑150 (MQ=25)
                                                                                                                       ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                                >  1:7461835/1‑150 (MQ=255)
                                                                                                                        cacaGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                               >  2:3812070/1‑150 (MQ=25)
                                                                                                                        cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                               >  1:6644383/1‑150 (MQ=255)
                                                                                                                         acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                              >  2:6805216/1‑150 (MQ=255)
                                                                                                                         acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                              >  2:4712229/1‑150 (MQ=255)
                                                                                                                           aGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTaa                            <  2:3338961/150‑1 (MQ=255)
                                                                                                                             gTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATg                          >  2:6000463/1‑150 (MQ=255)
                                                                                                                             gTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATg                          >  2:3496733/1‑150 (MQ=255)
                                                                                                                             gTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATg                          >  2:8225068/1‑150 (MQ=255)
                                                                                                                                cACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCtt                       <  1:6432896/150‑1 (MQ=255)
                                                                                                                                  ccGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCt                     >  2:7160357/1‑150 (MQ=255)
                                                                                                                                     cTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTAcc                  >  1:1362940/1‑150 (MQ=255)
                                                                                                                                      tGGTTCCAGAACCGGTGAAGAAACTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCg                 <  2:1268637/150‑1 (MQ=255)
                                                                                                                                        gTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGa               >  1:5854129/1‑150 (MQ=255)
                                                                                                                                          tCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTACGCCCCCGTGCATTAATGCTTCTACCGGAAt             >  2:1857425/1‑150 (MQ=255)
                                                                                                                                           ccAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtc            >  1:221583/1‑150 (MQ=255)
                                                                                                                                            cAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcc           >  1:7112630/1‑150 (MQ=255)
                                                                                                                                             aGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGACtcct          >  1:7653635/1‑150 (MQ=255)
                                                                                                                                               aaCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcctcc        <  1:1433455/150‑1 (MQ=255)
                                                                                                                                               aaCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcctcc        >  1:6605502/1‑150 (MQ=255)
                                                                                                                                               aaCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcctcc        <  1:1995513/150‑1 (MQ=255)
                                                                                                                                                aCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCTGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCg       <  1:794063/150‑1 (MQ=255)
                                                                                                                                                 ccGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCgg      <  1:347674/150‑1 (MQ=255)
                                                                                                                                                    gTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGcac   >  2:144054/1‑150 (MQ=255)
                                                                                                                                                    gTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCCTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGAGCATTAATGCTTCTACCGGAATCCTCCGGcac   >  2:6985891/1‑150 (MQ=255)
                                                                                                                                                     tGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGcaca  >  1:4013018/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACA  >  NC_002947/3075339‑3075637

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