Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,072,453 A→G 5.1% intergenic (+56/‑298) 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,072,4530AG 5.1% 589.9 / 12.4 175intergenic (+56/‑298)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base A (101/65);  new base G (6/3);  total (107/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 177 total aligned reads displayed.

GTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATC  >  NC_002947/3072309‑3072600
                                                                                                                                                |                                                                                                                                                   
gTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAcc                                                                                                                                                <  2:2510390/150‑1 (MQ=255)
   agaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                             <  1:1549280/150‑1 (MQ=255)
     aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                           <  1:3927922/150‑1 (MQ=255)
     aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                           <  1:4066789/150‑1 (MQ=255)
       ccaccaGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgcag                                                                                                                                         >  2:3492422/1‑150 (MQ=255)
       ccACACGATGATGCTTGACCATTATCAGCAGACCAAGCACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCGTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaag                                                                                                                                         <  2:2807530/150‑1 (MQ=255)
         acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaagaa                                                                                                                                       >  1:4019938/1‑150 (MQ=255)
         acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaagaa                                                                                                                                       >  1:4402907/1‑150 (MQ=255)
            aGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCt                                                                                                                                    >  2:2589061/1‑150 (MQ=255)
            aGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCt                                                                                                                                    >  2:1200392/1‑150 (MQ=255)
            aGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCt                                                                                                                                    >  2:1200391/1‑150 (MQ=255)
               tgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCtgtg                                                                                                                                 <  2:1861823/150‑1 (MQ=18)
                gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                                <  1:1814776/150‑1 (MQ=255)
                 atgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggg                                                                                                                               <  2:4295099/150‑1 (MQ=255)
                   gCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggggg                                                                                                                             >  1:3760947/1‑150 (MQ=255)
                    cTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggc                                                                                                                            <  2:3383931/150‑1 (MQ=255)
                     ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                           <  1:921308/150‑1 (MQ=255)
                     ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                           >  2:1641651/1‑150 (MQ=255)
                      tGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgg                                                                                                                          >  1:1460823/1‑150 (MQ=255)
                        aCCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgggc                                                                                                                        >  1:1736615/1‑150 (MQ=255)
                         ccATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTGCGGGCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGgagcgggcaa                                                                                                                      <  1:3744687/150‑3 (MQ=17)
                          cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCTCAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCGGTGGGggcgggcgt                                                                                                                      >  1:2199105/1‑149 (MQ=255)
                          cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                      >  2:2390480/1‑150 (MQ=255)
                          cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                      >  1:1015871/1‑150 (MQ=255)
                          cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                      <  1:2481002/150‑1 (MQ=255)
                          cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                      >  1:1316820/1‑150 (MQ=255)
                           aTTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAg                                                                                                                     >  2:3171006/1‑150 (MQ=255)
                            ttATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGc                                                                                                                    <  1:1307926/150‑1 (MQ=255)
                               tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCg                                                                                                                 >  1:852439/1‑150 (MQ=255)
                               tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCg                                                                                                                 <  1:1384214/150‑1 (MQ=255)
                                 agcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCa                                                                                                               >  1:1215949/1‑150 (MQ=255)
                                  gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                              >  2:1240794/1‑150 (MQ=255)
                                   cagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAACGCTGGATCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCACGCaaa                                                                                                             >  1:955015/1‑150 (MQ=255)
                                      aCCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAgg                                                                                                          <  1:809498/150‑1 (MQ=255)
                                       ccAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAggg                                                                                                         >  2:4425173/1‑150 (MQ=255)
                                                tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                                >  1:3632569/1‑150 (MQ=255)
                                                tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                                >  1:3632595/1‑150 (MQ=255)
                                                     gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgcatg                                                                                           >  1:328720/1‑150 (MQ=255)
                                                       gAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGc                                                                                         <  1:2865847/150‑1 (MQ=255)
                                                        aGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGcc                                                                                        <  2:3177722/150‑1 (MQ=255)
                                                          cTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTa                                                                                      >  1:199347/1‑150 (MQ=255)
                                                           tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                     >  1:2010814/1‑150 (MQ=255)
                                                            aCGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAcc                                                                                    <  1:1439811/150‑1 (MQ=255)
                                                               aCATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGt                                                                                 >  1:528561/1‑150 (MQ=255)
                                                               aCATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGt                                                                                 >  1:2784313/1‑150 (MQ=255)
                                                                   cAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCa                                                                             >  2:634290/1‑150 (MQ=255)
                                                                    aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCa                                                                                                                                                     >  2:1948374/1‑77 (MQ=255)
                                                                    aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCa                                                                                                                                                     <  1:1948374/77‑1 (MQ=255)
                                                                      tCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACg                                                                          <  2:1736615/150‑1 (MQ=255)
                                                                         gCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAc                                                                       >  1:4003976/1‑150 (MQ=255)
                                                                                ttcttcTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCaa                                                                <  1:3161312/150‑1 (MQ=255)
                                                                                ttcttcTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCaa                                                                >  2:4173685/1‑150 (MQ=255)
                                                                                ttcttcTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCaa                                                                <  1:2589061/150‑1 (MQ=255)
                                                                                ttcttcTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCGGGGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACAGGCATGGCCTACCCGTACCAACGCACCCCGCac                                                                >  2:2222222/1‑149 (MQ=255)
                                                                                   ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                             <  2:1413766/150‑1 (MQ=255)
                                                                                   ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                             <  1:413818/150‑1 (MQ=255)
                                                                                    tctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGc                                                                  <  2:3021439/144‑1 (MQ=255)
                                                                                    tctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGc                                                                  >  1:3021439/1‑144 (MQ=255)
                                                                                    tctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTc                                                            >  2:2701333/1‑150 (MQ=255)
                                                                                      tAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGt                                                          >  2:765602/1‑150 (MQ=255)
                                                                                       aaaCCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                         >  2:2606086/1‑150 (MQ=255)
                                                                                       aaaCCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACACGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                         >  2:4461953/1‑150 (MQ=255)
                                                                                        aaCCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTc                                                        >  1:2371621/1‑150 (MQ=255)
                                                                                          ccGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTTGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAAATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCa                                                      >  1:971518/1‑150 (MQ=255)
                                                                                           cgcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaa                                                     >  2:1126215/1‑150 (MQ=255)
                                                                                           cgcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaa                                                     <  1:4392608/150‑1 (MQ=255)
                                                                                                cgcgGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAg                                                <  1:4313538/150‑1 (MQ=255)
                                                                                                    gCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCAtt                                            <  1:4508143/150‑1 (MQ=255)
                                                                                                     ccTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTc                                           >  1:1323539/1‑150 (MQ=255)
                                                                                                       ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                         <  2:3158574/150‑1 (MQ=255)
                                                                                                        gagcgGGCAAGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCGGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCACGCTCGTTCCAAAGTAGCATTCGGt                                        >  2:2276030/3‑150 (MQ=255)
                                                                                                         cgcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTGCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                                       <  2:3991630/150‑1 (MQ=255)
                                                                                                           cgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCa                                     <  2:2037689/150‑1 (MQ=255)
                                                                                                           cgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCa                                     <  1:72795/150‑1 (MQ=255)
                                                                                                             ggCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGt                                   <  2:243275/150‑1 (MQ=255)
                                                                                                                ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                                >  1:3941894/1‑150 (MQ=255)
                                                                                                                ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                                <  1:2919762/150‑1 (MQ=255)
                                                                                                                ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                                >  1:3410063/1‑150 (MQ=255)
                                                                                                                ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                                >  1:1667126/1‑150 (MQ=255)
                                                                                                                    ccGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGcccc                            >  1:3661399/1‑150 (MQ=18)
                                                                                                                     cGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccccc                           <  2:2289067/150‑1 (MQ=255)
                                                                                                                         cccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGc                       <  2:1350016/150‑1 (MQ=255)
                                                                                                                          ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                                 <  1:3205266/123‑1 (MQ=255)
                                                                                                                          ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                                 >  2:3205266/1‑123 (MQ=255)
                                                                                                                           cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                     >  2:854385/1‑150 (MQ=255)
                                                                                                                           cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                     >  2:4295957/1‑150 (MQ=255)
                                                                                                                           cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                     >  1:2963530/1‑150 (MQ=255)
                                                                                                                           cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt                     >  1:2533684/1‑150 (MQ=255)
                                                                                                                            acaGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCGAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                    >  2:3128014/1‑150 (MQ=21)
                                                                                                                            acaGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCGAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                    >  1:2262481/1‑150 (MQ=21)
                                                                                                                            acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                    >  2:667582/1‑150 (MQ=255)
                                                                                                                            acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                    >  1:2243389/1‑150 (MQ=255)
                                                                                                                                 ttCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGc               >  1:1061955/1‑150 (MQ=255)
                                                                                                                                  tCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCt              >  1:2642143/1‑150 (MQ=255)
                                                                                                                                     ccGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCt           <  2:3781521/150‑1 (MQ=255)
                                                                                                                                       gCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTAc         >  1:3101745/1‑150 (MQ=255)
                                                                                                                                         tGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCg       <  1:798346/150‑1 (MQ=255)
                                                                                                                                           gTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGa     >  1:4311812/1‑150 (MQ=255)
                                                                                                                                              ccAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtc  >  2:1064563/1‑150 (MQ=255)
                                                                                                                                              ccAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACAAGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtc  >  2:68946/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATC  >  NC_002947/3072309‑3072600

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: