Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	708,774	0	G	A	78.0%	80.4 / 19.6	41	A147A (GCG→GCA)	chiP	chitooligosaccharide outer membrane channel
Reads supporting (aligned to +/- strand):  ref base G (7/2);  new base A (14/18);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGAC  >  NC_000913/708530‑709020                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       accTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCa                                                                                                                                                                                                                                                  <  1:472816/251‑1 (MQ=255)  ccTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAg                                                                                                                                                                                                                                                 <  2:51896/251‑1 (MQ=255)       ctcACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTACGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTa                                                                                                                                                                                                                                             >  1:61582/1‑250 (MQ=255)         cACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGTGCAGGTTAc                                                                                                                                                                                                                                            <  2:388856/249‑1 (MQ=255)         cACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCGATGACGAAGACTGGGCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTa                                                                                                                                                                                                                                             <  2:61582/248‑1 (MQ=255)          ccTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACAtt                                                                                                                                                                                                                                         <  2:338338/250‑1 (MQ=255)             ccACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCGCCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAg                                                                                                                                                                                                                                      <  2:258402/251‑1 (MQ=255)                   ggAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAAc                                                                                                                                                                                                                                 <  2:90007/250‑1 (MQ=255)                                  aTTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTa                                                                                                                                                                                                                 <  2:329108/251‑1 (MQ=255)                                    tttCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGGTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAg                                                                                                                                                                                                                <  1:164340/250‑1 (MQ=255)                                                     tgctgATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCttt                                                                                                                                                                                               >  1:378874/1‑250 (MQ=255)                                                              gTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGt                                                                                                                                                                                      <  1:154995/250‑1 (MQ=255)                                                               ttCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGAGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTAc                                                                                                                                                                                    >  1:168009/1‑251 (MQ=255)                                                                       tGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAAGTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGTAGCTTTATGCCAGGTACTTATCa                                                                                                                                                                              <  2:318359/249‑1 (MQ=255)                                                                        gATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGCAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAggg                                                                                                                                                                           >  1:422530/1‑251 (MQ=255)                                                                        gATATTGCTGCGTTTACGGCGATTGAAATGGAGGAAAACGGGGGCAGCTCCCACCAGAACGGAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGAATGGTCCGGCGACAAAAGCGGTATAAGCCTGAATAAAGCTGCGGCCAAATTTAAAGACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAAAGCTGTTAGCGCAGCACTGGAGCTTTATGCCAGGTACTTATCAggg                                                                                                                                                                           >  1:170867/1‑251 (MQ=255)                                                                                       aCGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGc                                                                                                                                                            <  1:341880/251‑1 (MQ=255)                                                                                           cGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAAt                                                                                                                                                        <  2:10434/251‑1 (MQ=255)                                                                                                                gCGACAGCTCCCACACGAACGAAATCGCGTTTTCAAAAAGTAATCAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAACCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGAAGGGGCGAATTTTGATTACGGCGATGCTGGt                                                                                                                                   >  2:46433/1‑251 (MQ=255)                                                                                                                  gACAGCTCCCACCCGAACGAAAGCGCGTTTACAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTGTAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGc                                                                                                                                 <  2:178811/251‑1 (MQ=255)                                                                                                                           cACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGTTGCGTTGAGttt                                                                                                                        <  1:358016/251‑1 (MQ=255)                                                                                                                              ccGAAGGAAATCGCGTTTTAAAAAAGTAATAAAGCCTATGACGCAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGCATTTTGATTACGGCGATGCTGGTGCGTTGAGTTtctc                                                                                                                     >  2:232733/1‑251 (MQ=255)                                                                                                                                           cgTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGAcc                                                                                                        >  1:334007/1‑251 (MQ=255)                                                                                                                                            gTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTGGGGTTTCTCCTCCATTTGGAca                                                                                                        >  2:483926/1‑249 (MQ=255)                                                                                                                                                         aataaAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATAAAAAgc                                                                                          >  1:258830/1‑251 (MQ=255)                                                                                                                                                               gCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTg                                                                                    <  1:230502/251‑1 (MQ=255)                                                                                                                                                                   aTGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGATGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCAAACTGGTCAAACGCTGTTAGCGCAGCACTGGAGCTTTATGCCAGGTACTTATAAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGATGGGGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGACGTggaat                                                                                >  1:42851/1‑248 (MQ=255)                                                                                                                                                                   aTGACGAAGACTGGGCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCAt                                                                                <  2:258830/251‑1 (MQ=255)                                                                                                                                                                                  gcGGCGACAAAAGCGGTATAAGTCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGTCAACTTGTCAGACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGTGACTTATCAGGGGGCTGAAGCCGTGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCGACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGa                                                                  <  1:329663/249‑1 (MQ=255)                                                                                                                                                                                    ggCGACAAACGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGATGGTGCGGTGAGTTTCTCCTACATGGGGACCAACGAATACAAAGCGCCGTGGCATCTGGGAATGGATGAGtt                                                               >  2:356272/1‑251 (MQ=255)                                                                                                                                                                                       gACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGtttt                                                             >  1:26725/1‑250 (MQ=255)                                                                                                                                                                                        aCAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTAt                                                           >  2:293548/1‑251 (MQ=255)                                                                                                                                                                                                         cTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCCGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAaccac                                          >  1:71849/1‑251 (MQ=255)                                                                                                                                                                                                         cTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAaccac                                          >  1:6625/1‑251 (MQ=255)                                                                                                                                                                                                             ataAAGCTGCGCCCAAATTTAAATACGGTCCGGTTTGGGCGAGGACAGGATACATTCAGCCAACTGGTCAAACGCTGTCCGCGCCGCACTGGAGCTTTATGCAAGGTACTTATCAGCGGGCGGAAGACGGGGCGAATTTTGCTTACGGCGATGCTGGTGCGTTGAGTTTCTCCGACAGGTGGACCAACGATTACAAAGCGCCGTGGTAACTGGAAATGGAAGAGCTTTTTCAGAACGATAAAACcaacaa                                       >  2:293276/1‑250 (MQ=255)                                                                                                                                                                                                               aaaGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCa                                                                                 >  1:150357/1‑206 (MQ=255)                                                                                                                                                                                                               aaaGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCa                                                                                 <  2:150357/206‑1 (MQ=255)                                                                                                                                                                                                                                        gTCCGGTTTGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCg           <  1:45550/251‑1 (MQ=255)                                                                                                                                                                                                                                            ggTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAAt       >  1:196649/1‑251 (MQ=255)                                                                                                                                                                                                                                              tttGGGCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATa      >  1:424461/1‑250 (MQ=255)                                                                                                                                                                                                                                                 gggCAAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGAc  >  1:348156/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ACCTTTCTCACTCCACCTGGAATGCCAACCTCGATTTTCAGTCCGGTTATGCTGCTGATATGTTCGGCCTTGATATTGCTGCGTTTACGGCGATTGAAATGGCGGAAAACGGCGACAGCTCCCACCCGAACGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCCAACTGGTCAAACGCTGTTAGCGCCGCACTGGAGCTTTATGCCAGGTACTTATCAGGGGGCGGAAGCCGGGGCGAATTTTGATTACGGCGATGCTGGTGCGTTGAGTTTCTCCTACATGTGGACCAACGAATACAAAGCGCCGTGGCATCTGGAAATGGATGAGTTTTATCAGAACGATAAAACCACCAAAGTTGATTATCTGCACTCCTTTGGGGCGAAATACGAC  >  NC_000913/708530‑709020

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