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

CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366381‑366673                                                                                                                                                     |                                                                                                                                                      ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGGGCC‑A‑GGGGGGgtt                                                                                                                                                  >  2:177612/1‑151 (MQ=255)      ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTCGCGCTGACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATTGGCCAATTCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑C‑GGGTGGTttttctt                                                                                                                                             <  2:593595/151‑1 (MQ=255)         gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGCGTCTTGGGGGCCAA‑GGGGGGttttttt                                                                                                                                              >  2:269200/1‑146 (MQ=255)             aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGATTTGGGGCC‑C‑GGGGGGGTTTTCTTTTT‑Acc                                                                                                                                      >  2:528011/1‑150 (MQ=255)               tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGATTGGGCGCC‑C‑CGGGGGGTTTTTTTTTC‑Acc                                                                                                                                      >  2:441050/1‑148 (MQ=255)                 agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTTCGTCTTGGGGCCC‑C‑GGGGGGTGTTTTTTTTT‑ACCAGTg                                                                                                                                  >  2:39300/1‑150 (MQ=255)                    gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTCTTTGGGGCC‑A‑CGGGGGGTTTTCTTTTT‑ACCAGGgaga                                                                                                                               >  2:354486/1‑150 (MQ=255)                    gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTgaga                                                                                                                               <  1:910510/150‑1 (MQ=255)                      gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACgg                                                                                                                            <  1:621488/151‑1 (MQ=255)                        tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTTGGCGCC‑A‑GGGGGGGTTTTTTTTTT‑ACCAGCGAGACggg                                                                                                                           >  1:460241/1‑150 (MQ=255)                          aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCaa                                                                                                                        <  1:712195/151‑1 (MQ=255)                             cacaTTAATTGCGTTGGGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGCGTCTTGGGCGGC‑G‑GGGGGGGTTTTTTTTTT‑CCCACAGAGACGGGcca                                                                                                                        >  1:272655/1‑146 (MQ=255)                                aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGATTTGGGGCC‑C‑GGGGGGTTTTTTTTTTT‑ACCCGTGAGACGGGGCACA‑GCTg                                                                                                                 >  1:690334/1‑151 (MQ=255)                                 ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGATTTGGCGCC‑C‑GGGGGGGTTTTTTTTTT‑ACCCGTGAGACGGGCAACA‑Gccg                                                                                                                 >  2:707269/1‑148 (MQ=255)                                   aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGAtt                                                                                                              <  1:113544/151‑1 (MQ=255)                                      tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGc                                                                                                            <  2:309019/150‑1 (MQ=255)                                        cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTTGGCGGC‑A‑GGGGGGCTTTTCTTTTT‑ACCAGTGAGACGGGCCACA‑ACTGATTTccc                                                                                                          >  1:38487/1‑150 (MQ=255)                                           tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGTATTGGGCGCC‑C‑GGGGGGGTTTTTTTTTT‑ACCAGTGAGAAGGGGCACA‑ACTGATTGCCCTTcc                                                                                                      >  2:337111/1‑150 (MQ=255)                                           tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTTGGGGCC‑A‑GGGGGGGTTTTTTTTTC‑ACCAGTGAGACGGGCAAAA‑GCTGATTGCCCtt                                                                                                        >  1:309240/1‑149 (MQ=255)                                                tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCCTTATTGGTGGCC‑C‑GGGGGTGTTTTCTTTTC‑‑CCCGAGGGGCGGGCGACAGGCGGATTGCCCtcccctcc                                                                                                  >  2:430689/1‑144 (MQ=255)                                                  aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTg                                                                                               <  2:813779/151‑1 (MQ=255)                                                    tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGGCTTGGGGGGC‑A‑GGGGGGGTTTTCTTTTC‑ACCACCGAGACCGGCCACA‑ACTGGTTGCCCCTCACccccc                                                                                                >  2:1033182/1‑147 (MQ=255)                                                        cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTGCGTATTGGGGGCC‑A‑GGGGGGGTTTTCTTTTT‑TCCAAGGCGACGGGGCACA‑GGCGAATGGCCTTCACCGCCCGGccc                                                                                           >  1:524452/1‑149 (MQ=255)                                                          cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGGTTTTTTTTTC‑ACCAGTGAGACGGGCAACA‑GCCGATTGCCCTTCACCGCCCGGCCCTg                                                                                         >  1:288148/1‑149 (MQ=255)                                                                aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTGGGCGCC‑C‑GGGGGGGTTTTTTTTTT‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                 >  1:857393/1‑151 (MQ=255)                                                                 gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑ATTGCGTGCTTGGGGGC‑A‑GGGTGGGTTTTCTTTTT‑CCCCGAGAGACCGGGACAA‑GCTGAATGCCCATCACCCCCTGGGCCAgacaggggtg                                                                                >  2:228136/1‑146 (MQ=255)                                                                  tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGCATTTGGCGCC‑A‑GGGTGGGTTTTTTTTTT‑ACCAGTGAGACGGGCAACA‑GCCGATTGGCCCTCCCCCCCTGGcccccggggaagt                                                                                 >  2:757647/1‑139 (MQ=255)                                                                    gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTTGGCGCC‑A‑GGGGGGGTTTTTTTTTT‑ACCCGGGGGACGGGCAACA‑GCCGATTGGCCCTCACCCCCCGGccccgggcgaggt                                                                                 >  2:607448/1‑137 (MQ=255)                                                                    gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcag                                                                             <  1:740225/151‑1 (MQ=255)                                                                      gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTATTGGGCGCC‑A‑AGGGGGTTTTTTTTTTT‑ACCAGTGAGACGGGCAAAA‑GCTGATTGTCCTTCACCCCCCGGCCCCGAGAGAGGTgcag                                                                             >  2:733755/1‑149 (MQ=255)                                                                         aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGATTTGGGGCC‑A‑GGGGGGTTTTTTTTTTT‑ACCCAGGAGACGGGCAACA‑GCCGATTGCCCTTCAACGCCCGGCCCCgagaggggtggagaaa                                                                          >  2:1022726/1‑138 (MQ=255)                                                                         aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                        <  1:370719/151‑1 (MQ=255)                                                                             gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTGTTTGGGGCC‑A‑GGGGGGGTTTTTTTTTC‑ACCACTGAGACGGGGAACA‑AGTGGTTTGCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCg                                                                       >  2:81357/1‑148 (MQ=255)                                                                                  gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑C‑GGGTGGGTTTTTTTTTT‑ACCAGTGGGACGGGCAAAA‑GCTGATTGCCCCTCACCGCCCGGCCCCGGGAGAGTTGCAGCAAGCgg                                                                      >  1:502197/1‑144 (MQ=255)                                                                                    cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTg                                                             <  1:27015/151‑1 (MQ=255)                                                                                     aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGATTGGGGGCC‑A‑GGGGGGGTTTTTTTTTCAACCAG‑GAAAAGGGCAAAA‑GCTGATTGCCCCTCCCCCCCCGGCCCCGAGCGGAGTGCCGCAAGCGGTCCAc                                                                 >  1:676372/1‑146 (MQ=255)                                                                                       cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGGGTATTTTGAGCC‑AGGGGGGGTTTTTTTTTCC‑ACCA‑AGGGACGGGGAACA‑GCCGATTGCCCCTCACCCCCCCTCCCCCTGAGAGGTGCCGCAAGCGGTCCACGCTgggt                                                          >  2:89672/1‑149 (MQ=255)                                                                                        tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTTGGCGCC‑C‑GGGGGGGTTTTCTTTTT‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACCCTGGttt                                                         >  1:577207/1‑151 (MQ=255)                                                                                               aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                  <  1:216050/151‑1 (MQ=255)                                                                                               aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcag                                                  <  2:527491/151‑1 (MQ=255)                                                                                               aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                                   <  2:484499/150‑1 (MQ=255)                                                                                                  aaTCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTATTTGGCGCC‑A‑GGGTGGTTTTTCTTTTT‑CCCCGTGAGAAGGGCAACA‑ACTGATTGGCCTTCCCCGCCCGGCCCCGAGAGGGGTTCAGCAAACGGGCCACCCTGGGTTGCCCccgc                                                 >  1:351957/1‑149 (MQ=255)                                                                                                   aTCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAgg                                              <  2:78803/151‑1 (MQ=255)                                                                                                     cGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                             <  1:618210/150‑1 (MQ=255)                                                                                                        ccAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                         <  2:573468/151‑1 (MQ=255)                                                                                                        caaaCGGGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                         <  2:490319/149‑1 (MQ=255)                                                                                                                cgGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                 <  1:1026542/151‑1 (MQ=255)                                                                                                                  gggAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                               <  1:343309/151‑1 (MQ=255)                                                                                                                    gagaGGCGG‑TTTGGGTATTGGGGGCC‑C‑GGGGGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACC‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGgtt                                >  2:623701/1‑148 (MQ=255)                                                                                                                       agGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           >  1:1066454/1‑150 (MQ=255)                                                                                                                          cGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                        <  2:837735/150‑1 (MQ=255)                                                                                                                            g‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                                      <  2:998169/134‑1 (MQ=255)                                                                                                                            g‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                                      >  1:998222/1‑134 (MQ=255)                                                                                                                              tttGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                     >  1:850342/1‑150 (MQ=255)                                                                                                                                 gCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                  <  1:362960/150‑1 (MQ=255)                                                                                                                                  cGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                <  2:523995/151‑1 (MQ=255)                                                                                                                                   gTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAccgcg               <  1:386710/151‑1 (MQ=255)                                                                                                                                     aTTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg              <  2:772126/150‑1 (MQ=255)                                                                                                                                            gcC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc      <  1:759561/151‑1 (MQ=255)                                                                                                                                                a‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg   >  1:180609/1‑151 (MQ=255)                                                                                                                                                  gggTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa  <  1:284300/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTATTGGGCGCC‑A‑GGGTGGTTTTTCTTTTC‑ACCAGTGAGACGGGCAACA‑GCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366381‑366673

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