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

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366378‑366674                                                                                                                                                       |                                                                                                                                                         aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATAGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGG‑GTA‑TTTGGGGCCCgggggg                                                                                                                                                        >  2:848674/1‑151 (MQ=255)    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGA‑ATTGGCGCCAGGGGGGttt                                                                                                                                                     >  2:500176/1‑151 (MQ=255)        gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGCTTG‑GTA‑ATTGGGGCCAGGGGGGgtttttt                                                                                                                                                 >  2:3235/1‑149 (MQ=255)            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTT‑GTC‑TTTGGGGGCTGGGGGGGttttcttttt                                                                                                                                             >  1:1126695/1‑150 (MQ=255)            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                              <  2:1100764/150‑1 (MQ=255)              cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGTT‑TTGGGGGCCGGGGGGGTtttttttt                                                                                                                                               >  2:795819/1‑148 (MQ=255)                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAg                                                                                                                                        <  2:1066134/151‑1 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                       <  1:668194/150‑1 (MQ=255)                    agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                     <  2:447097/151‑1 (MQ=255)                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑ATTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTg                                                                                                                                      >  1:1118768/1‑149 (MQ=255)                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                   <  2:1091497/150‑1 (MQ=255)                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGC‑GGA‑GTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTgag                                                                                                                                    >  2:741050/1‑149 (MQ=255)                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGT‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGCgaga                                                                                                                                   >  1:388163/1‑149 (MQ=255)                          cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGC‑TTTGGCGCCCGGGGGGGTTTTTTTTTTTTCCGCGAGACggg                                                                                                                               >  1:505006/1‑151 (MQ=255)                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGc                                                                                                                              <  1:887298/150‑1 (MQ=255)                               tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGA‑TTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAA‑c                                                                                                                          >  2:177940/1‑150 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑Ca                                                                                                                         <  2:256146/150‑1 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑Ca                                                                                                                         <  2:1029979/150‑1 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑Ca                                                                                                                         <  2:1041305/150‑1 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑Ca                                                                                                                         <  2:161768/150‑1 (MQ=255)                                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCaa                                                                                                                            >  1:566639/1‑144 (MQ=255)                                      aaTTGCGTTGCGCTCACTGCCGGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTAGC‑GGG‑TTAGGGGCCAGGGTGGGTTTTATTTTTTCCACTGAGACGGGcggaagctggt                                                                                                                    >  2:97010/1‑140 (MQ=255)                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑TTT‑TTTGGGGCCAGGGGGGGTTTTTTTTTCATCCATGAGACGGGCAA‑Cacgggat                                                                                                                   >  1:736160/1‑141 (MQ=255)                                         tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTG‑GTC‑TTTGGCGCCAGGGGGGTTTTTCTTTTCACCCCTGCGACGGGGAA‑CCG‑ccggtt                                                                                                                 >  2:564274/1‑143 (MQ=255)                                           cgttgcgCTCACTACCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTA‑GTA‑TTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGCGA‑AAT‑CGGATTGccc                                                                                                             >  1:471936/1‑150 (MQ=255)                                            gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCtt                                                                                                           <  1:911046/151‑1 (MQ=255)                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGTT‑TTGGGGGCCGGGGGGGTTTTT‑TTTTTACCAGGGGGAGGGGGAA‑AAG‑CTGAGTGCCCtt                                                                                                           >  2:840739/1‑149 (MQ=255)                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTGGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAA‑CAG‑CTGATTTcccct                                                                                                           >  1:385921/1‑147 (MQ=255)                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTT‑TTTGGGGCCAGGGGGGGGTTTCTTTTTACCAATGAGAAGGGGAA‑CAG‑CTGAATGCCCTTCA‑CCGcc                                                                                                   >  2:811775/1‑150 (MQ=255)                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGca                                                                                                   <  2:491078/150‑2 (MQ=255)                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTg                                                                                                 <  2:842763/151‑1 (MQ=255)                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTTGGCGCCCCGGGGGGTTTTTTTTTTACCCCTGGGACGGGCCA‑CAG‑CTTATTGCCCCTCACCCGCCGGGccccg                                                                                           >  1:341712/1‑149 (MQ=255)                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAATGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCCCCCGGCCCTgg                                                                                          >  1:989440/1‑150 (MQ=255)                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGt                                                                                    <  1:347548/150‑1 (MQ=255)                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGT‑GGT‑TTTGGCGGCAGGGGGGGTTTTTTTTTTCCCACGGAGACGGGGAA‑CCGCCGGATTTCCCTCCC‑CCCCCTGCCCTTggggag                                                                                      >  2:1060074/1‑137 (MQ=255)                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTG‑TTAGGCGGCAGGGGGGGTTTTTTTTTTTCCCGTGAAACGGGCAA‑CAG‑CTGATTGTCCCTCC‑CCCCCTGGCCCCGAGGGAGTTgcag                                                                               >  1:846957/1‑147 (MQ=255)                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCgg                                                                        <  1:903829/149‑1 (MQ=255)                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑Tc                                                                     <  2:989415/151‑1 (MQ=255)                                                                                 tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGCGGCAGGGGGGGTTTTTTTTTTACCCGGGAGACGGGCAA‑AAG‑CTGATTTGCCCTCC‑CCCCCTGGCCCCGGGGGAGTTGCAGCAAACGG‑Gcc                                                                    >  2:187382/1‑151 (MQ=255)                                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTg                                                              <  2:702690/151‑1 (MQ=255)                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGt                                                            <  1:48016/151‑1 (MQ=255)                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GTT‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAA‑CAA‑CCGATTGGCCCTCA‑CCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGgtt                                                          >  1:40381/1‑151 (MQ=255)                                                                                               ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑ATTGGGGCCACGGTGGGTTTTCTTTTCACCAGTGAGACGGGCAACCAC‑CTG‑TTGGCCTTCC‑CCGCCCTGCCCCGGGAGAGGTGCCGCAAGCGG‑Tcccacctgggtttccc                                                      >  2:344118/1‑137 (MQ=255)                                                                                                tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGA‑TTTGGCGCCAGGGGGGTTTTTCTTTTTACCCGTGAGACGGGGAA‑CAG‑CTGATTGGCCTTCA‑CCGCCCGGCCCTGAGAGAGTTGCAACCAGCGG‑GCCACCCTGGGTTGcccc                                                     >  1:1109513/1‑151 (MQ=255)                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGATTGC‑GCA‑GAGGGGGGCCGGGTGGGTTTTCTTTTCTCCCGCGAGAAGGGCAA‑CAA‑CTAATTGGCCTTCA‑CCCCCCGGCCCTGAGAGAGTTTGAGCAAGCGG‑Tcccacccgggttgcc                                                       >  2:422958/1‑131 (MQ=255)                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGc                                              <  1:773004/150‑1 (MQ=255)                                                                                                          gCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGa                                            <  2:250910/150‑1 (MQ=255)                                                                                                           ccAACGCGCGGGGAGAGGCGGTTTGT‑GTA‑TTGGGCGCCCGGGGGGTTTTTCTTTTCACCAGAGAGACGGGCAA‑CAG‑CTGATTGCCCCTCA‑CCGCCCGGCCCCGGGAGAGTTGCAGCAAGCGG‑TCCACGCCGGGTTGTCCCCGCAAGCg                                             >  1:112580/1‑148 (MQ=255)                                                                                                            cAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAGGGGCAA‑CCG‑CCGGTTGCCCTTCA‑CCCCCCGGCCCTGGGAGAGGTGCAGCAAACGG‑GCCACGCCGGGTTGCCCCCGCAGGCgga                                           >  2:498278/1‑147 (MQ=255)                                                                                                            cAACGCGCGGGGAGAGGCGGTTTGT‑GTA‑TTGGGCGCCCGGGGGGGTTTTCTTTTCACCCGGGAGACGGGCAA‑AAG‑CGGGTTGCCCTTCA‑CCCCCCGGCCCCGAGAGAGTTGCAGCAAGCGG‑GCCACGCTGGTTTTCCCccgcaagggca                                           >  2:945100/1‑143 (MQ=255)                                                                                                             aaCGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCAATGAGACGGGGCA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGGGGGGTGGAGCAAGGGG‑TCCACGCCGGGTTTCCCCCGCCGGGGgaaa                                         >  1:304018/1‑150 (MQ=255)                                                                                                              aCGCGCGGGGAGAGGCGGTTTGG‑GTA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTTTCAGTTGGACGGGGAA‑AAT‑GTGAATGTCCCTCA‑CCGCCTTGCCCTAAGAGAGTGGTAGCAAGCGGCTCTAC‑CTGGATTTCCCCATCAGGCGaaaa                                         >  1:235269/1‑149 (MQ=255)                                                                                                                gcgcgGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCt                                     <  1:429263/151‑1 (MQ=255)                                                                                                                  gcgGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGGCAGGGTGGTTTTTCTTTTCACCAGGGAGGCGGGGAA‑CAA‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGGGAGAGGTGCAGCAAGCGG‑GCCCCGCTGGGTTGCCCCCACCGGCGaaacac                                       >  2:774046/1‑144 (MQ=255)                                                                                                                    ggggAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                  <  2:137655/150‑1 (MQ=255)                                                                                                                       gagaGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                              <  2:417617/151‑1 (MQ=255)                                                                                                                          agGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                            >  2:891781/1‑150 (MQ=255)                                                                                                                            gCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                         >  1:1047697/1‑151 (MQ=255)                                                                                                                               gTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  1:498293/151‑1 (MQ=255)                                                                                                                                    c‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  <  2:736139/150‑1 (MQ=255)                                                                                                                                        a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               <  2:652364/150‑1 (MQ=255)                                                                                                                                        a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               <  2:688471/150‑1 (MQ=255)                                                                                                                                          ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCgggg             >  2:1102404/1‑150 (MQ=255)                                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          >  1:1085897/1‑149 (MQ=255)                                                                                                                                                 ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa        >  1:15199/1‑149 (MQ=255)                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa   <  1:146275/151‑1 (MQ=255)                                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  1:811800/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAA‑CAG‑CTGATTGCCCTTCA‑CCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGG‑TCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366378‑366674

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