Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,522	0	A	C	43.1%	64.3 / 13.4	55	L330R (CTG→CGG)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base C (22/0);  total (39/16)
Fisher's exact test for biased strand distribution p-value = 1.02e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670                                                                                                                                                    |                                                                                                                                                        aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGTA‑TTGGGGGCCAggggg                                                                                                                                                    >  2:710184/1‑149 (MQ=255) aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGggtg                                                                                                                                                    <  1:51546/151‑1 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTA‑TTTGGCGCCAGGGGGGgttt                                                                                                                                               >  1:155358/1‑151 (MQ=255)        ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTACAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttc                                                                                                                                             <  2:415572/151‑1 (MQ=255)          ggTGCCTAATGAGTGAGCTAACTCACATTGATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTTGGCGCCCggggggggtttttt                                                                                                                                           >  2:936422/1‑143 (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)                                                                                                                                                    |                                                                                                                                                        AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATT‑G‑CCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670

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