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

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAC  >  NC_000913/366377‑366669                                                                                                                                                       |                                                                                                                                                 aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGggt                                                                                                                                                  <  1:1356423/150‑1 (MQ=255)  aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGggtgg                                                                                                                                                <  1:960865/151‑1 (MQ=255)    gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGtt                                                                                                                                              <  1:586855/151‑1 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGT‑GTTTGGCGCCAGGGGGGgttt                                                                                                                                            >  1:870617/1‑151 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGCGCCAGGGGGGgttt                                                                                                                                            >  1:220479/1‑151 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGGGGCCGGGGGGgttt                                                                                                                                            >  1:309678/1‑151 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGTGG‑ATTTGGGGCCAGGGGGGgttt                                                                                                                                            >  2:203490/1‑151 (MQ=255)         gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑ATTTGGCGCCCGGGGGGgtttttt                                                                                                                                         >  1:1371591/1‑149 (MQ=255)          ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTttttctt                                                                                                                                        <  1:1137232/151‑1 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                      <  1:654475/150‑1 (MQ=255)               cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑ATTGGGCGGCCGGGGGGGttttcttt                                                                                                                                       >  2:978222/1‑147 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTCGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                             <  1:98128/151‑1 (MQ=255)                        gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                          <  2:1320213/151‑1 (MQ=255)                         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGCGCCAGGGGGGGtttttttttt                                                                                                                                     >  1:1266997/1‑138 (MQ=255)                         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                         <  2:403959/151‑1 (MQ=255)                            tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑ATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACggg                                                                                                                       >  1:1044334/1‑150 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑ATTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACggg                                                                                                                       >  1:676693/1‑149 (MQ=255)                                   caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑ATTTGGCGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACaac                                                                                                                >  1:336557/1‑148 (MQ=255)                                    aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGG‑ATTTGGGGCCAGGGGGGGTTTTTTTTTCACCCGTGAGAAGGGGAACAGCTg                                                                                                              >  1:1013457/1‑151 (MQ=255)                                       aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                           <  1:147220/151‑1 (MQ=255)                                       aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                           <  1:237276/151‑1 (MQ=255)                                         ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGAtt                                                                                                           >  1:1256167/1‑149 (MQ=255)                                         ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAACAGCTGAtt                                                                                                           >  1:1298725/1‑149 (MQ=255)                                               tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACACGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTc                                                                                                   <  1:1038698/150‑1 (MQ=255)                                                 cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑ATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTcccc                                                                                                >  1:1362231/1‑151 (MQ=255)                                                    tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑GTTTGGGGCCCGGGGGGGGTTTTTTTTTTCCCGGGAGACGGGGACCAGGCGAATGC‑CCCTCACcccc                                                                                             >  2:651139/1‑151 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGT‑CTTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGCGAGAAGGGCAACAACCGATTGC‑CCTTCACcccc                                                                                             >  1:651604/1‑150 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGc                                                                                              <  1:73343/149‑1 (MQ=255)                                                       cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑TTTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGGGAGGCGGGCAACAGCTGATTGG‑CCtttcacccc                                                                                             >  1:1330862/1‑141 (MQ=255)                                                          cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGccc                                                                                       <  1:605208/151‑1 (MQ=255)                                                             gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGT‑ATTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCCCCCGGccc                                                                                       >  1:672213/1‑148 (MQ=255)                                                                   cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑GTTGGGGGCCCGGGGGGGTTTTTTTTTTCCCCGCGAGAAGGGGGACCGCGGGTTGT‑TCTTTACCGCCTgg                                                                                          >  2:310502/1‑139 (MQ=255)                                                                      tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGT‑TTTGGGGGCCCGGGGGGGTTTTCTTTTTACCAGCGAGAAGGGCAACAACTGATTGG‑CCCTCTCCCCCCGGCCCCGGGAGAGGTgc                                                                           >  1:814658/1‑151 (MQ=255)                                                                       cGGGAAACCGGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACTGCCTGGCCCTGAGAGAGTTgc                                                                           <  2:298217/150‑1 (MQ=255)                                                                           aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTTC‑CTTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                      <  2:1368958/151‑1 (MQ=255)                                                                             aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                    <  1:311774/151‑1 (MQ=255)                                                                                 gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGGGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAACTGATTGC‑CCTTCACCGCCCGGCCCTGAGAGAGGTGGAGCAAACgg                                                                  >  1:1321954/1‑149 (MQ=255)                                                                                   cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTcc                                                               <  1:1063717/150‑1 (MQ=255)                                                                                    gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTTGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCCTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGTCCAc                                                             >  1:352619/1‑151 (MQ=255)                                                                                     tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑ATTGGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAAACGCTGCTTGC‑CCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAgggg                                                                  >  2:47233/1‑145 (MQ=255)                                                                                      gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGT‑ATTTGGCGCCCGGGGGGGTTTTTTTTTCTCCCGTGAGGAGGGGAACAACCGATTTC‑CCTTTCCCCCCTGGCCCTTAGAGAGTTGCCACAAGAgggcacc                                                             >  2:1149168/1‑144 (MQ=255)                                                                                       ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                          <  1:1302949/151‑1 (MQ=255)                                                                                          gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                       <  1:525483/151‑1 (MQ=255)                                                                                           cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGT‑ATTTGGCGCCCGGGGGGGGTTTTTTTTTACCCGGGAGGCGGGGAACAGCCGGTTGG‑CCCTCACCCCCTGGGCCCGAGGGAGGTGGAGGAAGCGGTCCACGCTgg                                                        >  2:532194/1‑149 (MQ=255)                                                                                           cTGCATTAATGAATCGGCCAAAGCGCGGGGAGAGGCGGCTTTCGTGTTTGGGGCCCAGGGTGG‑TTTTCTTTTTACAAATGCGACGGGCAACACCCGATTGCGCCTTC‑CCGCCTGGCCCCGCGAGAGTTgtaggacaaggtgcccac                                                           >  2:1401168/1‑132 (MQ=255)                                                                                              cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGtt                                                      <  1:809328/148‑1 (MQ=255)                                                                                              cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                   <  1:783247/151‑1 (MQ=255)                                                                                                 tAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTGT‑ATTTGGGGGCAGGGGGGGTTTTTTTTTCACCAGAGAGACGGGGAAAAACCGATTTC‑CCTTCCCCGCCCTGCCCCGAGAGAGGTGCCGCAAGCGGGTCACGCCGGGTTTccc                                                 >  1:1091776/1‑150 (MQ=255)                                                                                                  aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGcccc                                                <  2:1104700/150‑1 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGG‑ATTTGGCGCCAGGGGGGGTTTTCTTTTCACCAATGAGAAGGGGAACAACCTATTTC‑CCCTCACCCTCTTGCCCTGAGAGAGTTTGAGCAAACGGTCCACCCTGGGTTGCCCcagc                                             >  1:1154766/1‑150 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                            <  1:147262/151‑1 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                            <  1:995122/151‑1 (MQ=255)                                                                                                      aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATT‑GGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                          <  2:1287565/151‑1 (MQ=255)                                                                                                       aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑ATTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGAAGGGGAACAGCTGATTTC‑CCTTCACCGCCCGGCCCCGAGAGAGTTGCAGCAAGCGGGCCACGCTGGGTTGCCCCCGCAgg                                          >  1:256019/1‑151 (MQ=255)                                                                                                         cGGCCAACGCGCGGGGAGAGGCGGTTTGCGG‑ATTGGGCGCCCGGGGGGGTTTTTTTTTTCCCCGGGAGAAGGGGAAAAGCCGATTGC‑CCTTCACCGCCTGGGCCTTAGGGAGGTGGAGCAAGCGGGCCACGCCGGGTTGgccccggcgggg                                        >  2:252459/1‑143 (MQ=255)                                                                                                             cAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGataa                                    <  2:1062930/151‑1 (MQ=255)                                                                                                              aaCGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                   <  2:670465/151‑1 (MQ=255)                                                                                                               aCGCGCGGGGAGAGGCGGTTAGCGT‑ATTGGGCGCCAGGTTGCGTTTTCTTTTCGCCAGTAAGCGGGGCAGCAGCTGTTTGC‑CGTGCACCGCCAGGCCCTGAGAGGGTTGCAGCAAGCGGTCGACGGTGGTGTGCCCCAGCAGGCGaaa                                     <  2:337502/148‑1 (MQ=255)                                                                                                                     gttgagaGGCCGTTTGCGT‑TTTGGGCGCTAGGGTGGTCTTTCTTTTCCCCACTGAGACGGGCAACAGTTGATTAC‑CCTTTACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACTCTGGTTTTCCCCAGCAGGCGAAAATCCTGt                              <  2:138693/146‑1 (MQ=255)                                                                                                                       ggAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                          <  2:1202868/151‑1 (MQ=255)                                                                                                                       ggAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                          <  1:532826/151‑1 (MQ=255)                                                                                                                           agGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                      >  1:592379/1‑151 (MQ=255)                                                                                                                                gTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                 <  1:1339670/151‑1 (MQ=255)                                                                                                                                  ttGCGT‑ATTGGGCGCCCGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc               >  2:1074897/1‑151 (MQ=255)                                                                                                                                  ttGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc               >  1:311484/1‑151 (MQ=255)                                                                                                                                  ttGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc               >  1:629689/1‑151 (MQ=255)                                                                                                                                   tGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc               <  1:922421/150‑1 (MQ=255)                                                                                                                                    gCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg             <  2:673134/151‑1 (MQ=255)                                                                                                                                     cGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc            <  1:66064/151‑1 (MQ=255)                                                                                                                                       t‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACACTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg          >  1:191655/1‑151 (MQ=255)                                                                                                                                            gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata      <  1:1075997/151‑1 (MQ=255)                                                                                                                                            gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat       <  1:739377/150‑1 (MQ=255)                                                                                                                                             ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata      >  2:1402916/1‑150 (MQ=255)                                                                                                                                              gcgcCAGGGTGGTTTTTCTTTCCATCAGCGAGACGGGCAACAGCTGACTGC‑CCTTCACCTCCTGGCCCTGAGGGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGACGGTGGTTAACGGCGGGatat     >  2:84349/1‑150 (MQ=255)                                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat       >  2:524721/1‑147 (MQ=255)                                                                                                                                                 ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaa  >  2:367280/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                 AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGT‑ATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGC‑CCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAC  >  NC_000913/366377‑366669

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