Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,075,469 C→T 5.4% intergenic (+45/‑309) 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,4690CT 5.4% 914.3 / 26.4 278intergenic (+45/‑309)aldB‑II/pqqD‑IIaldehyde dehydrogenase/coenzyme PQQ synthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (138/124);  new base T (5/10);  total (143/135)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 282 total aligned reads displayed.

AGTCCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCT  >  NC_002947/3075321‑3075618
                                                                                                                                                      |                                                                                                                                                     
  tCCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGc                                                                                                                                                    >  1:6334933/1‑150 (MQ=255)
  tCCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGc                                                                                                                                                    >  2:242979/1‑150 (MQ=255)
      gtgtTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGt                                                                                                                                                <  2:793581/150‑1 (MQ=255)
        gtTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTc                                                                                                                                              <  1:3644769/150‑1 (MQ=255)
            gCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGa                                                                                                                                          >  2:5312254/1‑150 (MQ=255)
             ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGaa                                                                                                                                         <  1:725937/150‑1 (MQ=255)
              cGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAc                                                                                                                                        <  2:2383792/150‑1 (MQ=255)
               gTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAAcc                                                                                                                                       <  1:672592/150‑1 (MQ=255)
                  agaCCCACCAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                    <  2:2352270/150‑1 (MQ=255)
                  agaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                    <  1:6605934/150‑1 (MQ=255)
                  agaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGt                                                                                                                                    <  2:6808659/150‑1 (MQ=255)
                     cccACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaa                                                                                                                                 <  2:3574094/150‑1 (MQ=255)
                      ccACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTgaag                                                                                                                                >  2:4515590/1‑150 (MQ=255)
                           aGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCt                                                                                                                           <  1:5721428/150‑1 (MQ=255)
                             atgatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCtgt                                                                                                                         <  2:4054611/150‑1 (MQ=255)
                               gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGTAGAACCTGTgg                                                                                                                       <  1:2713876/150‑1 (MQ=255)
                               gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                       <  1:2158755/150‑1 (MQ=255)
                               gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                       <  1:1566229/150‑1 (MQ=255)
                                 tgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTggga                                                                                                                     <  2:4179658/150‑2 (MQ=18)
                                    ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTgg                                                                                                                       <  2:2517898/145‑1 (MQ=255)
                                    ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                  >  2:920408/1‑150 (MQ=255)
                                    ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                  >  2:688125/1‑150 (MQ=255)
                                    ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcg                                                                                                                  >  2:2606277/1‑150 (MQ=255)
                                    ttGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGgagcg                                                                                                                  <  2:6429769/150‑1 (MQ=18)
                                     tGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgg                                                                                                                 <  2:3180793/150‑1 (MQ=255)
                                      gACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgagcggg                                                                                                                <  1:5582630/150‑1 (MQ=255)
                                       aCCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCCTTTTGCGGGCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGgagcgg                                                                                                                 <  2:2050325/150‑1 (MQ=14)
                                        ccATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGggcgggcg                                                                                                              >  1:6808968/1‑150 (MQ=255)
                                         cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                             <  2:5194559/150‑1 (MQ=255)
                                         cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                             >  2:5844272/1‑150 (MQ=255)
                                         cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGa                                                                                                             >  1:6988621/1‑150 (MQ=255)
                                          aTTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAg                                                                                                            >  2:1748270/1‑150 (MQ=255)
                                           ttATCAGCAGACCACGACCCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGc                                                                                                           <  1:5327111/150‑1 (MQ=255)
                                             aTCAGCAGACCACGAACCGGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGccc                                                                                                         <  2:4384366/150‑1 (MQ=255)
                                                agcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCa                                                                                                      >  1:7296495/1‑150 (MQ=255)
                                                agcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCa                                                                                                      >  2:7341372/1‑150 (MQ=255)
                                                 gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTggg                                                                                                                      >  2:6574311/1‑133 (MQ=255)
                                                 gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                     <  1:4669369/150‑1 (MQ=255)
                                                 gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                     >  1:4440128/1‑150 (MQ=255)
                                                 gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                     >  1:4382971/1‑150 (MQ=255)
                                                 gcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTGCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCAGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaa                                                                                                     >  1:112472/1‑150 (MQ=255)
                                                   agACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCaaaa                                                                                                   >  2:3815454/1‑150 (MQ=255)
                                                    gACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAg                                                                                                  <  1:5860347/150‑1 (MQ=255)
                                                     aCCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAgg                                                                                                 <  2:7448382/150‑1 (MQ=255)
                                                       gaaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                               >  1:3366773/2‑150 (MQ=255)
                                                       cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGgggcgg                                                                                                                 <  2:4438240/132‑1 (MQ=255)
                                                       cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                               <  1:3712214/150‑1 (MQ=255)
                                                       cAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGc                                                                                               >  1:3528963/1‑150 (MQ=255)
                                                         aGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCg                                                                                             <  2:6997466/150‑1 (MQ=255)
                                                           aaCCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGc                                                                                           <  2:3318977/150‑1 (MQ=255)
                                                               tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                       >  1:4285501/1‑150 (MQ=255)
                                                               tgctgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcat                                                                                       >  2:2973597/1‑150 (MQ=255)
                                                                      gAGCTACGACACCAAACCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTCCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCAcatgcggtaa                                                                                <  2:7227931/150‑6 (MQ=255)
                                                                       aGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGGTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACCTGCGTGGcc                                                                               >  1:4126631/1‑150 (MQ=255)
                                                                             gACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCg                                                                         >  2:6107637/1‑150 (MQ=255)
                                                                                aTCAATCCGCGGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc                                                                      <  2:5760917/150‑1 (MQ=255)
                                                                                   aaTCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCaa                                                                   >  1:991941/1‑150 (MQ=255)
                                                                                     tCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACg                                                                 >  2:5580578/1‑150 (MQ=255)
                                                                                       cGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCa                                                               >  1:6907456/1‑150 (MQ=255)
                                                                                        gCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCgg                                                                                                                                     <  1:4460334/79‑1 (MQ=255)
                                                                                         cTGGGCTTCTTCTACACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAcc                                                             <  1:196147/150‑1 (MQ=255)
                                                                                          tGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAccc                                                            <  1:5939693/150‑1 (MQ=255)
                                                                                           gggCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTAGCCCGCTCCCACAGCTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGTGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCAcccc                                                           <  1:1079055/150‑1 (MQ=255)
                                                                                                tcttctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAg                                                      >  1:316480/1‑150 (MQ=255)
                                                                                                  ttctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                    >  2:1648373/1‑150 (MQ=255)
                                                                                                  ttctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCt                                                    <  1:4958697/150‑1 (MQ=255)
                                                                                                   tctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTc                                                   >  1:7375557/1‑150 (MQ=255)
                                                                                                    ctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                            >  2:1409396/1‑124 (MQ=255)
                                                                                                    ctAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTAc                                                                            <  1:1409396/124‑1 (MQ=255)
                                                                                                      aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                <  2:6268888/150‑1 (MQ=255)
                                                                                                      aaaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt                                                <  1:5120518/150‑1 (MQ=255)
                                                                                                       aaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTc                                               >  2:1349006/1‑150 (MQ=255)
                                                                                                       aaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTc                                               >  1:2890840/1‑150 (MQ=255)
                                                                                                       aaCCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTc                                               >  2:663346/1‑150 (MQ=255)
                                                                                                            gggagcgGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCGTGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg                                          >  2:5402982/5‑150 (MQ=255)
                                                                                                            cgTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg                                          <  1:1648373/150‑1 (MQ=255)
                                                                                                              tCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                        <  1:1702011/150‑1 (MQ=255)
                                                                                                              tCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTa                                        <  1:5763980/150‑1 (MQ=255)
                                                                                                                gcgGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCAGTGGGGGCGGGCGAGCCCGCAAAAGGGCCCGCAAATGCAAGGCCTACCCGTACCAAAGCACCCCCCAAGCTCGATCCAAAGTAGc                                      >  1:6300119/1‑150 (MQ=255)
                                                                                                                        ttCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCgg                                >  1:509049/1‑150 (MQ=255)
                                                                                                                          cgcgGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                              <  2:1877926/150‑1 (MQ=255)
                                                                                                                             gggCTCGCCCGCCCCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCAAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAg                           >  1:3050598/1‑150 (MQ=17)
                                                                                                                               gCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTc                         >  1:4837959/1‑150 (MQ=255)
                                                                                                                               gCTTGCCCGCGCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTc                         <  1:242979/150‑1 (MQ=255)
                                                                                                                                cTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTcc                        <  2:268012/150‑1 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTc                              <  2:5158944/143‑1 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                       >  1:6572763/1‑150 (MQ=255)
                                                                                                                                 tcgcccgcccCCACAGGTTCTTCACCGGTTCTGGAACCAGCGGTGAACCTGTGGGAGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCg                       >  1:854503/3‑150 (MQ=18)
                                                                                                                                    cccGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccc                    >  1:2273462/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGccccc                  <  2:6642471/150‑1 (MQ=255)
                                                                                                                                          cccACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGc              <  2:6907456/150‑1 (MQ=255)
                                                                                                                                            cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt            >  1:2584895/1‑150 (MQ=255)
                                                                                                                                            cacaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAt            >  1:4792314/1‑150 (MQ=255)
                                                                                                                                             acaGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCAtt           <  2:2180428/150‑1 (MQ=255)
                                                                                                                                              caGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTa          >  1:390461/1‑150 (MQ=255)
                                                                                                                                                 gTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATg       >  2:5042978/1‑150 (MQ=255)
                                                                                                                                                  ttCACCGCTGGGTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGc      <  2:7306715/150‑1 (MQ=255)
                                                                                                                                                    cACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCtt    <  2:4254091/150‑1 (MQ=255)
                                                                                                                                                      ccGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCt  >  1:6680270/1‑150 (MQ=255)
                                                                                                                                                      ccGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCt  >  1:1359269/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGTCCGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCC‑‑TTCGCGGGCTTGCCCGCTCCCACAGGTTCACCGCTGGTTCCAGAACCGGTGAAGAACCTGTGGGGGCGGGCGAGCCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGTCCAGTCCGCCCCCGTGCATTAATGCTTCT  >  NC_002947/3075321‑3075618

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: