Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,072,463 A→G 5.7% intergenic (+66/‑288) 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,4630AG 5.7% 449.0 / 9.5 140intergenic (+66/‑288)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base A (82/50);  new base G (3/5);  total (85/55)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 150 total aligned reads displayed.

ACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACAA  >  NC_002947/3072314‑3072611
                                                                                                                                                     |                                                                                                                                                    
aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTga                                                                                                                                                      >  1:4265055/1‑150 (MQ=255)
  ccACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaag                                                                                                                                                    <  1:1948601/150‑1 (MQ=255)
   cacaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaaga                                                                                                                                                   >  2:1460873/1‑150 (MQ=255)
    acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGACCCGGTgaagac                                                                                                                                                  >  2:4422034/1‑149 (MQ=255)
    acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaagaa                                                                                                                                                  >  2:783214/1‑150 (MQ=255)
     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAAc                                                                                                                                                 <  1:4470849/150‑1 (MQ=255)
     caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAAc                                                                                                                                                 <  2:4027939/150‑1 (MQ=255)
        gatgatGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtg                                                                                                                                              >  1:74937/1‑150 (MQ=255)
          tgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCtgtg                                                                                                                                            <  1:407889/150‑1 (MQ=18)
          tgatgCTTGACAATTATCAGCAGACCAAGCACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCtgtg                                                                                                                                            <  2:3848285/150‑1 (MQ=18)
           gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                                           >  1:4453507/1‑150 (MQ=255)
             tgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggtg                                                                                                                                         <  1:340494/150‑3 (MQ=255)
                ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                                      >  1:1654010/1‑150 (MQ=255)
                ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                                      >  1:5043949/1‑150 (MQ=255)
                ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                                      >  2:3806607/1‑150 (MQ=255)
                ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCAGTGGgggcg                                                                                                                                      >  1:2490645/1‑150 (MQ=255)
                 tGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgg                                                                                                                                     <  2:2740911/150‑1 (MQ=255)
                   aCCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgggc                                                                                                                                   <  1:1236059/150‑1 (MQ=255)
                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                                 <  2:4074700/150‑1 (MQ=255)
                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                                 >  1:812229/1‑150 (MQ=255)
                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                                                 <  2:220841/150‑1 (MQ=255)
                     cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCAGTGGGGGCGGGCGa                                                                                                                                 <  2:2490645/150‑1 (MQ=255)
                       ttATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGc                                                                                                                               <  1:4748434/150‑1 (MQ=255)
                        tATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGTACCTGTGGGGGCGGGCGAGcc                                                                                                                              >  2:4323483/1‑150 (MQ=255)
                        tATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTGCGGGCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGccc                                                                                                                             <  1:5341596/150‑1 (MQ=17)
                         aTCAGCAGACCACGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGccc                                                                                                                             <  1:1051661/150‑1 (MQ=255)
                         aTCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGccc                                                                                                                             >  1:2697839/1‑150 (MQ=255)
                          tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAACCCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCg                                                                                                                            <  2:5359744/150‑1 (MQ=255)
                             gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                                         >  1:2741039/1‑150 (MQ=255)
                             gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                                         >  2:5012285/1‑150 (MQ=255)
                              cagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaaa                                                                                                                        >  1:55805/1‑150 (MQ=255)
                              cagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaaa                                                                                                                        >  1:2289103/1‑150 (MQ=255)
                                gACCACGAACCTGCTGGTGAGCTACGACATCACTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAg                                                                                                                      <  2:2848164/150‑1 (MQ=255)
                                 aCCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAgg                                                                                                                     >  1:49293/1‑150 (MQ=255)
                                   cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                                                   >  2:1441611/1‑150 (MQ=255)
                                      gAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCgg                                                                                                                >  1:2047186/1‑150 (MQ=255)
                                       aaCCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGc                                                                                                               >  2:757094/1‑150 (MQ=255)
                                          ctgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGcaca                                                                                                            >  1:4546806/1‑150 (MQ=255)
                                          ctgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGcaca                                                                                                            >  1:4231364/1‑150 (MQ=255)
                                              tgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgca                                                                                                        >  2:765605/1‑150 (MQ=255)
                                                 tGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATgg                                                                                                     >  1:4551305/1‑150 (MQ=255)
                                                 tGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATgg                                                                                                     >  2:3958683/1‑150 (MQ=255)
                                                 tGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATgg                                                                                                     >  2:5233585/1‑150 (MQ=255)
                                                      tACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAACGGGCCGGCGCATGCATGGCCTAc                                                                                                >  2:2923786/1‑150 (MQ=255)
                                                      tACGACATCAATCCGCTGGGCTTCTGCGAAACCGCGTCGCGGCCTACGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACATGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                                                >  2:524036/1‑150 (MQ=255)
                                                           cATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATgg                                                                                                     <  1:4855396/140‑1 (MQ=255)
                                                           cATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTa                                                                                           <  1:765605/150‑1 (MQ=255)
                                                            aTCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                                                                          >  2:3580941/1‑150 (MQ=255)
                                                            aTCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                                                                          >  2:3580939/1‑150 (MQ=255)
                                                               aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCaa                                                                                       >  1:5411320/1‑150 (MQ=255)
                                                               aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCaa                                                                                       <  2:2741039/150‑1 (MQ=255)
                                                                   cGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCa                                                                                   >  1:2197129/1‑150 (MQ=255)
                                                                            tcttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCacg                                                                          >  2:4425458/1‑148 (MQ=255)
                                                                            tcttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAg                                                                          >  2:452985/1‑150 (MQ=255)
                                                                              ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCCTGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                                        >  1:938974/1‑150 (MQ=255)
                                                                              ttctAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                                        >  1:1482029/1‑150 (MQ=255)
                                                                                 tAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGt                                                                     >  1:3417166/1‑150 (MQ=255)
                                                                                      cgcgTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaa                                                                >  2:5131145/1‑150 (MQ=255)
                                                                                          tCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                                            >  2:3313141/1‑150 (MQ=255)
                                                                                          tCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                                            <  2:1302080/150‑1 (MQ=255)
                                                                                            gcgGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCCGTGAAGAACCGGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGc                                                          >  1:1062567/1‑150 (MQ=255)
                                                                                                 ttttgcgGGCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCCAGCTCGTTCCAAAGTAGCATTCg                                                     >  1:3807684/5‑150 (MQ=14)
                                                                                                 cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGCGCCCGCACAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCg                                                     >  1:2392426/1‑150 (MQ=255)
                                                                                                 cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAg                                                                                                                      <  1:3305283/85‑1 (MQ=255)
                                                                                                 cTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAg                                                                                                                      >  2:3305283/1‑85 (MQ=255)
                                                                                                  ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                                    <  1:1127245/150‑1 (MQ=255)
                                                                                                    cgcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                                                  <  1:3023815/150‑1 (MQ=255)
                                                                                                    cgcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                                                  <  1:5115669/150‑1 (MQ=255)
                                                                                                    cgcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                                                  <  1:517952/150‑1 (MQ=255)
                                                                                                     gcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGATCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATgg                                                                                                     <  1:768515/98‑1 (MQ=255)
                                                                                                           ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCCTGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTCGCATTCGGTCCAGTCCg                                           >  2:4243687/1‑150 (MQ=255)
                                                                                                           ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                                           >  1:210734/1‑150 (MQ=255)
                                                                                                            tGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGc                                          <  1:4929146/150‑1 (MQ=255)
                                                                                                                cGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccccc                                      >  2:5379792/1‑150 (MQ=255)
                                                                                                                 gCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCg                                     >  2:2248423/1‑150 (MQ=255)
                                                                                                                     ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgc                                                                                                         <  2:2037614/78‑1 (MQ=255)
                                                                                                                     ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgc                                                                                                         >  1:2037614/1‑78 (MQ=255)
                                                                                                                     ccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCa                                 >  1:4022390/1‑150 (MQ=255)
                                                                                                                       acaGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                               >  1:16925/1‑150 (MQ=21)
                                                                                                                       acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt                               >  2:377314/1‑150 (MQ=255)
                                                                                                                         aGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTaa                             >  2:1593022/1‑150 (MQ=255)
                                                                                                                         aGGTTCACCGCTGGTTCAAGAACCGGTGAAGAACCGGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGAACCAACGCACCCCGAAAGCTCGATCCAAAGTAGCATTCTGTCAAGTACGCCCCCGTGCATTaa                             >  1:1899787/1‑150 (MQ=255)
                                                                                                                           gTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATg                           >  1:3026447/1‑150 (MQ=255)
                                                                                                                             tCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCt                         >  1:679946/1‑150 (MQ=255)
                                                                                                                                 cGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTa                     <  1:5027739/150‑1 (MQ=255)
                                                                                                                                  gCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTAc                    <  2:543506/150‑1 (MQ=255)
                                                                                                                                    tGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCg                  <  1:3647964/150‑1 (MQ=255)
                                                                                                                                          cAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                                            >  1:2283395/1‑102 (MQ=255)
                                                                                                                                          cAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcc            >  2:588828/1‑150 (MQ=255)
                                                                                                                                          cAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcc            >  2:2313840/1‑150 (MQ=255)
                                                                                                                                            gAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcctc          <  2:4231364/150‑1 (MQ=255)
                                                                                                                                            gAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAAtcctc          <  2:4551305/150‑1 (MQ=255)
                                                                                                                                              aCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCg        <  2:4265055/150‑1 (MQ=255)
                                                                                                                                              aCCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCg        <  2:3529717/150‑1 (MQ=255)
                                                                                                                                               ccGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCgg       >  2:3432578/1‑150 (MQ=255)
                                                                                                                                               ccGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCAGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCgg       <  2:812229/150‑1 (MQ=255)
                                                                                                                                                 ggTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGca     <  1:632902/150‑1 (MQ=255)
                                                                                                                                                  gTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGcac    <  2:1246652/150‑1 (MQ=255)
                                                                                                                                                    gaagaaCCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACaa  >  1:3233209/1‑150 (MQ=255)
                                                                                                                                                    gaagaaCCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACaa  <  1:1460873/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCTACCGGAATCCTCCGGCACAA  >  NC_002947/3072314‑3072611

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: