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

CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366381‑366674                                                                                                                                                   |                                                                                                                                                      ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCGAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGttt                                                                                                                                                  <  1:940358/151‑1 (MQ=255)   tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGGGCCCGGGGGGgtttt                                                                                                                                                >  1:6652/1‑151 (MQ=255)    ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GTATTGGGGGGCCGGGGGGgtttt                                                                                                                                                >  2:595113/1‑150 (MQ=255)    ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGtttt                                                                                                                                                 <  2:351473/149‑1 (MQ=255)         gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTG‑GGATTTGGCGCCCGGGGGGGtttttttt                                                                                                                                            >  2:947113/1‑149 (MQ=255)         gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                           <  2:300031/150‑1 (MQ=255)            tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGCCAGGGGGGGTTTTTTTTTTAcc                                                                                                                                       >  2:393505/1‑151 (MQ=255)             aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGGCAGGGGGGGTTTTCTTTTcccc                                                                                                                                       >  2:965548/1‑150 (MQ=255)               tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGC‑GGGTTGGGGGCCAGGGGGGGTTTTCTTTTCAcc                                                                                                                                       >  2:374806/1‑148 (MQ=255)                gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                    <  1:552059/150‑1 (MQ=255)                  gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                 <  2:452786/151‑1 (MQ=255)                   tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGATTTGGCGCCCGGGGGGGttttctttttccc                                                                                                                                       >  2:856055/1‑140 (MQ=255)                     agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAc                                                                                                                               >  1:222760/1‑150 (MQ=255)                       cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGATTTGGGGCCCGGGGGGGTTTTTTTTTTAACAGCGAGGCggg                                                                                                                            >  2:404616/1‑151 (MQ=255)                            tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGC‑GGATTTGGCGCCCGGGGGGGTTTTTTTTTCCCCAGGGGGACggg                                                                                                                            >  2:55146/1‑146 (MQ=255)                             cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GTATTTGGGGGCAGGGGGGGTTTTTTTTTTACCCGTGAAACggggca                                                                                                                         >  2:920344/1‑145 (MQ=255)                             cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                                                                                                      <  1:539003/151‑1 (MQ=255)                                aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAAAAGCTg                                                                                                                   >  1:538712/1‑151 (MQ=255)                                   aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                                                                                                 <  2:423404/150‑1 (MQ=255)                                      tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGGCCGGGTGGTTTTTCTTTTCACCAGGGAGAAgggggacaaaacgtt                                                                                                                >  2:178270/1‑136 (MQ=255)                                       gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCCGGGGGGTTTTTTTTTTTACCAATGAGACGGGCAACAGCTGATTTccc                                                                                                            >  1:674013/1‑151 (MQ=255)                                            gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                                                                                                       <  2:186837/151‑1 (MQ=255)                                              gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGTTTTGGGGCCCGGGGGGGTTTTTTTTTTCCCCA‑TGAGACGGGGCACAACTGATTTGCCCTCAccc                                                                                                     >  2:837190/1‑150 (MQ=255)                                              gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                                                                                                     <  2:597837/151‑1 (MQ=255)                                                  aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGATTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCTTCACcccc                                                                                                   >  1:894350/1‑149 (MQ=255)                                                   cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGCACAGCTGATTGCCCTTCACcccc                                                                                                   >  2:1038640/1‑148 (MQ=255)                                                            ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                         <  1:980423/149‑1 (MQ=255)                                                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGGGAGAAGGGGAACAGCCGATTTCCCCTCACCCCCCGGccccggg                                                                                         >  2:414463/1‑143 (MQ=255)                                                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTG‑GTATTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGGGAGAAGGGCAAAAGCTGATTGCCCCTTACCGCCCGGGCCCgag                                                                                         >  2:609781/1‑147 (MQ=255)                                                                     ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGCGGCAGGGGGGGTTTTTTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGTTgcag                                                                               >  1:687239/1‑150 (MQ=255)                                                                     ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGC‑GTATTGGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACGGGTGATTGCCCCTCCCCCCCTGGCCCCgagaggg                                                                                     >  2:473388/1‑142 (MQ=255)                                                                      gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  1:51092/151‑1 (MQ=255)                                                                       aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                            <  2:399168/151‑1 (MQ=255)                                                                          ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCg                                                                         <  2:878693/151‑1 (MQ=255)                                                                            tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGGGCCCGGGGGGGTTTTTTTTTCACCCGGGAGAAGGGCAACAGCTGGTTGCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAACg                                                                         >  2:549074/1‑149 (MQ=255)                                                                             gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                        <  2:453702/149‑1 (MQ=255)                                                                               cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGATTTGGCGCCACGGGGGGTTTTCTTTTTACCAGCGAGACGGGCAACAGCTGATTGCCCCTCACCCCCCCGCCCCCTGAGGGGTGCAGCAAGCGGGcc                                                                     >  2:73349/1‑150 (MQ=255)                                                                               cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                                     <  2:910786/150‑1 (MQ=255)                                                                               cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                    <  1:800712/151‑1 (MQ=255)                                                                               cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGGTTTTCTTTTTACCAGGGGGACGGGCAACAGCTGATTTGCCTTTCCCGCCCGGCCCTGGGAGAGGTGCAGCAAACGGGcc                                                                     >  2:620082/1‑150 (MQ=255)                                                                                  gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCAGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGc                                                                 <  1:554937/151‑1 (MQ=255)                                                                                    cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑g                                                              <  1:861895/151‑1 (MQ=255)                                                                                        tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGATTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCCACAGCTGATTGCCCCTCCCCCCCCGGCCCCGAGAGAGGTTCAGCAAACGGGCCCCGCT‑GGttt                                                          >  2:145395/1‑151 (MQ=255)                                                                                            ttAATGAATCGGCAATCGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCTTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGttt                                                          <  2:932868/147‑1 (MQ=255)                                                                                              aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGGCAGGGGGGTTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTTCCCTTCACCCCCTGGCCCTGGGAGAGTTGCCGCAAGCGGGCCACCCT‑GGGTTTcccc                                                     >  1:329845/1‑150 (MQ=255)                                                                                               aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTATTGGGCGCCCGGGGGGGTTTTCTTTTCCCCAAGGGGACGGGCAACAGCTGATTGCCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACCCTGGGTTTGccc                                                      >  2:845606/1‑149 (MQ=255)                                                                                                 gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCcagca                                                 <  1:841331/151‑1 (MQ=255)                                                                                                 gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCcagc                                                  <  1:711794/150‑1 (MQ=255)                                                                                                 gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCcagc                                                  <  1:809020/150‑1 (MQ=255)                                                                                                  aaTCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTATTTGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAATGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGGTTGCCCccgc                                                  >  1:777370/1‑149 (MQ=255)                                                                                                   aTCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAgg                                               <  1:613163/151‑1 (MQ=255)                                                                                                     cGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTATTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCg                                             >  1:309531/1‑151 (MQ=255)                                                                                                     cGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGGTTGCCCCAGCAgg                                               >  1:474953/1‑149 (MQ=255)                                                                                                      ggCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGCGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGGAACAGCTGATTGGCCCTCACCCCCCGGCCCCGGGAGAGTTGGAGCAAAAGGGCCACGCC‑GGTTTGCCCCCGCAcgcg                                             >  2:810360/1‑150 (MQ=255)                                                                                                         cAACGCGCGGGGAGAGGCGGTTTTC‑GTATTGGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCCACAGCTGATTGCCCTTCACCGCCCGGCCCTGGGAGAGGTGCAGCAAGCGGGCCACGCT‑GGTTTGCCCCCGCAGGCGaaa                                          >  2:129416/1‑150 (MQ=255)                                                                                                               gcgGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGt                                   <  1:36238/151‑1 (MQ=255)                                                                                                                       agGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                            >  2:344246/1‑150 (MQ=255)                                                                                                                            gTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  2:481318/151‑1 (MQ=255)                                                                                                                             tttGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      >  2:721598/1‑150 (MQ=255)                                                                                                                                gC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  <  1:1068890/151‑1 (MQ=255)                                                                                                                                   gTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                <  1:1020206/151‑1 (MQ=255)                                                                                                                                   gTATTGGGCGCCAGGGTGGTTTTTCCTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                >  2:219287/1‑151 (MQ=255)                                                                                                                                    tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCAACGCT‑GGATTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               >  2:242451/1‑151 (MQ=255)                                                                                                                                     aTTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              <  2:373609/151‑1 (MQ=255)                                                                                                                                      ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCgggg             >  2:578241/1‑150 (MQ=255)                                                                                                                                      ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              <  1:779007/150‑1 (MQ=255)                                                                                                                                           cgcCAGGGTGG‑TTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa        >  1:1004113/1‑150 (MQ=255)                                                                                                                                               tgggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      >  2:316362/2‑149 (MQ=255)                                                                                                                                               aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  1:1008931/151‑1 (MQ=255)                                                                                                                                                 ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  1:466707/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCT‑GGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366381‑366674

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