Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084636,352,9930AT23.8% 891.8 / 72.2 378intergenic (+73/‑46)dgcB/PA14_RS29035dimethylglycine demethylation protein DgcB/electron transfer flavoprotein subunit alpha/FixB family protein
Reads supporting (aligned to +/- strand):  ref base A (213/66);  new base T (0/87);  total (213/165)
Fisher's exact test for biased strand distribution p-value = 2.40e-42
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 409 total aligned reads displayed.

CGGAAGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGGCGGCGAGGTGC  >  NC_008463/6352844‑6353141
                                                                                                                                                     |                                                                                                                                                    
cGGAAGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCaggga                                                                                                                                                      <  2:9506751/150‑1 (MQ=255)
 ggAAGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCagggag                                                                                                                                                     >  1:8324593/1‑150 (MQ=255)
   aaGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCagggaggg                                                                                                                                                   <  1:685648/150‑1 (MQ=255)
   aaGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCagggaggg                                                                                                                                                   <  2:8335730/150‑1 (MQ=255)
     gcgcTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAggggg                                                                                                                                                 <  2:5638509/150‑1 (MQ=255)
      cgcTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGc                                                                                                                                                <  1:943481/150‑1 (MQ=255)
      cgcTGAGAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGTGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGTCAGCACGGATAGCGGCATCTCCATCCAGGGAGGGGGc                                                                                                                                                <  2:2607263/150‑1 (MQ=255)
             agagCAGGCGTAAGCGCCGGCAAAATAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGa                                                                                                                                         <  2:9783555/150‑1 (MQ=255)
               agCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGAcc                                                                                                                                       >  2:4978189/1‑150 (MQ=255)
               agCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGAcc                                                                                                                                       >  2:8119728/1‑150 (MQ=255)
                   ggCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGt                                                                                                                                   >  2:1985800/1‑150 (MQ=255)
                       gAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGACTCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGcc                                                                                                                               >  2:10670268/1‑150 (MQ=255)
                       gAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGcc                                                                                                                               >  2:4853630/1‑150 (MQ=255)
                           cgccgcAAAAAACGCGGGGGCGGCGTGACGCGGTGGGGGTGGTGCATTGATGTTTGTGTCTAGCTTCACAGCGATGCCGCTGCTGCCCGGGGCCGCCCGGATCGCGGCCCCCCCCCCCCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAg                                                                                                                           <  2:10290678/150‑1 (MQ=255)
                             cgggCAAAAACGCCGGTGGGGCGTGAAGCCGGGGCGGAGTGGCTTTGATGAGATTTTCGAGCTTACAATCTAATCCGATGCTGACCGGGGCAGCACGGATTGCGGCCTCTCCCTCCCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGcc                                                                                                                         <  2:7058204/148‑1 (MQ=255)
                              cGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCca                                                                                                                        >  1:7857979/1‑150 (MQ=255)
                              cGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCca                                                                                                                        >  1:3987010/1‑150 (MQ=255)
                               ggCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCcac                                                                                                                       >  1:10682928/1‑150 (MQ=255)
                                gCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCcaca                                                                                                                      >  1:8609851/1‑150 (MQ=255)
                                gCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCcaca                                                                                                                      >  2:5491430/1‑150 (MQ=255)
                                  aaaaaaGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAgg                                                                                                                    >  2:7700181/1‑150 (MQ=255)
                                  aaaaaaGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAgg                                                                                                                    >  2:4709713/1‑150 (MQ=255)
                                   aaaaaGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGa                                                                                                                   >  2:9098977/1‑150 (MQ=255)
                                        cccGGTGCGGCGTGAAGCCGTGGCGGTGGTGCTTTGATGAGAGTGTCGAGCTCCACAGCGAATCCGCTGCTGACCGGGGCAGCACGGTTAGCGGCCTCCCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCAt                                                                                                              <  2:1854064/149‑1 (MQ=255)
                                           ggggcggcgTGAAGCCGGGGCGGTGGTGCATTGATGAGAGTGTCGAGCTCCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCACCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTcc                                                                                                           <  2:7313296/147‑1 (MQ=255)
                                           ggTGCGGCGTGAAGCCGTGGCGGCGGGGCATTGATGAGAGGGACGCGCTACCCATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCACTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTcc                                                                                                           <  2:9132160/150‑1 (MQ=255)
                                                 gcgTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACCTCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCAGTCCGCATGa                                                                                                     >  1:7106116/1‑150 (MQ=255)
                                                     gAAGCGGTGGCGGAGGTGCTTTGATTCGAGTGTCTAGCTTCACATCGCATCCGCTGGTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGa                                                                                                 <  2:9281976/150‑1 (MQ=255)
                                                         ccGTGGGGGAGGTGCATTGATGAGTGTGGCGCGGGACACATCGAATCCGGTGCTGCCCGGGGCAGCACGGATTGCGCCCTCGCCCCCCTGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATc                                                                                             <  1:10851062/150‑1 (MQ=255)
                                                         ccGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATc                                                                                             >  2:7543966/1‑150 (MQ=255)
                                                            tGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtc                                                                                          >  1:3957722/1‑150 (MQ=255)
                                                            tGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACAGCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtc                                                                                          >  1:161067/1‑150 (MQ=255)
                                                             ggCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtcc                                                                                         >  1:9157632/1‑150 (MQ=255)
                                                              gCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccg                                                                                        >  2:2672064/1‑150 (MQ=255)
                                                              gCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccg                                                                                        >  1:6814153/1‑150 (MQ=255)
                                                              gCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccg                                                                                        >  2:8376964/1‑150 (MQ=255)
                                                              gCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccg                                                                                        >  2:2796325/1‑150 (MQ=255)
                                                              gCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccg                                                                                        >  2:4309704/1‑150 (MQ=255)
                                                               gggtggTGCATTTTTGAGAGGTTCGAGATTCACACCAAATCGGCTGCTGACGGGGCCAGAACGGATCGCGCCCTCTCCCTCCTGGTTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccgt                                                                                       <  1:11159832/146‑1 (MQ=255)
                                                                gttggTGTATTGTTGTGTTTGTCTTGCTTCACAGCTATACCGTTGCTTACCGGGGCTGCACGTTTAGCGGCCTCTCCCCCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccgtc                                                                                      <  2:920050/147‑1 (MQ=255)
                                                                ggAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATcgtccgtc                                                                                      >  1:2750617/1‑150 (MQ=255)
                                                                  cggTGCATTTATGAGAGTGTCGCGCTTCACATCTAATCCGCTGCTGACCGGGGCAGCACGGATCGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTcgc                                                                                    <  2:1824242/149‑1 (MQ=255)
                                                                   ggTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTcgcg                                                                                   >  1:55141/1‑150 (MQ=255)
                                                                    gTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGa                                                                                  >  1:9979592/1‑150 (MQ=255)
                                                                      gCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATc                                                                                >  1:7587854/1‑150 (MQ=255)
                                                                       cATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATcc                                                                               >  1:4400235/1‑150 (MQ=255)
                                                                        aTTGATGTGTGTTTCTAGCTTCAATTCGATTCCGCTGATGACCGGGGCAGATCGTTTAGGGGCCTCTCCCCACAGGGTGTGGGCGTCTGGACCGTGTAGCCGAATCCTCAGGAAGCATTACGCATGAACGACATCGTCCGTCGCGATCCa                                                                              <  2:6748525/150‑1 (MQ=255)
                                                                         ttGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCAc                                                                             >  1:6344534/1‑150 (MQ=255)
                                                                          tgatgtGAGTGTCGAGCTTCAACTCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCAcg                                                                            <  2:6854080/150‑1 (MQ=255)
                                                                            atgaGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCAcgcg                                                                          >  2:7983832/1‑150 (MQ=255)
                                                                            atgaGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCAcgcg                                                                          >  2:9755562/1‑150 (MQ=255)
                                                                               agagTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCg                                                                       >  2:10760802/1‑150 (MQ=255)
                                                                               agagTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCg                                                                       >  1:5193767/1‑150 (MQ=255)
                                                                                  gtgtCGCGCGTCCCATCGACTCCGCTGCTGCCCGGGGCAGCAGGGATTGCGGCCAGTCCCCCCAGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAAt                                                                    <  2:10969163/150‑1 (MQ=255)
                                                                                  gtgtCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGTATTGCGGCCTCTCCCTCCTGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAAt                                                                    <  2:1535837/150‑1 (MQ=255)
                                                                                  gtgtCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAAt                                                                    >  2:6052698/1‑150 (MQ=255)
                                                                                   tgtCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATg                                                                   >  2:11071084/1‑150 (MQ=255)
                                                                                    gtCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATgg                                                                  >  1:10548785/1‑150 (MQ=255)
                                                                                    gtCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATgg                                                                  >  1:1368549/1‑150 (MQ=255)
                                                                                     gcgcgCTACACAGCGAACCCGCGGCTCCCGGGGGCCGCAGGGATAGCGGCCTCCCCCCCCTGGGGGGGGGCCCCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGa                                                                 <  1:9888737/146‑1 (MQ=255)
                                                                                       cccctccacaCCGAAGCCGCTCCTGACCGGGGCAGCCCGGCTTGCGCCCTCCCCCACCCGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATc                                                               <  2:3352174/144‑1 (MQ=255)
                                                                                         gCTTCACATCGAATCCGCGGCTGACCGGGGCAGCACGGATAGCGGCCACTCCCTCCAGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGc                                                             <  2:8229406/150‑1 (MQ=255)
                                                                                           ttCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGccc                                                           >  1:9607875/1‑150 (MQ=255)
                                                                                              acaTCTAATCCGCTGCTGACCGGGGCCGCACGTATCGCGGCCACCCCACCCTGGGTGGGGGCCCCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCa                                                        <  1:2384719/150‑1 (MQ=255)
                                                                                                   gAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGc                                                   >  2:9803610/1‑150 (MQ=255)
                                                                                                    aaTCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCCGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGcc                                                  >  2:10442535/1‑150 (MQ=255)
                                                                                                    aaTCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCCACCGcc                                                  >  2:3740201/1‑150 (MQ=255)
                                                                                                    aaTCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGcc                                                  >  2:8752030/1‑150 (MQ=255)
                                                                                                     aTCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCt                                                 >  1:7981539/1‑150 (MQ=255)
                                                                                                      tCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTg                                                >  1:2530420/1‑150 (MQ=255)
                                                                                                       ccGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGc                                               >  2:5246329/1‑150 (MQ=255)
                                                                                                         gctgctGACGGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCAt                                             <  2:8917877/150‑1 (MQ=255)
                                                                                                         gctgctGACCGGGGCAGCACGTATAGCGGCCACTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCAt                                             <  2:2466406/150‑1 (MQ=255)
                                                                                                         gctgctGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCAt                                             >  1:5080469/1‑150 (MQ=255)
                                                                                                         gctgctGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCAt                                             >  1:8522599/1‑150 (MQ=255)
                                                                                                          ctgctgACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCCGGGAGGGGGCCGCGGGACCGGGGAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATc                                            >  2:8878858/1‑150 (MQ=255)
                                                                                                             ctgACGGGGGCACCACGGATCGCGGCCTCTCCCACCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGc                                         <  2:10702398/150‑1 (MQ=255)
                                                                                                              tgACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCt                                        >  1:4119007/1‑150 (MQ=255)
                                                                                                               gACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTg                                       >  1:5982003/1‑150 (MQ=255)
                                                                                                                   ggggCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCCACCGCCTGCATCCGCTGCATg                                   >  2:1275435/1‑150 (MQ=255)
                                                                                                                   ggggCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATg                                   >  2:2971898/1‑150 (MQ=255)
                                                                                                                      gcccccgggATCGCGGCCTCCCCCCCCCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCgg                                <  1:10930704/143‑1 (MQ=255)
                                                                                                                       cagcaCGTATAGCGGCCTCACCCTCCCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGc                               <  1:10054947/150‑1 (MQ=255)
                                                                                                                        agcaCGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGcc                              >  1:4179889/1‑150 (MQ=255)
                                                                                                                        agcaCGGATAGCGGCCACTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGCTCGCCCGCAACCGCCTGCATCCGCTGCATGCGGcc                              >  1:6625976/1‑150 (MQ=255)
                                                                                                                            cgtaTAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGc                          <  1:10986767/147‑1 (MQ=255)
                                                                                                                               aTAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACt                       >  2:6530013/1‑150 (MQ=255)
                                                                                                                                ttgcgcCCTCCCCCCCCTGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTc                      <  1:2879053/148‑1 (MQ=255)
                                                                                                                                tcgCGGCCCCCCCCACCGGGGTGGGGCCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTc                      <  1:9687524/148‑1 (MQ=255)
                                                                                                                                 aGCGGCCTCTCCCTCCAGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTcc                     <  2:2256808/150‑1 (MQ=255)
                                                                                                                                  gCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCg                    >  1:865984/1‑150 (MQ=255)
                                                                                                                                         ctccctccCGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAgg             <  2:59140/150‑1 (MQ=255)
                                                                                                                                         ctccctccAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGAGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAgg             >  1:129941/1‑150 (MQ=255)
                                                                                                                                         ccccctccTGGGGGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAgg             <  1:1929508/148‑1 (MQ=255)
                                                                                                                                          tccctccAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAggc            >  1:8998463/1‑150 (MQ=255)
                                                                                                                                             ctccTGGGTGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAggcggc         <  1:3274823/150‑1 (MQ=255)
                                                                                                                                              tccAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGgcggcg        >  1:10220046/1‑150 (MQ=255)
                                                                                                                                                 cgggagggGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGGCGGCGAgg     <  2:6298250/149‑1 (MQ=255)
                                                                                                                                                  gggtgggGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGGCGGCGAGGt    <  1:6704152/150‑1 (MQ=255)
                                                                                                                                                  gggagggGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGGCGGCGAGGt    >  1:1763513/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGAAGCGCTGATAGAGCAGGCGGAAGCGCCGGCAAAAAAGCCGGTGCGGCGTGAAGCCGTGGCGGAGGTGCATTGATGAGAGTGTCGAGCTTCACATCGAATCCGCTGCTGACCGGGGCAGCACGGATAGCGGCCTCTCCCTCCAGGGAGGGGGCCGCGGGACCGGGTAGCCGCAGCCACAGGAAGCATTCCGCATGAACGACATCGTCCGTCGCGATCCACGCGCCGAATGGATCGCCCGCAACCGCCTGCATCCGCTGCATGCGGCCATGCACTCCGCCCAAGGCGGCGAGGTGC  >  NC_008463/6352844‑6353141

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