Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	46.4%	53.9 / 30.5	56	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (19/11);  new base C (0/26);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 2.41e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357730‑358023                                                                                                                                                          |                                                                                                                                                   aCTCGTATGTCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCGGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACTCCATGTTCAAAAATGCCTT‑C‑GCGACGGCTGTCTGCCAGCTAttgtt                                                                                                                                                    <  2:793562/151‑1 (MQ=14) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                    >  1:957946/1‑151 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                 <  1:1116491/150‑1 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                 <  2:177564/150‑1 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                                  <  2:808279/149‑1 (MQ=255)     gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgc                                                                                                                                               >  1:1070631/1‑151 (MQ=255)        ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                            <  2:172441/151‑1 (MQ=255)         gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                            >  1:1172195/1‑150 (MQ=255)             gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                         <  1:1232236/149‑1 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaa                                                                                                                                     <  1:1301364/151‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑aaaa                                                                                                                                    <  1:1205811/150‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAc                                                                                                                                   >  1:1303713/1‑151 (MQ=255)                   cGTAGCGCATCAGGCAATTCCACCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAAcc                                                                                                                                  <  2:1219159/149‑1 (MQ=255)                         gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtc                                                                                                                           >  2:554059/1‑150 (MQ=255)                            tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑TCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTtctctc                                                                                                                       >  1:1057746/1‑151 (MQ=255)                              aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                                      >  1:1185853/1‑150 (MQ=255)                                    ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCCCTCCAGGAGTg                                                                                                               >  2:674692/1‑151 (MQ=255)                                      ccACCCGCAGCCTTGTGTCGGCGGCTCCCGTGTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCAACCGCCTGGCGAAAAAATCCCTC‑C‑GGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCCCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                             <  1:1324225/151‑1 (MQ=255)                                       cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                           >  2:715539/1‑151 (MQ=255)                                          ccGCAGCCTTGTCTCGGCGCCTCCTGTTATTTATTCCCGCAAACAC‑CCGCGCATTAATCCAACGCCTTGTG‑AAAAATGTCTTCC‑GCGGCGCCTGTCTGCAATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:171586/143‑1 (MQ=255)                                           cccagccatgtgCCGGCGGCTGCCGTGATCCATCCCCCCCAACAC‑CCGCGCTTTTACCCACCCCGATGTAAAAAAGTTCGTC‑G‑GCGGCGCCTGCCGCCCCGCTATTTTC‑CCCGCCAAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                         <  2:1456042/141‑1 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                        >  1:547115/1‑150 (MQ=255)                                               gcccTGTGTCAGCGCCTACCGTTTTTCTTTCCCGCCAACAC‑CCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTT‑C‑GCGGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAag                                                                                                    <  1:657477/148‑1 (MQ=255)                                                  cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                        >  1:585558/1‑143 (MQ=255)                                                     tttgggCGTCTACCGTGATTCTTCCCCCCCAACAC‑CCGCGCTTTTATCCAGCCCCATGGGAAAAAGTCCCTC‑C‑GCGGCGGCTGTGTCCCAGTTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaaa                                                                                               <  1:664317/144‑1 (MQ=255)                                                     tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggagggggaa                                                                                              >  2:692297/1‑140 (MQ=255)                                                          gcgcCTCCCGTGTTTCATTCCCCCCACCAC‑CCGCGCTTTTCCCCACCCCTTTGTAAAAAAAT‑CCTT‑CGGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg                                                                                            <  2:98300/147‑1 (MQ=255)                                                            ggCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGaaaaaaattt                                                                                        >  1:858129/1‑142 (MQ=255)                                                            ggCTACCGTAATTCTTTCCCGCCAACAC‑CCGGGCTTTTACCCACCCCCCTGTAAAAAAACGCCTG‑C‑GCGGGGGCTGCCTCCCAGTTTTTTTC‑CCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAt                                                                                              <  1:15446/143‑1 (MQ=255)                                                            gcctcccGTGTTTCTTTCCCGCCACCAC‑CCCCGCTTTTACCAACCGCCAGGT‑CAAAAAAGCCTTCC‑GCGGCGGTTGCCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGaa                                                                                                <  2:230456/136‑1 (MQ=255)                                                            gcctcccGTGATTAATTCCCGCAACCAC‑CCGCGCTTTTACCCACCCCCATGTGAAAAAAGCCCTT‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCt                                                                                          <  1:949156/142‑1 (MQ=255)                                                             gCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGgaaaaaattt                                                                                        >  1:181921/1‑141 (MQ=255)                                                              cccccGTGTTTCTTTCCCGCCAACAC‑CCCCGATTTTACCCACCGCTATGTCAAAAAATCCTTC‑C‑GCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCACAAA‑AAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAg                                                                                     <  2:1217260/147‑1 (MQ=255)                                                              cTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGAAGaaaaaaactt                                                                                        >  2:1323704/1‑140 (MQ=25)                                                                ccggttaTTAATTCCCGCCAACAC‑CCCCCCATTTACCCACCCCATTGTGAAAAAATCCCTC‑G‑GCGGCGGCTTTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc                                                                                  <  2:1346563/145‑1 (MQ=255)                                                                aCCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGAAAAAACTTTAAAccc                                                                                  >  1:598663/1‑151 (MQ=255)                                                                  cggtttttaTTTCCCCCAACAACCCGCGCATTTACCCACCCCCTTGTAAAAAATTCCCTC‑C‑GGGGCGGTTGTCGGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCtt                                                                                         <  2:1254825/135‑1 (MQ=255)                                                                           ttccccccAAAAC‑CCGCCCTTTTACCCACCCCAATGTAAAAAAATCCCTG‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCCCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                           <  1:592878/141‑1 (MQ=255)                                                                           tccccGCCACCAC‑CCGCTTTTTTACCCCACGCCATGTAAAAAAATCCTTC‑C‑GCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCACACA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGACCGGCAtgttg                                                                         <  2:608206/146‑1 (MQ=255)                                                                              cccccAACAC‑CCGCGCTTTTACCCACCGCCTTGTGAAAAATTGCTTT‑G‑GGGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                      <  1:252997/146‑1 (MQ=255)                                                                               ccccaaccacccccgcTTTTACCCACCCCCATGGGAAAAATTCCCTG‑C‑GGGGCGGTTGTC‑GCCAATTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                           <  2:954952/131‑1 (MQ=255)                                                                                   aacac‑cCGCGCATTTACCCACCCCCTTGTGAAAAAATTCCTT‑C‑GGGGCGGCTTTCTGCCATTTTTTTTC‑CCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                <  1:1368108/150‑1 (MQ=255)                                                                                         ccgcacttttccccACCGCCATGTAAAAAAGTGCATG‑G‑GCGGGGGCGGTCGGCGGTTTTTTTGC‑CCGCCACAAA‑AAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                               <  2:688549/134‑1 (MQ=255)                                                                                           gcgcATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAAAAAAAAACTTTAA‑CCCCTCCTGTTTTGTTGGAATccc                                                             >  2:1309730/1‑144 (MQ=25)                                                                                                   aTCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGAAAAAAATTTTAACCCCTCCATTTGTGTTGCGAACGCTGCTggggaaaa                                                 >  2:774083/1‑144 (MQ=25)                                                                                                     ccaccgccagggaaaaaaaTGCCTG‑C‑GCGGCGCCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  2:1276166/132‑1 (MQ=255)                                                                                                         ccccttgtaaaaaaaTGCTTC‑C‑GCGGCGGTTGCCGGCCAGTTTTTTTC‑CCCCCGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                               <  2:852803/136‑1 (MQ=255)                                                                                                           ccagggaaaaaaaTGCCTC‑C‑GGGGCGGTTTTCGGCCAGTTTTTTTC‑CCCCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                           <  2:751548/140‑1 (MQ=255)                                                                                                               gtcaaaaaaTCCCTG‑C‑GGGGCGGTTGTCTCCCATTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                       <  1:62683/144‑1 (MQ=255)                                                                                                               gtaaaaaaaGGCCTC‑C‑GCGGGGCTTTTCTGCCAGTTTTTTTC‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                       <  2:75230/144‑1 (MQ=255)                                                                                                                taaaaaaaTCCTTC‑C‑GCGGCGCTTGTCTGCCAGTTTTTTTC‑CCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                      <  2:476208/145‑1 (MQ=255)                                                                                                                taaaaaaaTCCCTC‑C‑GCGGCGGCTTTCTCCCATTTTTTTTT‑CCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATAGCTGCTGGTAACAATGTGGAAAGGGa                                     <  1:97284/146‑1 (MQ=255)                                                                                                                aacAAAAATGCCTC‑C‑GCGGCGGTTGTCTGCCGGTTATTTTC‑CCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                   <  2:1143154/149‑1 (MQ=255)                                                                                                                     aaaTGCCTG‑C‑GGGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                              <  2:588244/151‑1 (MQ=255)                                                                                                                      aaTGCCTT‑CGGGGGCGCCTTTCTCCCGTCT‑TTTTT‑CCCCCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                              <  2:1513976/151‑1 (MQ=255)                                                                                                                           cTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATTGCTGCTGGGAACAAAGGGGGAAGGGGAgggggggg                            >  2:1343207/1‑140 (MQ=255)                                                                                                                                   ggcggcTGTCTGCCAGTTTTTTTC‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                   <  1:437978/150‑1 (MQ=255)                                                                                                                                   ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGACGAGGAAAAGAATATTACCCGACATATTGTGTGGGGCTCGGCGCTGGGGACAATATTGAAAAGGGAGGGCGGGGAGTCAGc                     >  2:166503/1‑148 (MQ=255)                                                                                                                                                 atcttttttc‑CCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa     <  1:190109/145‑1 (MQ=255)                                                                                                                                                    ttttttc‑CCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  2:165898/148‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                   ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357730‑358023

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