Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	606,612	0	G	A	44.6%	106.7 / 32.2	74	Q258* (CAG→TAG)	ybdK	weak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (35/6);  new base A (0/33);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 3.09e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGA  >  NC_000913/606464‑606760                                                                                                                                                           |                                                                                                                                                      aGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGaa                                                                                                                                                     <  2:724328/151‑1 (MQ=255)  gCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGaa                                                                                                                                                     <  1:139923/150‑1 (MQ=255)   cGGTCGCCGATCTCCAGTGTGCGGATCGGTGATGACGCATTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCA‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAAt                                                                                                                                                    >  1:444385/1‑150 (MQ=255)          cgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCg                                                                                                                                           >  2:246949/1‑151 (MQ=255)          cgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCg                                                                                                                                           >  1:851317/1‑151 (MQ=255)          cgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCg                                                                                                                                           >  1:833527/1‑151 (MQ=255)           gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGc                                                                                                                                          >  2:980312/1‑151 (MQ=255)           gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGc                                                                                                                                          >  2:827450/1‑151 (MQ=255)           gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGc                                                                                                                                          >  1:144704/1‑151 (MQ=255)            aTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGcc                                                                                                                                         >  1:524983/1‑151 (MQ=255)                ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCa                                                                                                                                        <  1:718793/148‑1 (MQ=255)                            tgggTGATGCCCCCCCCAACCCAATAGCGACAGGCCTGAAAGCGTTTAAT‑TTTATACAGCAGTTATTTTTTTTCT‑TG‑TT‑TTTAAACGGGGCTTCCC‑T‑TCAGAAACCGGGGGGGGTTACCTTAAATA‑AACCCCGCCATATTTACTGCGTgg                                                                                                                          <  1:938322/148‑1 (MQ=255)                              ggTGATGCCGCCTTCAACCCCAAAGCACCAGG‑CAGGAAAACGTTGAA‑ATTTAACAGGAGGTTATTTTTTTTC‑TTGTT‑TTTAAAG‑CGGCGGTCC‑GTTAATTAACCAGGGGGCGGTACCGTAAATT‑AACCCCGCCATATTTACTGCGTGGCTa                                                                                                                       <  2:1045464/151‑1 (MQ=255)                                   tgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGt                                                                                                                  >  2:34670/1‑151 (MQ=255)                                    gaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGtt                                                                                                                 >  2:340328/1‑151 (MQ=255)                                      cGCCTTCAACCCCATAGCGACAGGCCTGAAAACGTTTAAT‑TTTATACAGCAGTTATTTTTTTTCC‑TGATT‑TTTAACG‑GGGCGTCCC‑G‑TCATTACCCGGTGGGGGGTAGCCTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGt                                                                                                                  <  2:1159738/148‑1 (MQ=255)                                             agccccATAGCGCCAGGCCTGAAAAGGGTTGAT‑TTTATAAAGAAGTTATTTTTTTTCT‑TG‑TT‑TTTAAACGGGGCTTTCC‑T‑TAAGTACCGGGGGGGGGGTAGCCGGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTa                                                                                                        <  1:290047/148‑1 (MQ=255)                                               gCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATc                                                                                                      >  1:1144944/1‑151 (MQ=255)                                                ccaatggCGACGGGCCTGAAAACGTTTAAA‑TTTATACAGCGGTTATTTTTTTCCC‑TGTTT‑TTTAACC‑GGGCGTCCC‑G‑TAAGTACCCGGGGGGGGGTAGCCTAATTT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTAt                                                                                                       <  2:398219/143‑1 (MQ=255)                                                   aaaGCGACGGGCGAGAAAACGGTTAAT‑TTTATACAGAAGGTAATTTTTTTCT‑TGTTG‑TTAAACC‑GGGCTTTCC‑G‑TAAGTAACAGGTGGGGGGTAG‑CTAAATTAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCAt                                                                                                   <  2:873657/148‑1 (MQ=255)                                                     aGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAc                                                                                                >  2:709679/1‑151 (MQ=255)                                                        gACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCg                                                                                             >  2:1148152/1‑151 (MQ=255)                                                         acgggCCGGGAAAGGTTTGAT‑TTTATACAGCAGTTATTCTTTTCCC‑TGTTT‑TTAAAAG‑GGGCGTCCC‑G‑TCAGTACCCGGTGGGCGGTAGCTTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCa                                                                                                 <  2:642621/143‑1 (MQ=255)                                                             ggccgaaaaGGGTTTAATTTTATCCAGCAGTTAATTTTTTTCT‑TGTTG‑TTTAAGG‑GGGCGTCCC‑G‑TCAGTAACCGGGGGGGGGTACCCTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGaa                                                                                           <  1:505996/143‑1 (MQ=255)                                                              ccTGGAAGCGTTTAAA‑TTTATACAGCAGTTAATTTTTTTCC‑TGATT‑TTAAAGG‑GGGCTTTCC‑G‑TAAATAACAGGGGGGGGGTAGCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTc                                                                                       <  2:1104085/151‑1 (MQ=255)                                                              ccTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTc                                                                                       >  2:145132/1‑151 (MQ=255)                                                                 gaaaaCGTTTGAA‑TTTATACAGAAGTTATTTTTTTCCC‑TTATT‑TTAAACC‑GGGCGTCCC‑G‑TCAGTAACCGGTGGGCGGTACCATAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCAc                                                                                    <  2:412508/149‑1 (MQ=255)                                                                  gAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCAcc                                                                                   >  1:698640/1‑151 (MQ=255)                                                                    agcGGTTAAT‑TTTATACAGCAGTTATTTTTTTTCT‑TGATT‑TTAAACC‑GGGCTTCCC‑G‑TCAGTAACCGGGGGGCGGTGGCTTAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGt                                                                                 <  2:554900/149‑1 (MQ=255)                                                                       cgtttaattttAAACGGCAGTTAAT‑TTTTTCC‑TTGTGTTTTAACC‑GGGTTTCCC‑G‑GCATTACCGGGGGGGCGGTCCCCTGATTT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGt                                                                                 <  1:1150242/141‑1 (MQ=255)                                                                        gTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGcc                                                                              >  1:985676/1‑150 (MQ=255)                                                                             attttATACGGCAGTAATTTTTTTTCT‑TAATG‑TTAACGG‑GGGCTTCCC‑G‑GCAGAAACGGGGGGGGGGAACCCAAAATA‑ACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaa                                                                            <  1:254237/146‑1 (MQ=255)                                                                             attttATACAGCAGTTATTCTTTTTCC‑TGTTG‑TTAAACG‑GGGCGTTCC‑G‑TCAGTAACCGGGGGGCGGTACCTTAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaaa                                                                           <  1:432141/147‑1 (MQ=255)                                                                               tttATACAGCAGTTATTTTTTTTCC‑TGTTT‑TTTAAAG‑GGGCGTCCC‑G‑TCAGAACCCGGTGGGGGGTAGCTTAAATA‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATg                                                                        <  2:1151547/150‑1 (MQ=255)                                                                                ttATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATg                                                                        >  2:1248926/1‑149 (MQ=255)                                                                                   taCAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGAc                                                                   >  1:648664/1‑151 (MQ=255)                                                                                     cagaagttattttttttCC‑TTTTT‑TTTAAGG‑GGGCTTTCC‑G‑TAAGTAACCGGGGGGCGGTAGCCTAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGAc                                                                   <  2:972832/137‑1 (MQ=255)                                                                                      agcagGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgg                                                                >  1:39026/1‑151 (MQ=255)                                                                                          ggTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                            >  2:552762/1‑151 (MQ=255)                                                                                          ggTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                            >  1:433797/1‑151 (MQ=255)                                                                                            tAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGa                                                           >  1:954983/1‑150 (MQ=255)                                                                                               ttttttcCC‑TGTTT‑TTAAAAG‑GGGCTTTCC‑G‑TCAGTACCCGGTGGGGGTTAGCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAat                                                         <  2:363861/147‑1 (MQ=255)                                                                                                 tttttttctttttttttAAGC‑GGGCTTCCC‑G‑TAATTACCGGGGGGGGGGTGGCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgggg                                                              <  2:33545/130‑1 (MQ=255)                                                                                                 tttttCT‑TGTTT‑TTTAAGC‑GGGCTTTCC‑T‑TCAGTACCAGGTGGGGGGGGCCCAAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAat                                                         <  1:567415/147‑1 (MQ=255)                                                                                                 tttttCC‑TGAT‑‑TTTAAACGGGGCTTCCC‑G‑TCAGTACCGGGGGGGCGGTAGCCTAAATT‑ATTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAatat                                                       <  1:1027436/149‑1 (MQ=255)                                                                                                   tttCCTTGATG‑TTTAACC‑GGCCGTCCC‑G‑TCAGAACCCGGGGGGGGGTGGCATAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATcc                                                     <  2:777683/150‑1 (MQ=255)                                                                                                    ttCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc                                                    >  2:811634/1‑149 (MQ=255)                                                                                                      cc‑TGATT‑TTTAAGC‑GGGCGTTCA‑G‑TGAGAAACGAGGGGGGGGAAGTCTAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGt                                                 <  1:733443/150‑1 (MQ=255)                                                                                                       cttgttG‑TTTAACG‑GGGCTTTCC‑G‑TCAGAAACCGGGGGGGGGTAGCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCa                                                   <  2:147067/143‑1 (MQ=255)                                                                                                         tGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                               >  1:1134296/1‑150 (MQ=255)                                                                                                            gattttAAAC‑GGGGCTTCCGG‑TCAGTACCGGGGGGGCGGTAGCGTAAAGT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTg                                                <  1:657574/145‑1 (MQ=255)                                                                                                               tttAACC‑GGGCGTTCC‑G‑TCAGTAACCGGTGGGCGGTAGCCTAAATA‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATc                                          <  2:148502/150‑1 (MQ=255)                                                                                                               ttaaaaG‑GGCCGTTCC‑G‑ACATTACCCAGGGGGGGGTAGCCTAAATT‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATc                                          <  2:66358/147‑1 (MQ=255)                                                                                                               ttaaaaC‑GGCCTTTCC‑G‑TCAATAACCGGGGGGCGGTAGCGTAATTT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGt                                                 <  2:1093703/140‑1 (MQ=255)                                                                                                                ttAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCtt                                        >  2:1068944/1‑151 (MQ=255)                                                                                                                 tAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCttt                                       >  2:372014/1‑151 (MQ=255)                                                                                                                 tAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCttt                                       >  1:637329/1‑151 (MQ=255)                                                                                                                  aagc‑GGGCTTCCC‑G‑TAAGAACCCGGGGGGGGGTGGCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCtt                                        <  2:546043/146‑1 (MQ=255)                                                                                                                  aagc‑GGGCGTCCC‑G‑TAATTACCGAGGGGGGGGTAGCTTAAATA‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTa                                      <  2:188751/148‑1 (MQ=255)                                                                                                                  aaaC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCttt                                       >  1:475466/1‑150 (MQ=255)                                                                                                                   aaC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTa                                      >  2:417423/1‑150 (MQ=255)                                                                                                                       ggccttccc‑G‑CCAGTACCTAGGGGGCGTTACCTTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATg                                   <  2:650639/143‑1 (MQ=255)                                                                                                                           gTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTc                              >  2:853675/1‑151 (MQ=255)                                                                                                                            tccc‑G‑ACAGTAACCGGGGGGCGGTACCCTAAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTc                              <  2:163145/148‑1 (MQ=255)                                                                                                                                   tCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatcat                        >  2:335053/1‑151 (MQ=255)                                                                                                                                   tCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatcat                        >  1:364537/1‑151 (MQ=255)                                                                                                                                        acccggggggCGGTAGCTTAAATT‑ACTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatcatc                       <  2:321603/140‑1 (MQ=255)                                                                                                                                          cgggTGGGGGGTAGCATAAATT‑ACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGCGCAGATCTTTAATGCTGTCGATCATCgtgg                   <  1:1079748/146‑1 (MQ=255)                                                                                                                                            aGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGgtgt                >  1:999712/1‑150 (MQ=255)                                                                                                                                             gTGGGGGGTAGCCTAAATT‑ACTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGgtgt                <  2:462917/149‑1 (MQ=255)                                                                                                                                             gTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTa               >  2:127295/1‑150 (MQ=255)                                                                                                                                              ggggCGTTACCATAAATT‑ATACCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatca                         <  2:473668/138‑1 (MQ=255)                                                                                                                                                       gCCTAAATA‑AACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACa      <  1:301528/149‑1 (MQ=255)                                                                                                                                                          tGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGa  >  2:912684/1‑150 (MQ=255)                                                                                                                                                           |                                                                                                                                                      AGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TGATG‑TTTAAAC‑GGGCGTTCC‑G‑TCAGTAACCAGTGGGCGGTAGCCTGAATT‑AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGA  >  NC_000913/606464‑606760

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