Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	41.0%	111.7 / 13.8	61	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base G (25/0);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 7.09e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366377‑366674                                                                                                                                                         |                                                                                                                                                      aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGggtg                                                                                                                                                      <  2:953475/151‑1 (MQ=255)  aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGC‑G‑GCAgggggg                                                                                                                                                     >  1:430561/1‑151 (MQ=255)  aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTT‑TTTGGC‑G‑CCAgggggg                                                                                                                                                     >  2:34896/1‑151 (MQ=255)     ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGG‑G‑CCCGGGGGGgtt                                                                                                                                                  >  1:46578/1‑151 (MQ=255)        ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGC‑G‑GCCGGGGGGgtttt                                                                                                                                                >  1:814780/1‑150 (MQ=255)              ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCa                                                                                                                                         <  1:540010/151‑1 (MQ=255)                 aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCAc                                                                                                                                        <  2:437911/149‑1 (MQ=255)                  aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGT‑TTTGGG‑G‑GCAGGGGGGttttttttttttcc                                                                                                                                       >  1:195659/1‑145 (MQ=255)                    gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTTCCAGTg                                                                                                                                   >  1:303805/1‑151 (MQ=255)                      gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGG‑TTTGGC‑C‑CGAGGGGGGGtttttttttttac                                                                                                                                       >  1:864982/1‑141 (MQ=255)                       tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGT‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCAggg                                                                                                                                   >  2:409192/1‑146 (MQ=255)                        gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                               <  1:179489/151‑1 (MQ=255)                           cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑CTCTGCGG‑CCGGGGGGGtttttttttttataaagagcccgg                                                                                                                             >  1:497751/1‑137 (MQ=255)                           cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCCGTGGGAAggg                                                                                                                            >  1:657373/1‑151 (MQ=255)                                  acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGG‑G‑CCGGGGGGGTTTTT‑TTTTTACCCGTGAGACGGGCaa                                                                                                                         >  1:723165/1‑147 (MQ=255)                                           gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGC‑G‑CCAGGGGGGGTTTTTTTTTCTCCCGAGAGACGGGGAAAAGCTGATTGccc                                                                                                            >  1:761955/1‑151 (MQ=255)                                               tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGG‑G‑CCAGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAAAAGCCGATTGTCCtt                                                                                                          >  2:528019/1‑149 (MQ=255)                                               tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑GCCCGGGGGGGTTTTTTTTTTCCCCCG‑GAGACCGGCAACAGCTGATTGccc                                                                                                            >  1:740316/1‑147 (MQ=255)                                                gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTT‑TTTGGG‑G‑CCACGGGGGTTTTTTTTTTTCTCCGGGGGACAGGGAACAGCTGATTTCCCCTCAc                                                                                                       >  1:17842/1‑151 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  2:812613/151‑1 (MQ=255)                                                          cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑G‑CCCGGGGGGGTTTTTTTTTTTCCAGGGAGACGGGGAACAGCTGATTTGCCTTTCCCCCCCGGccc                                                                                             >  1:222954/1‑151 (MQ=255)                                                               tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCCGATTGCCCTTCACCGCCTGGgcc                                                                                             >  2:481320/1‑146 (MQ=255)                                                                    aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                   <  2:743272/151‑1 (MQ=255)                                                                      tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTT‑TTGGGG‑G‑CCCGGGGGGGTTTTCTTTTTTCCAGTGGGGCGGGCAACCGCTGAGTGCCCCTCACCCCCCGGgcccggggggggt                                                                                   >  1:481018/1‑139 (MQ=255)                                                                       cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                <  1:898509/151‑1 (MQ=255)                                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGG‑G‑CCCGGGGGGGTTTTCTTTTCACCAATGAGAAGGGCAACAACTGATTGCCCTTCACCCCCCGGCCCCGGGAGAGGTgc                                                                                 >  1:485951/1‑149 (MQ=255)                                                                          gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  1:726479/151‑1 (MQ=255)                                                                           aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTTTTTTTACCAGTGGGAAGGGGAACAGCTGATTGCCCCTTACCCCCTGGCCCTGcgggaggt                                                                                   >  2:898951/1‑142 (MQ=255)                                                                              ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                         <  1:871811/151‑1 (MQ=255)                                                                                 gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTA‑TTGGGG‑G‑GCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCCGCCCCGGGAGAGGTGCAGCAAGCgg                                                                        >  2:118792/1‑149 (MQ=255)                                                                                 gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                       <  2:854369/150‑1 (MQ=255)                                                                                   cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGGA‑TTGGGC‑G‑GCCGGGGGGGTTTTTTTTTTTCCCGTGAGACGGGCAACAGCTGGTTGCCCCTCACCGCCCGGCCCCGAGAGAGGTGGAGCAAACGGGcc                                                                     >  1:746727/1‑150 (MQ=255)                                                                                   cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                    <  1:682159/151‑1 (MQ=255)                                                                                       ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGC‑TTTGGC‑G‑CCCGGGGGGTTTTTTTTTTCCCCACTGAGACGGGCAAAAGCGGATTGGCCTTCACCGCCCGGCCCTGCGAGAGTTGGACCAAACggggcac                                                                   >  1:756118/1‑143 (MQ=255)                                                                                             gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTT‑TTGTGG‑G‑CCAGGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCCGATTGCCCTTCTCCCCCCTGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGC‑TGGgtt                                                          >  1:40829/1‑150 (MQ=255)                                                                                             gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGttt                                                          <  2:820699/150‑1 (MQ=255)                                                                                              cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑ATTAGC‑G‑CCCGGGGGGTTTTTCTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCCTCCCCGCCCGGCCCTGGGAGAGTTGCAGCAAGCGGGCCACGC‑CGGgtt                                                          >  1:612797/1‑149 (MQ=255)                                                                                                 tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGccc                                                      <  2:662561/150‑1 (MQ=255)                                                                                                  aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCca                                                    <  2:959895/151‑1 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGC‑G‑CCAGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAGCTGATTGCCCTTCACCGCCCGGCCCTGGGAGAGGTGCAGCAAGCGGGCCACGC‑TGGTTTGCCCcagc                                                  >  2:45/1‑150 (MQ=255)                                                                                                        tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCCGGGGGGTTTTTTTTTTTACCAGTGAGACCGGGAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCACGCTTGGTTTCcccccccg                                                 >  2:772864/1‑145 (MQ=255)                                                                                                              aaCGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAAt                                        <  2:41457/151‑1 (MQ=255)                                                                                                                cgcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATcc                                      <  2:250927/151‑1 (MQ=255)                                                                                                                cgcgcgGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATcc                                      <  1:504728/151‑1 (MQ=255)                                                                                                                          gagGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                            <  1:481306/151‑1 (MQ=255)                                                                                                                          gagGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                            <  2:212537/151‑1 (MQ=255)                                                                                                                            ggCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                          <  1:117401/151‑1 (MQ=255)                                                                                                                            ggCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                          <  2:743299/151‑1 (MQ=255)                                                                                                                            ggCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                          >  2:644108/1‑151 (MQ=255)                                                                                                                                gTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  2:566010/151‑1 (MQ=255)                                                                                                                                gTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑AGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                       >  1:402977/1‑150 (MQ=255)                                                                                                                                   tGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                   >  2:871894/1‑151 (MQ=255)                                                                                                                                    gCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  >  1:797515/1‑151 (MQ=255)                                                                                                                                    gCGTA‑TTGGGC‑G‑CCAGGGTCGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  <  2:323815/151‑1 (MQ=255)                                                                                                                                      gTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                <  2:557315/151‑1 (MQ=255)                                                                                                                                      gTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                <  1:34898/151‑1 (MQ=255)                                                                                                                                       tA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               >  2:666747/1‑151 (MQ=255)                                                                                                                                           tGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  2:131142/1‑149 (MQ=255)                                                                                                                                           tGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  2:389333/1‑149 (MQ=255)                                                                                                                                           tGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa             <  1:56495/150‑1 (MQ=255)                                                                                                                                            gggC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata           <  2:755041/151‑1 (MQ=255)                                                                                                                                                    cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     >  2:474162/1‑151 (MQ=255)                                                                                                                                                    cAGGGTGGTTTTTCTTTTCACCAGGGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGGGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAGATaaaa      >  1:159669/1‑148 (MQ=255)                                                                                                                                                     aGGGTGGTTTTTCTTTTCACCAGTGGGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  1:450745/151‑1 (MQ=255)                                                                                                                                                       ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  2:222958/151‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                      AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGC‑TGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366377‑366674

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