Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,060,280	0	T	.	74.1%	108.8 / 70.0	58	intergenic (+50/‑167)	bipA/yihL	50S ribosomal subunit assembly factor BipA/putative DNA‑binding transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base T (5/10);  new base . (20/23);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGCCGG  >  NC_000913/4060036‑4060525                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      aCCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTAGCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTc                                                                                                                                                                                                                                                <  1:621863/251‑1 (MQ=255)  cctatatgcggcccgccGGTAATGACGAATACGTTGTTCTGATACCGCCTTTCCGAAAGACGGTGGAACAAGCTCTGGAGTACATCGATGATGACGAACTGGTAGAAGTGACTCCGGCCTCTATATGTATTCGTAAACGTCGCCTGACTTAAAACGATAGTCGCCACGCCAACCGCGCACCTAAAGACGATTAATTTCGTTTTTTGGTCATAAAAAACCTGCCAGCGATGGCAG‑TTTCTTTTTGCCTGTCg                                                                                                                                                                                                                                               <  2:593366/236‑1 (MQ=255)         acGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTGT‑CCTGTCGTGAgcgc                                                                                                                                                                                                                                        >  2:7968/2‑251 (MQ=255)           gTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCAGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGcac                                                                                                                                                                                                                                      >  1:425090/1‑251 (MQ=255)           gTGCTTCCGGAACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGcac                                                                                                                                                                                                                                      >  2:423176/1‑251 (MQ=255)                   ggTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACTATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTTTCGTGAGCGCACACCAATGc                                                                                                                                                                                                                              <  2:395107/251‑1 (MQ=255)                           cGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                                                                                                                                                                                       <  1:40708/251‑1 (MQ=255)                           cGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                                                                                                                                                                                      <  1:375980/251‑1 (MQ=255)                              aGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                                                                                                                                                                                   >  1:510577/1‑251 (MQ=255)                                 cGTTGTTCTGGTTACGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCttt                                                                                                                                                                                                                >  2:95162/1‑250 (MQ=255)                                              ccGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                                                                                                                                                                                    <  2:317902/250‑1 (MQ=255)                                              ccGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGGCTGTGGtt                                                                                                                                                                                                   <  1:215136/251‑1 (MQ=255)                                               cGCCTAACCGCATGACTCTGGAACAAGCTCTGGAGTTCAGCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                                                                                                                                                  <  1:98231/251‑1 (MQ=255)                                                      ccGCATGACTCTGGAACAAGCGCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGGAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGACTGTGGTTACGTCAAt                                                                                                                                                                                           <  2:522661/251‑1 (MQ=255)                                                       cGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAGTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                                                                                                                                                                                          <  1:423176/251‑1 (MQ=255)                                                         cATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                                                                                                                                           >  1:42035/1‑249 (MQ=255)                                                                  ggAACAAGCTCGGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                                                                                                                                                 <  2:98833/249‑1 (MQ=255)                                                                              tgAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTTAGCGCACACCAAGGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc                                                                                                                                                                     <  1:85509/249‑1 (MQ=255)                                                                               gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg                                                                                                                                                                   >  2:347847/1‑251 (MQ=255)                                                                                   tCATCGATGATGACGAACTGGTAGAAGTTACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACTGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaaa                                                                                                                                                               <  1:82853/250‑1 (MQ=255)                                                                                       cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCa                                                                                                                                                           <  1:420063/250‑1 (MQ=255)                                                                                        gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAACAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                                                                                                                                          >  2:305880/1‑250 (MQ=255)                                                                                           gatgaCGATCTTGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACgt                                                                                                                                                      <  2:357313/251‑2 (MQ=255)                                                                                               aCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                                                                                                                                                   <  1:438884/251‑1 (MQ=255)                                                                                                  aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAACGTCATACGCCGTTAAt                                                                                                                                               >  2:485603/1‑251 (MQ=255)                                                                                                     tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTCCCAGCGAAGGCAGGTTTTTTTTTGCTTGTCGTGAGCGCACACAAATACGAGTATCTCCCTTATAAGTCTGTGATTACGTCAATAGATAGCTTTCTCACAATTATCTGCACAGTCCTACGCCGTTACTTg                                                                                                                                             >  1:33138/1‑251 (MQ=255)                                                                                                      ggTAGAAGTTACTTCGACCTCTATCCTTATTTGTAAACGTCATCGGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGGGAGTGCACACCAATGCTAGGATCTCCCTTATAATTCTGGGGTTACGTCTATAGAGAGCTTTGTCATAATTATCTGCAAAGTCATACGGCGTTAATTGc                                                                                                                                            <  2:401544/251‑1 (MQ=255)                                                                                                      ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTg                                                                                                                                             <  1:7857/250‑1 (MQ=255)                                                                                                          gAAGTGACTCCGACTTCTATTCGTATTCGTAAAAGTCACCTGACGGAAAACGATCGTCCCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAATCACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACACTTAACTGCACAGCCATACGCCTTTAATTGctttct                                                                                                                                       >  2:8506/1‑251 (MQ=255)                                                                                                          gAAGTGACTACGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttct                                                                                                                                       >  2:48389/1‑251 (MQ=255)                                                                                                                  tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGCTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGCCACCATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc                                                                                                                               >  2:231087/1‑251 (MQ=255)                                                                                                                     gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAATACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt                                                                                                                             <  1:384284/251‑1 (MQ=255)                                                                                                                     gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgta                                                                                                                            <  1:467366/251‑1 (MQ=255)                                                                                                                             tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta                                                                                                                      >  1:476509/1‑249 (MQ=255)                                                                                                                                  aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATg                                                                                                                 <  2:49601/250‑1 (MQ=255)                                                                                                                                     cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCCGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGGAAGACGCTTc                                                                                                             <  2:164853/250‑1 (MQ=255)                                                                                                                                      gTAAACGTCCCCTGACGGCAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAACAAACCCGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGCGCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCCTACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCAt                                                                                                           >  2:253147/1‑251 (MQ=255)                                                                                                                                       tAAACGGCACCTGACGGAAAACGATTGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGATTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTACCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATc                                                                                                          >  2:530044/1‑251 (MQ=255)                                                                                                                                                    aCGGCAAACGATCGTCGCCGCGCCAACCGCGCACCGAATGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTTGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCTTAAGCGTAAGATGCTTCATCTGGTTTAAACCaa                                                                                             <  1:38721/251‑1 (MQ=255)                                                                                                                                                     cGGAAAACGCTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGTTTCATCTGGTTTAAACCaa                                                                                             <  2:52301/251‑1 (MQ=255)                                                                                                                                                       gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAg                                                                                          <  2:84096/251‑1 (MQ=255)                                                                                                                                                          aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAATATGCTTCATCTGGTTTAAACCATAAGGAt                                                                                       <  2:563561/251‑1 (MQ=255)                                                                                                                                                            cGATCGTCGCCGCGCCAACCGCGCACCGCAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTAATCTGGTTTAAACCAAAAGGAtt                                                                                      >  2:400888/1‑250 (MQ=255)                                                                                                                                                                cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACa                                                                                 <  1:618212/251‑1 (MQ=255)                                                                                                                                                                 gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTCTAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGAAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACa                                                                                 >  1:1528/1‑250 (MQ=255)                                                                                                                                                                   cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAAt                                                                               <  1:243926/250‑1 (MQ=255)                                                                                                                                                                   cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATg                                                                              <  1:350596/251‑1 (MQ=255)                                                                                                                                                                                 cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaatcaatc                                                                >  2:195578/1‑251 (MQ=255)                                                                                                                                                                                 cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGaatcaatc                                                                >  1:525750/1‑251 (MQ=255)                                                                                                                                                                                   cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATccc                                                              >  1:136402/1‑250 (MQ=255)                                                                                                                                                                                             cGATTCATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGaaa                                                     >  1:506444/1‑250 (MQ=255)                                                                                                                                                                                                        gtcttttctgcttAAAAAACCTGCCAGAGATGGCAGGTTTTCTTTTGCCTGGCGTGAGCCCACACCCATGCGCGTATCTACCTTATACGTCTGTGGTTACGTCCATAGAGAGCTTTGTCACAATTATCTGCAAAGTCTTACGCCGTTAATTGCTTTCTTTTTTGGAGTAAGCGTAAGATGTTTCATCTGGTTTAAACCAAAAGGTTTCAACAATGGCGGAGAATCAATCCACCTTCGAAAATGCaaaaa                                            >  2:367726/13‑248 (MQ=255)                                                                                                                                                                                                                  cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCg                               >  2:224619/1‑251 (MQ=255)                                                                                                                                                                                                                     aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCggtg                            <  2:448108/251‑1 (MQ=255)                                                                                                                                                                                                                          acatGCCAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTg                        >  2:167096/4‑250 (MQ=255)                                                                                                                                                                                                                               ccAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGAt                  <  2:107053/251‑1 (MQ=255)                                                                                                                                                                                                                                cAGCGATGGCAGGTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATg                 <  1:478957/251‑1 (MQ=255)                                                                                                                                                                                                                                    cgatggccggtattttttGCCTGTCGTGAACGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTACGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATACAAAAGAGAAACTGGATCGGTGGTTGAAAGAGGGc               >  2:102389/13‑250 (MQ=255)                                                                                                                                                                                                                                        gCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGTGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAAGCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATcac          <  2:42035/251‑1 (MQ=255)                                                                                                                                                                                                                                                ttttttGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGGGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCCAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCCTCACCACGCCgg  >  2:355379/1‑251 (MQ=255)                                                                                                                                                                                                                                                ttttttGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTGCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAAGGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGCCg   <  2:29529/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ACCAACATGCGTGCTTCCGGTACTGACGAAGCCGTTGTTCTGGTTCCGCCTATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAAGATGCTTCATCTGGTTTAAACCAAAAGGATTAAACAATGGCGGAGAATCAATCCACCGTAGAAAATGCAAAAGAGAAACTGGATCGGTGGTTGAAAGATGGCATCACCACGCCGG  >  NC_000913/4060036‑4060525

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