Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,534	0	C	T	42.0%	65.7 / 18.4	53	R326K (AGA→AAA)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base T (21/0);  total (38/15)
Fisher's exact test for biased strand distribution p-value = 3.60e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGT  >  NC_000913/366384‑366678                                                                                                                                                         |                                                                                                                                                    ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTACAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttc                                                                                                                                                     <  2:415572/151‑1 (MQ=255)     tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTGGGCGCCCGGGGGGGttttttttt                                                                                                                                                 >  2:697953/1‑151 (MQ=255)      gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                                 <  2:27743/150‑1 (MQ=255)      gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                                 <  2:30992/150‑1 (MQ=255)       ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGCATTTTGGGGCCGGGGGGGTtttttttt                                                                                                                                                  >  2:259185/1‑149 (MQ=255)          aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                            <  1:896561/151‑1 (MQ=255)           aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                             >  1:296935/1‑149 (MQ=255)                       aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCaa                                                                                                                              <  2:341218/151‑1 (MQ=255)                         tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGGA‑TTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGAGA‑CGGGCaa                                                                                                                              >  2:859052/1‑149 (MQ=255)                              ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCGTGAGA‑CGGGCAACa                                                                                                                            >  2:104448/1‑146 (MQ=255)                              ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGA‑CGGGCAAc                                                                                                                             >  2:759837/1‑145 (MQ=255)                                  ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGAACAGCTGAtt                                                                                                                     >  2:353014/1‑149 (MQ=255)                                      gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGA‑TTGGGGGGCCGGGGGGGTTTTTTTTTTTCCAGGGAGA‑AGGGGGACAGCTGATT‑G‑tcc                                                                                                               >  1:748805/1‑149 (MQ=255)                                       ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGA‑CGGGGAACAACTGATT‑G‑cccct                                                                                                             >  1:567202/1‑148 (MQ=255)                                              cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGc                                                                                                       <  2:885320/149‑1 (MQ=255)                                               aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGGGGT‑TTTGGCGCCCGGGGGGGTTTTCTTTTTCCCCGTGGGA‑CGGGCAACAGCTGATT‑G‑CCCCTCACcccc                                                                                                      >  1:195724/1‑149 (MQ=255)                                                 tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTG‑TTGGGGGGCAGGGGGGGTTTTTTTTTTACCAGGGAGG‑AGGGGAACAGCTGGTT‑G‑CCCTTCACcccc                                                                                                      >  1:764185/1‑147 (MQ=255)                                                  gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTCACCCAGGAGA‑CGGGGCACAGCCGGTTGG‑CCCTTCccccccgg                                                                                                    >  1:486255/1‑147 (MQ=255)                                                    ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGTG‑TTGGGGGGCAGGGGGGGTTTTTTTTTTCCCCGGGAGA‑CGGGCAAAACCCGATT‑G‑CCCCTCACCCCCCGGccc                                                                                                >  2:781833/1‑150 (MQ=255)                                                    ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAGA‑CGGGCAACAGCTGATT‑T‑CCCTTCACCCCCCTGccc                                                                                                >  2:224325/1‑150 (MQ=255)                                                    ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGCAAAAGCTGATT‑G‑CCCTTCACCCCCTGGccc                                                                                                >  2:221220/1‑150 (MQ=255)                                                    ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCt                                                                                               <  2:106325/151‑1 (MQ=255)                                                      gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGGACAGCTTGTT‑G‑CCCCTCCCCCCCCGGccc                                                                                                >  1:915757/1‑148 (MQ=255)                                                                cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGG‑CGGGCAACAACTGATT‑G‑CCCTTCACCGCCCGGCCCCGagagaggt                                                                                      >  2:878026/1‑146 (MQ=255)                                                                cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGCCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGA‑CGGGGCACAGCTGATT‑G‑CCCTTCACCCCCCGGCCCTGGGAGAGGTg                                                                                     >  2:713345/1‑149 (MQ=255)                                                                 gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑AGGGGAACAACTGATT‑G‑CCCTTCACCCCCCGGCCCCGagagaggt                                                                                      >  2:244512/1‑145 (MQ=255)                                                                    aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCCGGGGGGTTTTTTTTTTTACCAGTGAGA‑CGGGCAAAAGCTGATT‑T‑CCCTTCACCCCCTGGCCCTGAGAGAGGTgcagca                                                                                >  2:1116175/1‑150 (MQ=255)                                                                     aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGCAAAAGCGGATT‑G‑CCCTTCCCCCCCCGGCCCTGAGAGAGGTGCAGCaa                                                                               >  2:123952/1‑150 (MQ=255)                                                                      aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTTACCCGTGAGA‑CGGGGAACAGCTGATT‑G‑CCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCaa                                                                               >  2:258877/1‑149 (MQ=255)                                                                         tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTTCTAT‑TTTGGCGCCGGGGGGGTTTTT‑TTTTTACCCGTGAGA‑CGTGCACCAGCACCTT‑GACCCTCCACC‑CCTGTTCCTCagagaggtttatcaatccg                                                                           >  1:209794/1‑137 (MQ=255)                                                                         tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGA‑AGGGCAACAGCTGATT‑G‑GCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAGCggg                                                                          >  1:1025718/1‑150 (MQ=255)                                                                          gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTc                                                                         <  1:218461/151‑1 (MQ=255)                                                                            cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGCT‑TTGGGGGCCGGGGGGGTTTTT‑TTTTTACCCCTGGGA‑CGGGGCACAGCCTAAT‑G‑TCCTTCACCCCCCCGCCCCGTGAGAGtgggcgccagccggcc                                                                        >  1:280624/1‑135 (MQ=255)                                                                            cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                       <  1:136403/151‑1 (MQ=255)                                                                                 cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCCTGAGA‑CGGGGAACAGCTGATT‑T‑CCCCTCCCCCCCCGGCCCTGAGAGAGGTGCAGCCAACGGTCCAc                                                                      >  1:563343/1‑147 (MQ=255)                                                                                       cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                            <  1:598537/151‑1 (MQ=255)                                                                                        aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGGG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGGGA‑CAGGCAAAAACCGATT‑G‑CCCTTCCCCGCCAGGCCCTCGGAGAGGTGCAGCAAGAGGGCCACGCCGGgtt                                                              >  2:760125/1‑148 (MQ=255)                                                                                        aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGGAAAAGCTGATT‑T‑GCCTTCACCGCCTGGGCCTGAGAGAGTTGCAGCAAACGGGCCACGCCGGGTTTcc                                                           >  2:740777/1‑151 (MQ=255)                                                                                          tAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTCCCCGTGAGA‑CGGGCAACAACCGATT‑G‑CCCCTTACCGCCCGGCCCTGAGAGAGGTGGAGCAAACGGGCCACGCCGGGTTTcccc                                                         >  2:1015672/1‑151 (MQ=255)                                                                                            aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                       <  2:395138/151‑1 (MQ=255)                                                                                              gAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                     <  1:416480/151‑1 (MQ=255)                                                                                               aaTCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                    <  1:1152805/151‑1 (MQ=255)                                                                                                   ggCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                                <  2:995191/151‑1 (MQ=255)                                                                                                   ggCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                                <  1:1119601/151‑1 (MQ=255)                                                                                                     ccAACGCGCGGGGAGAGGCGG‑TTTGGCTC‑ATAGGGGCCCGGGGGGTTTTTCTTTTCATCAGTGAGACCGGGCAA‑AGCGGACT‑G‑GCCTTTTCCGCCTTGCCCCGCGGGAGGTGCAGCAAGCCGGTCACCCTGGGGTGCCCCAGCAGGcaca                                               >  2:611104/1‑147 (MQ=255)                                                                                                                        gg‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                           >  2:560599/1‑151 (MQ=255)                                                                                                                           tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                          >  1:1105291/1‑150 (MQ=255)                                                                                                                                  a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                  >  2:523826/1‑151 (MQ=255)                                                                                                                                         cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg                  >  2:577047/1‑145 (MQ=255)                                                                                                                                         cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata             >  1:155657/1‑150 (MQ=255)                                                                                                                                            cAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa          >  1:803836/1‑148 (MQ=255)                                                                                                                                             aGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaaaa          >  1:524799/1‑147 (MQ=255)                                                                                                                                               ggtggtTTTTCTTTTCACCATTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTGCACCGCTTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg        <  2:751533/149‑1 (MQ=255)                                                                                                                                               ggtggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg      <  1:581916/151‑1 (MQ=255)                                                                                                                                                gtggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATAACATGa       >  1:333607/1‑149 (MQ=255)                                                                                                                                                gtggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg        >  1:1117874/1‑148 (MQ=255)                                                                                                                                                  ggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt    <  1:759800/150‑1 (MQ=255)                                                                                                                                                   gtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGt  >  2:307292/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGT  >  NC_000913/366384‑366678

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: ‑