Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	U00096	935,584	0	A	T	20.0%	196.3 / 11.7	92	P787P (CCA→CCT)	ftsK	cell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (44/28);  new base T (2/16);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCG  >  U00096/935437‑935731                                                                                                                                                    |                                                                                                                                                    tGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCca                                                                                                                                                     <  1:623675/148‑1 (MQ=255)   tACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                   <  2:623676/148‑1 (MQ=255)    aCAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcag                                                                                                                                                  <  1:623677/148‑1 (MQ=255)         agcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAgcagca                                                                                                                                             <  1:623678/148‑1 (MQ=255)            cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCCCAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAgcagca                                                                                                                                             >  1:623680/1‑145 (MQ=255)               aacaacCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAGTATCAGCAGCAGCAACAACAGGTTGAGCCACAGCAGCAATATca                                                                                                                                       >  2:623710/1‑148 (MQ=255)                    ccGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCACTATCCGCAGCCGCAACAACCGGATGCGCCTCAGCCGGAAATGCagcagc                                                                                                                                  >  1:623813/1‑148 (MQ=21)                            tcggcagcagcaATATCATGAGCCTCATCAACCTTTTCCGATTCATCATCAGTTTCTTCAGGCACATCAGCGGGTTGCGCCGGAGCAACAATTTCAGCAGCCGCAACAACCGGTTGCGCCTCAGCAGCATTATCAGCAGCCGcaacaa                                                                                                                          <  2:623704/145‑1 (MQ=9)                               gcagcagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTGGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGGTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCg                                                                                                                       <  1:678202/148‑1 (MQ=255)                                     gcaATATCAGCAGCCGCAACAACCAGTTCCTCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTgcg                                                                                                                 <  2:623682/148‑1 (MQ=255)                                        atatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                   <  1:623683/111‑1 (MQ=255)                                        atatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                   >  2:623683/1‑111 (MQ=255)                                        atatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGcgccg                                                                                                              <  1:623684/148‑1 (MQ=255)                                        atatCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCACATCCCCACTATCCGCAGCATCATCAACAGTTTACGCCTCATCAGCACAAGCGGCAGCAGCATCAACCGAATGcgccg                                                                                                              >  2:623754/1‑148 (MQ=12)                                         tatCAGCAGCCGCAACAACCAGTTCCTCCGCAGCAGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGcgccgc                                                                                                             <  2:623685/148‑1 (MQ=255)                                            catcatacGCACCAACCCGTTCCGCCGCAGCCGCAGTATCATCACCCACAGCAGACGAATGCGCCGCTGCAACATGTTCAGCAGCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagc                                                                                                          <  2:623692/141‑1 (MQ=18)                                            cagcagCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCAGCAGCCGCAACCGCCGGTTGCGCCAAAGCCGGAAATACAGCTGCGGCCACCGCCGGTTGggccaaagc                                                                                                          >  2:623728/1‑143 (MQ=17)                                              gcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCTCAGCCGCACTATCCGCAGCAGCATCAACCGGTTGCGCCACAGCCCGAAAAGCATCAGCAGCAACAACCGCATCCGCCgcagccg                                                                                                        >  2:623824/1‑146 (MQ=18)                                               cagcCTCACCAACCAGTTCCGCCGCTGCCGCAGTATCATCAGGAACATCAGAGGTTTGCGCTGCAGGAACATTATCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcagc                                                                                                       <  2:623687/148‑1 (MQ=18)                                               cagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcagc                                                                                                       <  1:623686/148‑1 (MQ=255)                                                  ccGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCCGCAGCCTCATCAGCCGGTTGAGCCAAATCAGCAAAAGCATCCGACGCATCAACCGCATTCGCCgcagcagactg                                                                                                    >  2:623829/1‑144 (MQ=18)                                                   cGCATCAACCTGATCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGCTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTa                                                                                                   <  1:623689/148‑1 (MQ=255)                                                      aacaacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACACCAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATca                                                                                                <  2:623690/148‑1 (MQ=255)                                                        caacCATTTCCGCCGCAGCCGCAGTATCATCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATcagc                                                                                              <  2:623691/148‑1 (MQ=255)                                                        caacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAata                                                                                                                                          >  1:623753/1‑104 (MQ=255)                                                            cAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCTCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGcc                                                                                          >  2:623781/1‑148 (MQ=255)                                                                 ccgccgCAGCCGCAGTATCAGCAGCCACATCAGCCGCTTGCGCCGCAGCAACATTATCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCacaac                                                                                     <  1:623696/148‑1 (MQ=18)                                                                 ccgccgCAGCAGGAGTATCAGCACCAACACCAGCCGTCTGCGCCGCAGCAACATTTTCAGCAGCCGCAAGAACCGGGTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCacaac                                                                                     <  1:623695/148‑1 (MQ=18)                                                                  cgccgcAGCAGGATTATCAGCAGGAACACCATCATTTTGCTCTGCTGGAACATTTTCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCacaaca                                                                                    <  1:623717/148‑1 (MQ=9)                                                                    ccgcagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACATCAGCAATATCAGCAGCCGCAACAACCGGATGCGCCGCAGCCGCAGAATCAGCAGCCACAACAGc                                                                                  >  2:623780/1‑148 (MQ=255)                                                                    ccgcagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGc                                                                                  >  1:623773/1‑148 (MQ=255)                                                                         gcagcatTATCATCAGCCTCATCATCCGTATGCGCTGCAGCATGATTATCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCATTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                                             <  2:623711/145‑1 (MQ=14)                                                                             cagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAGCCGGTTGCGCCACATCAGCAATATCAGCAGCCGCAACAACCGGATGCGCCGCAGCAGCAGTATCATCAGCAACAACAGCCAGATgcgc                                                                         >  2:623803/1‑148 (MQ=18)                                                                                 aTCAGCAGGCACATCAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcac                                                                      <  1:623697/147‑1 (MQ=255)                                                                                  tCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                                                    <  2:623699/148‑1 (MQ=255)                                                                                     gcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCCGCAGCCGCAACAACCGGATGCGCCGCAGCAGCAGTATCAGCAGGAACAACAGCCAGTTGCGCcccacacg                                                                 >  1:2893791/1‑144 (MQ=21)                                                                                     gcagcCACAACAGCCGGCTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                                                 <  2:623702/148‑1 (MQ=255)                                                                                        gcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAg                                                              >  2:623779/1‑148 (MQ=255)                                                                                           ttaacaGCCTCTGGCCCCCCCGCCACATCAACTGCAGCCGTAACATCCGGTCGTCCCACATCATCATTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGAt                                                           <  2:623755/146‑1 (MQ=17)                                                                                           acaacaGCCGGTTGCGCCGCAGCCACATTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCg                                                                                                                       <  1:623703/88‑1 (MQ=255)                                                                                            caacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGAt                                                           >  1:623834/1‑147 (MQ=255)                                                                                              acaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcag                                                                                                           >  1:623805/1‑97 (MQ=255)                                                                                                  ccGGTTGCGCCGCAGCCCCAATATCAGCAGCCGCAACACCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCCGCAGCAGTATCAGCCGCCACAACAGCCAGTGGCGCCACAACCGCAGGATACCCTg                                                     >  1:623807/1‑147 (MQ=21)                                                                                                     gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTc                                                 >  2:623802/1‑148 (MQ=255)                                                                                                       cgcgcCGCAGCCCCATTTTTAGCATCCGCATCCTCCGCCGTCGGCCCAGCAGCAATATCAGCAGCCGCAACATCCTGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCAt                                               <  2:623737/147‑1 (MQ=18)                                                                                                         cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACACCCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                                             >  2:623800/1‑148 (MQ=255)                                                                                                         cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTa                                                                                                   <  1:623706/94‑1 (MQ=255)                                                                                                         cgccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGGTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCGTCATc                                              <  1:623705/147‑1 (MQ=255)                                                                                                             gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcag                                                                                                                                                  <  1:623707/43‑1 (MQ=255)                                                                                                                gCCACATTTTCAGCAGCCGCAACATCCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCTCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg                                      <  2:623708/148‑1 (MQ=21)                                                                                                                gCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg                                      >  1:623747/1‑148 (MQ=255)                                                                                                                 ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg                                      >  1:623859/1‑147 (MQ=255)                                                                                                                 ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa                                     >  1:623809/1‑148 (MQ=255)                                                                                                                 ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa                                     >  2:623808/1‑148 (MQ=255)                                                                                                                 ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa                                     >  1:623736/1‑148 (MQ=255)                                                                                                                 ccACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa                                     >  1:623778/1‑148 (MQ=255)                                                                                                                  cacaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCGGTTGAt                                    >  2:623828/1‑148 (MQ=255)                                                                                                                   tcatttacagcagCTGCAACTACCGGTTGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                   <  2:623719/141‑1 (MQ=21)                                                                                                                   acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                   >  1:623845/1‑148 (MQ=255)                                                                                                                   acaATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg                                   >  2:623799/1‑148 (MQ=255)                                                                                                                      atatCAGCAGCCGCAACAACCGGTTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGt                                <  1:623710/148‑1 (MQ=255)                                                                                                                      atatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGt                                >  2:623806/1‑148 (MQ=255)                                                                                                                       tttCAGCAGCCGCATCACCCGGTGGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                               <  2:623716/146‑1 (MQ=21)                                                                                                                       tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                               >  1:623796/1‑148 (MQ=255)                                                                                                                       tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATACGCGGTTGATGCGTa                               >  2:623838/1‑148 (MQ=255)                                                                                                                        atCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTaa                              <  1:1301779/148‑1 (MQ=255)                                                                                                                          catcagcCGCATCATCCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg                            <  2:623723/145‑1 (MQ=21)                                                                                                                          cagcagCCGCAACAACCGGTTGCTCCACAGCAGCAATATCagcagc                                                                                                                                  <  1:623720/46‑1 (MQ=255)                                                                                                                          cagcagCCGCAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGcaaca                                                                                                                           <  1:623722/53‑1 (MQ=255)                                                                                                                          cagcagCCGCAACAACCGGTTGCGCCACAGCAGCaa                                                                                                                                            <  1:623718/36‑1 (MQ=255)                                                                                                                          cagcagCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCTGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg                            <  1:623721/148‑1 (MQ=255)                                                                                                                           aggagcCGCAACAACTGGTTGCTCCCCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATgg                           <  2:623724/145‑1 (MQ=255)                                                                                                                            gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGt                                                                                                                     <  1:623725/57‑1 (MQ=255)                                                                                                                            gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGc                          >  2:1530420/1‑148 (MQ=255)                                                                                                                             cagcCGCAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTg                                                                                                                   <  1:623726/58‑1 (MQ=255)                                                                                                                               gcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGAc                       >  2:623865/1‑148 (MQ=255)                                                                                                                                ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                      >  2:623825/1‑148 (MQ=255)                                                                                                                                ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                      >  2:623857/1‑148 (MQ=255)                                                                                                                                ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                      >  2:623841/1‑148 (MQ=255)                                                                                                                                ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                      >  2:623826/1‑148 (MQ=255)                                                                                                                                ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                      >  2:623893/1‑148 (MQ=255)                                                                                                                                 ctcaacaaCCGGCGGCGCCTCAGCCGCACTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCATAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAg                     <  1:623728/146‑1 (MQ=25)                                                                                                                                 cGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAg                     >  2:623758/1‑148 (MQ=255)                                                                                                                                   caacaaCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCG‑‑ATGGCGACaaaaaa                 >  2:1229687/1‑143 (MQ=255)                                                                                                                                    aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccg                  >  2:623812/1‑148 (MQ=255)                                                                                                                                    aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccg                  >  1:623794/1‑148 (MQ=255)                                                                                                                                        acCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg              >  1:623819/1‑148 (MQ=255)                                                                                                                                        acCGGGTGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg              <  1:623731/148‑1 (MQ=255)                                                                                                                                         ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccgt             >  2:623811/1‑148 (MQ=255)                                                                                                                                          cttttGCTCCTCAGCATCATTATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGtt            <  2:623757/145‑1 (MQ=37)                                                                                                                                             ttGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa         <  2:623733/148‑1 (MQ=255)                                                                                                                                                cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATaa      <  2:623734/148‑1 (MQ=255)                                                                                                                                                cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATaa      <  2:623735/148‑1 (MQ=255)                                                                                                                                                cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATaa      >  1:623816/1‑148 (MQ=255)                                                                                                                                                  ccACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGCCAGCCGTCCGTTGCATAAAc    >  1:623889/1‑148 (MQ=255)                                                                                                                                                  ccACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc    <  2:623736/148‑1 (MQ=255)                                                                                                                                                   cacaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAcc   <  2:1839336/148‑1 (MQ=255)                                                                                                                                                    acaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcg  >  2:623818/1‑148 (MQ=255)                                                                                                                                                    acaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcg  >  2:623830/1‑148 (MQ=255)                                                                                                                                                    acaGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACcg  >  2:623854/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAACCG  >  U00096/935437‑935731

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: ‑