Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U00096935,5840AT20.9% 164.4 / 17.5 89P787P (CCA→CCT) ftsKcell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (42/26);  new base T (9/9);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACC  >  U00096/935443‑935730
                                                                                                                                             |                                                                                                                                                  
gcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagcagc                                                                                                                                             >  1:650651/1‑149 (MQ=255)
 cagcCGCCACAACCGGATGCTCAGCAGCAGCATTATCAGCATCCTCAACAACCTCATCCGCCGCAGCAGCAGTATCCCCAGCCACACCAGCCTGTTCCTACGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAgcagca                                                                                                                                            <  2:650505/149‑1 (MQ=17)
  agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCaa                                                                                                                                           <  1:650506/149‑1 (MQ=255)
   gcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAat                                                                                                                                          <  1:650507/149‑1 (MQ=255)
    cctccacaacaGGTCGCTCCGCAGCAGCATGATCAGCAGCCTCAACAACCCCATCCGCCGCAGCAGGAGTCTCAGCAGCCACCCCAGCCGCCGGCTCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata                                                                                                                                         <  2:650509/144‑1 (MQ=16)
     cGCAACTACAGGTTGCTGCGCAGCAGCAATATCAGCAGCCGCATCTTCCAGTTCCGCTGCTGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAatat                                                                                                                                        <  2:1944861/149‑1 (MQ=21)
      tcaacaaCCGTATGCTCCGCAGCAGCTATATCATCACCCTCAACAAACTCTTCCGCCGCAGCAGCAGTATCAGCAGCCACATCAGCCGGCTGCTCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATc                                                                                                                                       <  2:650520/148‑1 (MQ=18)
       cttaaacCTGTTGCACCGCAGCTTGATTATCAGCAGCCGCATCAACGATTTCCTACGCAGCCTCATTTTCATCAGCAACAACTGCCGTTGGCTACGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATca                                                                                                                                      <  2:650535/145‑1 (MQ=12)
                        gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCCGCAGCAGCAACAACCGGTTCCGCCACAGCAGCAAAATCATCAGCAGCAACAGCCgg                                                                                                                     >  2:650629/1‑149 (MQ=25)
                        gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATAAGCAGCCACAACAGCCGGTTGCGCCGCAACCCCAGTATCAGCAGAAGCATCAACAGTTTTCGCCACAGCCTCAAAATCAGCAGCAGCAACCACCga                                                                                                                     >  2:650664/1‑148 (MQ=18)
                           gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCACTATCCGCAGCATCATCCACCGGTGACGCCTCATCCGGAAAGGCATCAGATGCAACAACCga                                                                                                                     >  2:650626/1‑145 (MQ=12)
                           gcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTGGCGCCTCATCCCCAAGATCAGCAGCAGCATCAACCGGAGACGCCTCAGCCGGAAATGCATCAGCTGCAACAGCCgccga                                                                                                                  >  2:650622/1‑145 (MQ=12)
                            cagcaATAACATCATCCGCACCCTCCGCTTCCGCTGGAGGAGCATTATCTTCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATTTCAGCAGCCTCAACTACCTGTTGCTCTACAGCAGCAATATCAGCAGCCGCAACAACCGGTTgc                                                                                                                 <  2:650528/149‑1 (MQ=12)
                             agcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCAGCAACAACCGCATgcg                                                                                                                >  1:650656/1‑149 (MQ=255)
                              gcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCACTATCCGCAGCAGCAACAACCGGTTACGCCTCATCCGCAAAAGCCGCAGCAGCATCAACCgcagacgc                                                                                                               >  2:650640/1‑142 (MQ=18)
                                  tatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCATCCACAATATCAGCAGCAGCATCAACAGGTTACGCCTCATCCCCAATAGCATCAGCAGCAACAACCGAATAcgccgct                                                                                                           >  2:650620/1‑148 (MQ=21)
                                   atCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCAATATCAGCAGCAGCAACAACCGGTTCCGCCTCAGCCGCAAAATCATCAGCAGCAACAACCGCATACGCCgcag                                                                                                          >  1:650657/1‑149 (MQ=21)
                                     cagcacCCGCCACCACCCCTTCCGCCGCAGCCGCAGTATCATCAGCCACCTCAGACGTTTGCTCTGCTGGAACAATATCAGCAGCCTCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagca                                                                                                        <  2:650519/149‑1 (MQ=18)
                                          gcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAGTATCAGCAGCAGCAACAACCGGTTGCGCCACAGCAGCAATAACAGCAGCAGCATCAACCGAGTGCGCCGCAGCAGCAGt                                                                                                   >  2:650612/1‑149 (MQ=37)
                                           ccGCATCAACCTTTTCCGCTGCAGCCGCATGATCAGCAGCCACATCAGCTGTTTGCGCCGCAGCCACAATTTCAGCAGCCGCATCATCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTa                                                                                                  <  2:650516/149‑1 (MQ=18)
                                           ccGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCCGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTa                                                                                                  <  2:650515/149‑1 (MQ=255)
                                           ccGCAACAACCAGTTACGACGCAGCAGCATTATCAGCAGCAACAATATCCGTTTGCGACGCAGCAACTATATCAGCAGCCGCAACAACCGTTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTa                                                                                                  <  2:650514/149‑1 (MQ=18)
                                                acaacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCCGCAGCCGCAACAGCCGGTTGCGCCTCATCCCCACTGTCATCCGCAGCATCAACCGGATGCGccgcagccgcagatgcagc                                                                                             >  2:650671/1‑142 (MQ=21)
                                                acaacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAAGATCACCAGCCTCAACCGCCGGTTGAGCCACAGGACCACTAACAGCCGCAGCAAAAACCGCTTGCGCCGCAGCAGGAGTATc                                                                                                >  1:650630/1‑146 (MQ=25)
                                                  aacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCGCAACAGCCGGTTGCGCCGCAGCCGCAAGATCAGCAGCCGCAACAGCAGGTTGCGCGACATCAGGAATATCAGCCGCCGCATCAACCGGTTGCGCCGCAGCCGAAGGATcagctg                                                                                           >  1:650644/1‑147 (MQ=255)
                                                  aacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCAGCAGCATCAACAGCAGGTTGCGCCTCATCAGCAAAAACATCAGAAGCAACAACCGGATGCGCCGCTGCAGCAGGATcagaag                                                                                           >  2:650665/1‑149 (MQ=25)
                                                   acCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATACGCAGCCGCAACAGCCGGTGGCGCCACATCAGCACTATCCTCAGCAGCATCAACCGGCGGCGCCGCAGCCGGAGCTTCagcagc                                                                                          >  1:650632/1‑149 (MQ=21)
                                                    cctgTTGCGCCTCTGCCTAAGTATCAGCAGCGTAAACAGCCGGCTGGGCCTCATGCTTATTATCAGCTGTCTCAACTTGCGGTTGCGACACAGTTGAAATATGAGCAGCCTGAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGcc                                                                                         <  2:1663850/146‑1 (MQ=255)
                                                                cagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAGTATCAGCAGCCGCAACCGCCGGTGGAGCCTCAGCAGCACAGACGGCAGCCGCATCAACCGGATGCGCCGCAGCCGGCGAATCAGGAGCCGCCACCGccgca                                                                             >  2:650680/1‑146 (MQ=9)
                                                                 agcctcagTATCAGCAGCCACATCAGCTGTTTGCTCCGCAGCAACATTATCAGCAGCCTCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                                            <  2:650522/149‑1 (MQ=21)
                                                                  gccgcagTATCAGCAGCCACAACAGCCTGTTGCGCCGCAGCCACAATATCAGCAGCAGCAACAACCGGTTGAGCCACAGCAGCAAAAACAGCAGAAGCATCAACCGGTTGCGCCGCAGCCGCAGAATCAGCAGCCTCAACAGCCAGATg                                                                           >  2:650643/1‑149 (MQ=255)
                                                                  gccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAGGATCACCAGCAGCATCAGCCGGTTGAGCCTCATCGCCAAAAACGGCCGACGCATCAACCGAATACGCCGCAGCCGACGGAGCAGAAGCCGCTACAGCGAGATg                                                                           >  2:650621/1‑149 (MQ=255)
                                                                  gccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAAAATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGAATCAGCAGCCACAACAGCCacatg                                                                           >  2:650642/1‑145 (MQ=255)
                                                                          aTCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcac                                                                                                                                                   <  1:650525/69‑1 (MQ=255)
                                                                           tCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGATGCGCCGCAGCAGCAGTATCAGCAGCCACCACAGCCACATGCGCCACAAc                                                                  >  2:650641/1‑149 (MQ=255)
                                                                               cagcAACATCAGCCGGTTGCGCCGCAGGCACATTATCAGCAGCCGCAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                              <  2:650529/149‑1 (MQ=21)
                                                                                 gcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCAGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGACGCAACAACAGCCAGTTGCGCCACAACCGCCgg                                                            >  2:650666/1‑149 (MQ=25)
                                                                                      aacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCCGCAATATCCGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCATCAGCCACAACAGCCAGATACGCCACAACCGCAGGATAcc                                                       >  2:650653/1‑149 (MQ=21)
                                                                                       acttccGGTACCGACGCTGCCCCCATCTCAGCCTCCGCAACATCCGTTTTCGCCCCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                                                      <  2:650556/145‑1 (MQ=18)
                                                                                          gCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                                                   >  2:650615/1‑149 (MQ=255)
                                                                                          gCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                                                   <  2:650530/149‑1 (MQ=255)
                                                                                          gCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAG‑CAGTTGCGCCACAACCGCAGGATACCCTGCt                                                  >  1:650685/1‑149 (MQ=255)
                                                                                             ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                                                    >  2:650531/1‑129 (MQ=255)
                                                                                             ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                                                    <  1:650531/129‑1 (MQ=255)
                                                                                             ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                >  1:1739922/1‑149 (MQ=255)
                                                                                             ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                >  2:650647/1‑149 (MQ=255)
                                                                                             ggTTGCGCAGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAAcc                                                                                                                       >  1:1739923/1‑78 (MQ=255)
                                                                                              gTTGCGCCGCAGCCACAATATCAGCATCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                                               <  1:2507274/149‑1 (MQ=255)
                                                                                              gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCTGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCCGTTGCGCCGCAGCAGCAGTATCATCAGCCACAACAGCCAGTTGCGCCACAACCGCAGCATACCCTGCTTCa                                               >  1:650678/1‑149 (MQ=255)
                                                                                              gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCACTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACACCCGAAGGATACCCTGCTTCa                                               >  1:650648/1‑149 (MQ=255)
                                                                                              gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACAGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                                               >  2:2585390/1‑149 (MQ=255)
                                                                                              ctcgcgcCGCAGGAACATTCTCAGCACCAGCAACAACCGGTTGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCTCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCa                                               <  2:650559/146‑1 (MQ=18)
                                                                                                    ccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATcagc                                                                                                                                    <  1:650534/58‑1 (MQ=255)
                                                                                                      gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGc                                          >  1:650645/1‑146 (MQ=255)
                                                                                                      gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt                                       >  2:650659/1‑149 (MQ=255)
                                                                                                       cAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                                                     <  1:650537/118‑1 (MQ=255)
                                                                                                        aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg                                     >  1:650679/1‑149 (MQ=255)
                                                                                                          ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt                                   >  1:2334107/1‑149 (MQ=255)
                                                                                                           cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                  >  1:650662/1‑149 (MQ=255)
                                                                                                           cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                  >  1:650638/1‑149 (MQ=255)
                                                                                                           cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                  >  1:650646/1‑149 (MQ=255)
                                                                                                            acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATcagc                                                                                                                                    <  1:650539/50‑1 (MQ=255)
                                                                                                             caATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCAGCTTCATCCGCTGTTGATGCg                                >  2:1682494/1‑149 (MQ=255)
                                                                                                               atatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                              >  1:3177325/1‑149 (MQ=255)
                                                                                                               agatCAGCAGCCGCAACAACCGGTTGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                              <  2:650540/147‑1 (MQ=255)
                                                                                                                 atCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt                            >  1:650661/1‑149 (MQ=255)
                                                                                                                 atCAGAAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCAGCTGATGATGCGTAAt                            >  1:650693/1‑149 (MQ=255)
                                                                                                                  tCAGCAGCCGCAACAACCGGTACCGCCACAGAAGCAATATCAGCAGCAGCAACAACCGGTTGCGCCGCAGCAGCAGTATAAGCAGCCACAAAAGACAGTTGCGCCACAACCGCAGGATACCCTGCATCATCCGCTGTAGATGCCTaaac                           >  2:650684/1‑147 (MQ=37)
                                                                                                                   cagcagCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGcaac                                                                                                                           <  1:650553/52‑1 (MQ=255)
                                                                                                                   cagcagCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGcaac                                                                                                                           <  1:650554/52‑1 (MQ=255)
                                                                                                                     gcagcCGCAACAACCGGGTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGcgccgc                                                                                                            <  1:650555/65‑1 (MQ=255)
                                                                                                                      cagcCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTgc                                                                                                                 <  1:1982428/59‑1 (MQ=14)
                                                                                                                       agcagcAACTACCGTTTGCTCATCAGCATCATTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGAc                      <  2:650585/149‑1 (MQ=25)
                                                                                                                         ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAg                    <  2:650557/149‑1 (MQ=255)
                                                                                                                           gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGGTACCCGGCTTCATCCGCTGTTGATGCGTAATGGCGACAGc                   >  1:650698/1‑148 (MQ=255)
                                                                                                                           gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGcc                  <  1:650558/149‑1 (MQ=255)
                                                                                                                             aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgt                >  1:650634/1‑149 (MQ=255)
                                                                                                                               catccGCTTGCCCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCTGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtcc              <  1:650560/146‑1 (MQ=255)
                                                                                                                               caacCGTTTGCGCCTCTGCAGCTATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtcc              <  2:650561/149‑1 (MQ=255)
                                                                                                                                aacCGGTTGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg             <  2:650562/149‑1 (MQ=255)
                                                                                                                                aacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcag                                                                                                       >  1:650627/1‑59 (MQ=255)
                                                                                                                                 acCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccgt            >  1:650639/1‑149 (MQ=255)
                                                                                                                                  ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCaacaac                                                                                                                        <  1:650563/40‑1 (MQ=255)
                                                                                                                                   cGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTg          <  1:3294297/149‑1 (MQ=255)
                                                                                                                                    ggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACCGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGc         >  2:650660/1‑149 (MQ=255)
                                                                                                                                     gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGc         >  1:650712/1‑148 (MQ=255)
                                                                                                                                     gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa        >  2:650690/1‑149 (MQ=255)
                                                                                                                                      tcgcgcCCCAGCAGCAATATCAGCAGCCGCAACAACCGTTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCAt       <  2:650564/147‑1 (MQ=255)
                                                                                                                                         cgcCCCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATaaa    <  2:650565/149‑1 (MQ=255)
                                                                                                                                           ccACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCGTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAcc  >  2:650669/1‑149 (MQ=255)
                                                                                                                                            cacaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAcc  >  1:650655/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACC  >  U00096/935443‑935730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑