Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	48.3%	86.3 / 20.6	59	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base G (28/0);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366380‑366676                                                                                                                                                   |                                                                                                                                                       gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGggtgggt                                                                                                                                                    >  2:802674/1‑149 (MQ=255) gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTGGTATTGGGCGGCCgggggg                                                                                                                                                      >  1:344580/1‑149 (MQ=255)   cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGgttt                                                                                                                                                  >  2:780397/1‑151 (MQ=255)      gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                               <  2:560829/151‑1 (MQ=255)          gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGTTTGGGCGCCCGGGTGGGttttcttttt                                                                                                                                           >  2:400934/1‑150 (MQ=255)               aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCCGGGGGGGTTTTCTTTTTAcc                                                                                                                                        >  1:88947/1‑149 (MQ=255)               aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCAGGGGGGGTTTTTTTTTTAcc                                                                                                                                        >  2:369514/1‑149 (MQ=255)                 gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGTTTTTTTTTTCACCAGTg                                                                                                                                    >  1:536590/1‑151 (MQ=255)                 gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                     <  1:695123/150‑1 (MQ=255)                 gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                     <  1:664859/150‑1 (MQ=255)                   gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGt                                                                                                                                     >  2:273707/1‑148 (MQ=255)                    tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAATg                                                                                                                                    >  1:371775/1‑148 (MQ=255)                         tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                             >  2:136624/1‑150 (MQ=255)                         tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCCGGGGGGGTTTTTTTTTCCCCAGGGAGACggg                                                                                                                             >  1:577351/1‑150 (MQ=255)                          aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTCTTTGGCGCCCGGGGGGGTTTTCTTTTTACCCGTGAGGCgg                                                                                                                              >  1:433486/1‑148 (MQ=255)                           aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGGCggg                                                                                                                             >  1:784388/1‑148 (MQ=255)                             tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGCATTTGGGGCCCGGGGGGGGTTTTTTTTTACCAGGGGGACGGGGaa                                                                                                                          >  1:746777/1‑149 (MQ=255)                              cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                       <  2:502551/151‑1 (MQ=255)                                 aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGGGAGACGGGGAACAACTg                                                                                                                    >  2:63042/1‑151 (MQ=255)                                    aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                                  <  2:271375/150‑1 (MQ=255)                                       tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCAGGGGGGTTTTTTTTTTTACCAGGGAGACGGGCAACAGCTGATTGcc                                                                                                              >  2:839605/1‑151 (MQ=255)                                         cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                                                                                                            <  2:30551/151‑1 (MQ=255)                                                  cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                                                                                                     <  1:373447/149‑1 (MQ=255)                                                         cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGccc                                                                                              >  2:585533/1‑149 (MQ=255)                                                            tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                          >  1:368457/1‑150 (MQ=255)                                                               ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGATTTCCCTTTACCGCCTGGCCCCgag                                                                                          >  2:166680/1‑147 (MQ=255)                                                                cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGCGCCAGGGGGGTTTTTTTTTTCACCACTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCTTGCCCTGggagag                                                                                      >  2:535649/1‑150 (MQ=255)                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagcgaggt                                                                                    >  1:674445/1‑146 (MQ=255)                                                                    cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                 <  1:164786/151‑1 (MQ=255)                                                                     gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcag                                                                                <  1:911502/151‑1 (MQ=255)                                                                       gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                              <  2:387453/151‑1 (MQ=255)                                                                         aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGGGAGAGGTg                                                                                   >  2:192513/1‑144 (MQ=255)                                                                                  tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGCGCCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCCTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGTCCAc                                                                    >  2:873091/1‑150 (MQ=255)                                                                                    ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAAAAACTGATTGCCCTTTACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGc                                                                  >  2:411864/1‑150 (MQ=255)                                                                                     cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGGCAGGGGGGTTTTTCTTTTCACCAGTGAGACGGGGAACAGCTGGTTGCCCTTCACCGCCTGGCCCCGAGAGAGGTGGAGCAAGCGGTCCACGCTg                                                                >  1:32912/1‑151 (MQ=255)                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGCGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTGGAGCAAGCGGGCCACGCTGGGTTggc                                                         >  1:158552/1‑149 (MQ=255)                                                                                            aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcc                                                         >  2:773248/1‑151 (MQ=255)                                                                                             ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTGccc                                                        >  2:844931/1‑151 (MQ=255)                                                                                               aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGTTTTGGGGCCCGGGGGGGTTTTCTTTTCACCAGGGAGACGGGGAACAGCTGATTGCCCCTCCCCCCCCGGCCCTGGGAGAGTTGCCGCAAACGGGCCCCGCCGGTTTTcccc                                                       >  1:207476/1‑150 (MQ=255)                                                                                                aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTGcccc                                                       >  2:908246/1‑149 (MQ=255)                                                                                                  gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                   <  2:724694/151‑1 (MQ=255)                                                                                                   aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGCGGCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCCTCACCCCCCGGCCCTGGGAGAGTTGCAGCAAGCGGGCCCCGCCGGGTTTcccc                                                       >  2:770978/1‑146 (MQ=255)                                                                                                   aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                   <  1:14697/150‑1 (MQ=255)                                                                                                      cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                                <  2:546532/150‑1 (MQ=255)                                                                                                           aaCGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                          <  2:344614/151‑1 (MQ=255)                                                                                                                 cgGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                    <  1:962779/151‑1 (MQ=255)                                                                                                                         ggCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                            <  1:414938/151‑1 (MQ=255)                                                                                                                           cGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                          >  2:916648/1‑151 (MQ=255)                                                                                                                            ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                         >  1:346718/1‑151 (MQ=255)                                                                                                                             gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGATGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                         >  2:784030/1‑150 (MQ=255)                                                                                                                              tttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                        >  2:318608/1‑150 (MQ=255)                                                                                                                               ttGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                       <  1:192504/150‑1 (MQ=255)                                                                                                                               ttGCGTATTGGGCGCCAGGGAGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGACAATCCTGTTTGATGGTGGTTaa                       >  1:787490/1‑150 (MQ=255)                                                                                                                                tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                     >  2:53986/1‑151 (MQ=255)                                                                                                                                tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                     >  2:202868/1‑151 (MQ=255)                                                                                                                                     aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                 <  1:967222/150‑1 (MQ=255)                                                                                                                                     aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                 <  1:166673/150‑1 (MQ=255)                                                                                                                                     aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                >  1:663201/1‑151 (MQ=255)                                                                                                                                           cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa          <  2:43308/151‑1 (MQ=255)                                                                                                                                             ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa        <  2:784466/151‑1 (MQ=255)                                                                                                                                             ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa        <  1:369477/151‑1 (MQ=255)                                                                                                                                                   tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt  >  2:338886/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366380‑366676

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