Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	49.2%	62.5 / 28.0	67	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (28/6);  new base C (0/33);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT  >  NC_000913/357730‑358021                                                                                                                                                             |                                                                                                                                                aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAtttt‑t                                                                                                                                                 >  1:313906/1‑151 (MQ=255)       aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCCACCGCGCTTCATTTCCGCCAACAA‑CCGCGCA‑TTGATCCAACGCCATGCGC‑AACAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGTTATTTT‑T‑ccgcc                                                                                                                                           >  1:755528/1‑150 (MQ=17)             gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGcaacaa                                                                                                                                    <  2:98430/151‑1 (MQ=255)             gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGcaaca                                                                                                                                     >  2:201156/1‑150 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGcaacaa                                                                                                                                    >  1:517970/1‑150 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAcc                                                                                                                                >  1:303448/1‑151 (MQ=255)                   cgtcccgcATCAGGCAATTCCAGCCGCCGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACaaa                                                                                                                                   <  1:365893/141‑1 (MQ=255)                   cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAccgc                                                                                                                              >  2:737622/1‑151 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAccgcc                                                                                                                             <  1:344105/151‑1 (MQ=255)                     tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAAccgccg                                                                                                                            >  2:765100/1‑151 (MQ=255)                         gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                         >  1:679287/1‑150 (MQ=255)                           aTCAGGCAACTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                                      >  2:822031/1‑151 (MQ=255)                               ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                    >  1:791289/1‑149 (MQ=255)                                  aaTTCCAGCCGCAGCCCTGTGTCAGGGCCTCCCGTGTTTCTTTCCCCCCAACAC‑CCGCGCA‑TTTACCCAACGCCTTGTCA‑AAAAATGCCTT‑C‑GCGGCGGCTTTCTCCCAGTTTTTTT‑C‑CCCCCAAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                               <  2:477983/151‑1 (MQ=255)                                    ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                             >  1:391535/1‑151 (MQ=255)                                    ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                             >  2:282914/1‑151 (MQ=255)                                      cccgCCGCAGACTTGTGCCGGCGCCTCCCGTTATTCTTTCCCGCCAACAC‑CCGCGCT‑TTTATCAACCCCCTTGTGA‑AAAAATGCCTC‑C‑GCGGCGGCTGTCTGCAAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                            <  1:1260429/147‑1 (MQ=255)                                       cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                          >  1:1009703/1‑151 (MQ=255)                                       cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                          >  2:778704/1‑151 (MQ=255)                                          ccGCAGCCTTGTGTCAGCGGCTACCGTGATTCTTTCCCCCCAACAC‑CCGCGCT‑TTTACCCACCGCCATGTCA‑AAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTATTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                        <  2:404631/150‑1 (MQ=21)                                            gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGCAAAAAATGCCTT‑C‑GCGGCGGTTGTCTGCCACCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCCCTCCAGTAGTGGCgg                                                                                                         >  2:1001077/1‑148 (MQ=255)                                             cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  2:632860/1‑148 (MQ=255)                                             cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  2:1283696/1‑148 (MQ=255)                                                 ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggg                                                                                                   >  2:754680/1‑144 (MQ=255)                                                   tgtgTCGGCGCCCACCGGGTTTCTTTCCCGCCACCAC‑CCGCCCT‑TTTCCCCACCGCCATGTCC‑AAAAAGCCCTC‑C‑GGGGCGGTTGTCGCCCAGCTTTTTT‑C‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑ga                                                                                               <  2:1121869/149‑1 (MQ=255)                                                    gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑Gaaa                                                                                             >  2:1297419/1‑150 (MQ=255)                                                         aGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑Gaaaaaaa                                                                                         >  1:738791/1‑145 (MQ=255)                                                           cGGCTCCCTTGATTC‑TTCCCCCCCACAACCCCCGCA‑TTTACCCACCACCTGGTGA‑AAAAATGCGTC‑C‑GGGGGGGCTGTCGGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTa                                                                                     <  1:1095010/151‑1 (MQ=255)                                                             cctccgttGATTC‑TTTCCCCCCACCA‑CCGCCCATTTTACCCACCGCCTTGTAA‑AAAATTGCTTC‑C‑GCGGCGGCTTTCTGCCAGTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAg                                                                                   <  1:53838/144‑1 (MQ=255)                                                              cTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggacaagggaaa                                                                                            >  1:868446/1‑131 (MQ=255)                                                                 ccGTGATTCATTCCCGCCACCAC‑CCGCGCT‑TTTACCCAACGCCATGTAC‑AAAAAGCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCg                                                                               <  2:514545/151‑1 (MQ=255)                                                                 ccGTGATTCATTCCCGCCAACAC‑CCCCGTT‑TTTACCCACCGCTATGTGA‑AAAAAGGCCTG‑C‑GGGGCGGCTGTCTGCCAGCTTTTTT‑C‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCtt                                                                                       <  1:1195444/143‑1 (MQ=25)                                                                 ccGTGATTCATTCCCCCCAACAC‑CCCCGCT‑TTTACCCACCCCCTTGTCA‑AAAATTGCTTC‑C‑GCGGCGCCTGTCTGCCACCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGcc                                                                                 <  2:232740/149‑1 (MQ=255)                                                                  cGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑Gaaaaaaacttaaaccc                                                                                >  2:304735/1‑136 (MQ=25)                                                                     gATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑aaaaaaaattttaaccc                                                                                >  2:454318/1‑133 (MQ=25)                                                                          aTTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAAGTTTTAAgtcgccccgt                                                                          >  1:162546/1‑139 (MQ=25)                                                                            tccccccaccaC‑CCGCGCA‑TTTACCAACCGCCATGTAA‑AAAATGGCCTC‑C‑GGGCCGGTTGTCGCCAAGTTTTTTT‑C‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgt                                                                         <  1:181078/141‑1 (MQ=255)                                                                             ccccccaccaC‑CCGCGCA‑TTTACCCACCGCCATGTAA‑AAAAATCCCTC‑C‑GCGGCGGCTGTCGGCCACCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                   <  2:264930/147‑1 (MQ=255)                                                                              cccccAAAAC‑CCGCGCA‑TTTACCCAACGCAATGTGA‑AAAAAGCCCTG‑C‑GCGGCGGCTGTCTGCCAGTTATTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtg                                                                          <  1:418768/140‑1 (MQ=255)                                                                                    aaaacccGCGCT‑TTTACCCACCCCCTGGTG‑‑AAAAAAGCCTTCC‑GGGGCGGCTTTCTCCCGGTTTTTTT‑T‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                             <  1:768875/146‑1 (MQ=255)                                                                                        cccccgcA‑TTTACCCACCCCCGTGTAA‑AAAAAGGCCTT‑C‑GGGGCGGTTGTCGCCCAGCTATTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                            <  2:313943/144‑1 (MQ=255)                                                                                          ccgcatttttATCAAACGCCTTGCGA‑AAAATTGCCTC‑G‑GCGGCGGCTGTCCGCCAGCTATTTT‑C‑CCCCGCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                        <  2:285873/144‑1 (MQ=255)                                                                                          ccgcaattttACCCAACGCCATGCAA‑AAAATTCCCTG‑C‑GGGGCGGCGGTCT‑CCCGCTTTTTTGC‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                          <  2:64696/142‑1 (MQ=255)                                                                                                 tttACCCACCGCCATGTGA‑AAAAATCCTTC‑C‑GCGGCGGCTGTCTCCCAGCTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                     <  2:371960/147‑1 (MQ=255)                                                                                                  ttATCCAACCCCTTGTGA‑AAAAATCCCTC‑C‑GCGGCGG‑TGTCCGCCAGCTTTTTT‑TCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                 <  1:69644/150‑1 (MQ=255)                                                                                                       caaccccttggtaaaaaaaTCCTTT‑C‑GGGGCGTTTTTCTGCAAGTTTTTTT‑C‑CCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                              <  1:682827/135‑1 (MQ=255)                                                                                                        aaCGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGCAAGA‑GAAAAAACATTCATACCGCCCTGTTTTGTTGCCATCCCTGCTggcgtaaaacag                                           >  1:153577/1‑140 (MQ=25)                                                                                                        aaCGCCATGTGA‑AAAAAGCCCTC‑C‑GCGGCGGCTGTCGGCCAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                          <  2:457509/151‑1 (MQ=255)                                                                                                         aCGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAACATTTAAACCC‑CCATGTTGTGTTGGGAACGCCGCTggggaaaaagagg                                          >  2:896539/1‑139 (MQ=19)                                                                                                            ccATGTGA‑AAAAATCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                       <  2:323203/150‑1 (MQ=255)                                                                                                             cATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAACTTTAATCCCCTCCTGTTGTGTGGGAATCGGTGCTGGGGACAAtaggggaa                                       >  2:283337/1‑142 (MQ=37)                                                                                                              ctggtaaaaaacggCCTC‑C‑GGGGCGGGTTTCTCCCATTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                         <  1:1115258/134‑1 (MQ=255)                                                                                                                    aaaaaaTGCTTT‑C‑GGGCCGCTTCTCTCCAAGTTTTTTT‑C‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                               <  1:422278/150‑1 (MQ=255)                                                                                                                     aaaaaTGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTT‑C‑CCGCCGCAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAggga                                    <  1:791945/145‑2 (MQ=255)                                                                                                                     aaaaaTCCCTT‑CGGCGGGGGCTGTCCGCGAGTTTTTTT‑C‑CCGCCGCAACAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                <  1:321296/150‑1 (MQ=255)                                                                                                                      aaaaGGGCTG‑C‑GGGGCGTTTGTCTGCCAGCTTTTTT‑C‑CCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                             <  1:941114/151‑1 (MQ=255)                                                                                                                           atccttccggggcggcTGTCTCCCACTTTTTTT‑C‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                          <  1:453656/141‑1 (MQ=255)                                                                                                                            ccTC‑C‑GGGGCGGTTGTCTGCCAGTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                        <  2:337091/150‑1 (MQ=255)                                                                                                                                tgggcggggGTTTTCCGCCGGTTATTTT‑C‑CCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                       <  1:1039559/146‑1 (MQ=255)                                                                                                                                 c‑gcgGCGGTTTTGTGCCACTTTTTTT‑C‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                      <  2:726631/148‑1 (MQ=255)                                                                                                                                 c‑gcgGCGGCTGTCTGCCAGCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                           <  1:429640/143‑1 (MQ=255)                                                                                                                                     gggggcTGTCTCCCAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                   <  2:943521/145‑1 (MQ=255)                                                                                                                                     ggcgtcTGTCCGCCAGTTTTTTT‑C‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                       <  1:355790/144‑1 (MQ=255)                                                                                                                                     ggcggcTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GGAATAAAATTTAGCGAGCCCTGTTGTGGGGCAAGAGCTGCTGGGTAAAATGGGGGAAGGGAAAGgcgggggactacaccc                >  2:183611/1‑139 (MQ=255)                                                                                                                                              cTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAACCCCTCCTGTTTTGTGGGAATCGCCGCTGGGTACAAAGTGGGAAGGGGAGGGGGGGGGATCAACTccccaa            >  2:1143956/1‑142 (MQ=255)                                                                                                                                                 ccATCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt     <  2:901376/150‑1 (MQ=255)                                                                                                                                                  cAGCTATTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt     <  2:416592/149‑1 (MQ=255)                                                                                                                                                   aGCTTTTTT‑C‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt  <  2:1032907/151‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                                ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCA‑TTTATCCAACGCCATGTGC‑AAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTT‑T‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT  >  NC_000913/357730‑358021

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