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%	61.6 / 29.6	65	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (21/12);  new base C (0/32);  total (21/44)
Fisher's exact test for biased strand distribution p-value = 7.97e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026                                                                                                                                                          |                                                                                                                                                       aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                        >  1:695833/1‑151 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                     >  2:109064/1‑151 (MQ=255)     gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                     <  2:38136/150‑1 (MQ=255)       aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                  <  1:459279/151‑1 (MQ=255)        ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                  <  2:261586/150‑1 (MQ=255)        ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCcgccg‑c                                                                                                                                                >  1:118607/1‑151 (MQ=255)        ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCcgccg‑c                                                                                                                                                >  1:176179/1‑151 (MQ=255)          ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑caa                                                                                                                                              <  1:942631/151‑1 (MQ=255)          ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGGCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCGT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑ca                                                                                                                                               <  2:448622/150‑1 (MQ=255)           cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑caac                                                                                                                                             >  2:894672/1‑151 (MQ=255)           cTGATAAGCGTAGCGCATCAGGCAATCCCAGCCGCAGACCTATGTCAGCGCCTACCGCGATTCATTTCCGCCACCAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑caa                                                                                                                                              <  2:1198898/150‑1 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACaaa                                                                                                                                          <  2:1080937/151‑1 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTTCCAGCTATTTTTCCGCCG‑CAACAAAAcc                                                                                                                                       <  2:1093957/151‑1 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAcc                                                                                                                                       <  2:725644/151‑1 (MQ=255)                   cGTAGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccgc                                                                                                                                     <  2:497664/151‑1 (MQ=255)                       gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGtt                                                                                                                                 >  2:831014/1‑151 (MQ=255)                          cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctc                                                                                                                              >  2:1059713/1‑151 (MQ=255)                           aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctct                                                                                                                             >  2:1248067/1‑151 (MQ=255)                            tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctctc                                                                                                                            >  1:726310/1‑151 (MQ=255)                             cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGCGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCCATTTTTCCGCCG‑CAACCAAACCGCCGTTCTCTcc                                                                                                                           >  2:747198/1‑151 (MQ=255)                              aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTcc                                                                                                                           >  1:579912/1‑150 (MQ=255)                                  aaTCCCAGCCGCAGCCTTGTGTCAGCGCCTCCCTTGTTTCATTCCCGCCAACACCCCCGCTTTTATCAAACGCCA‑TGGGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCagta                                                                                                                       <  2:694596/150‑1 (MQ=255)                                      ccAGCCGCAGCCCTGTGTCGGCGGCTCCCGTTTTTCTTTCCCGCAAACACCCGCGCTTTTACCCACCGCCA‑TGTAAAAAAATCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCCCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                    <  1:1303261/149‑1 (MQ=255)                                      ccACCCCCAGCCTTGTGTCAGCGCCTCCCTTGATCCATTCCCCCCAACCCCCGCCCTTTTACCCACCCCCA‑TGTAAAAAAATGCTTC‑C‑GCGGCGCTTTTCTGCAAGTTTTTTTCCCCCCG‑AAACAAAACCGCCGTTCTCTCCagta                                                                                                                       <  2:534347/146‑1 (MQ=255)                                       cccccGCAGCCTTGTTTCCGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGTTTTTACCCACCCCCA‑TGCACAAAAATCCTGC‑C‑GCGGCGGTTGTCTGCAAGTTTTTTTTCCGCCG‑CAAAAAAACCGCCGTTCTCTCCagta                                                                                                                       <  2:479296/142‑1 (MQ=255)                                       cAGCTGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑GGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                 >  2:782552/1‑151 (MQ=255)                                       cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                 >  1:892113/1‑151 (MQ=255)                                        aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                >  1:283682/1‑151 (MQ=255)                                          ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCA‑TGTGCAAAAATGCCTC‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              <  1:433027/151‑1 (MQ=255)                                           cGCAGCCTTGTGTCGGCGGCCCCCGTAATTCATTCCCCCCAACACCCGCGCATTTACCCAACGCCA‑TGTCAAAAAGTCCCTT‑C‑GCGCGGGCTGTCTCCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                    <  2:772590/144‑1 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                             >  2:691697/1‑150 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                               >  1:415606/1‑149 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                >  2:1226613/1‑148 (MQ=255)                                             cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                >  2:1063172/1‑146 (MQ=255)                                               ggccTGTGTCGGCGCTTGCCGTGATTCTTTCCGGCAAACACCCGCGCATTTCTCCAACGCAA‑TGT‑CAAAAAAGCCTTCC‑GGGGCGGTTGTTGGCCGGTTTTTTTCCCGCCC‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                         <  1:966288/149‑1 (MQ=255)                                                 cttgtgTCAGCGGCTCCCGTGATTCTTTCCCCCAAACACCCGCGCTTTTTCCCAACTCCA‑TTCAAAAAATTGCTTC‑C‑GCGGCGCCTCCCTCCCCGTTTTTTCCCCCCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                       <  1:895282/148‑1 (MQ=255)                                                   tgtgTCAGCGCCTCCCGTGTTTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCA‑TGTGAAAAAATCCCTT‑C‑GCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                       <  1:916213/148‑1 (MQ=255)                                                           cgcctcccgtttttCATGCCCCCCAACACCCGCCCTTTTCCCCACCGCCA‑TGAGAAAAAATGCCTC‑C‑GCGGGGGCTGTCTGCCGGTTATTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATa                                                                                                 <  2:356490/133‑1 (MQ=255)                                                            ggCTCCCGTGTTTCTTTCCCCCCAAAACCCGCGCTTTTATCCAACGCCT‑TGTGAAAAAATGCTTT‑C‑GGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCttt                                                                                            <  1:1197288/149‑1 (MQ=255)                                                                       ttCATTCCCCCCAACACCCCCGCATTTCCCCACCGCCA‑TGTGCAAAAATGCTTT‑C‑GCGCCGCCTTTCTCCCAGCTTTTTTCCCCCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCa                                                                                  <  1:417763/148‑1 (MQ=255)                                                                         cttgcccGCCAACACCCGCGCATTCCCCCACCGCCG‑GGGGCAAAAAGGCCGC‑G‑GGGGGGGCGGGCGGCGGGTTATTTTCCCGCCG‑AAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCg                                                                                     <  2:78909/138‑1 (MQ=255)                                                                         cttccccGCCAACACCCGCGCTTTAATCCACCGCCA‑TGTG‑AAAAATGCCCT‑CGGCGGCGCTTGTCTCCCAGCTATTTTCCCGCCA‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCa                                                                                  <  2:1103724/141‑1 (MQ=255)                                                                           ttccccccaccaCCCGCGCTTTTA‑CCAACCCCATTGTAAAAAAATGCCTT‑C‑GCGCGGCCTGTCTGCCGTTTATTTGCCCGCCA‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgt                                                                               <  2:824138/141‑1 (MQ=255)                                                                           ttccccccAACACCCCCGCTTTTACCCACCCCAT‑TGTGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttgt                                                                            <  1:109085/144‑1 (MQ=255)                                                                            ttCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AAAAAACTTTTAGCCCCTCCtgtt                                                                              >  1:37287/1‑147 (MQ=25)                                                                             tCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGAGG‑AAAAAACTTTaaaaccctaccttttgt                                                                           >  2:601505/1‑135 (MQ=21)                                                                             ccccccAACACCCGCTCTTTCACCAAACGCCA‑TGTCCAAAAATCCCTC‑C‑GGGGCGCTTTCCGGCCAGTTTTTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                             <  2:39567/143‑1 (MQ=255)                                                                                 ccaccaccccctcATTTTCCCACCCCCG‑TGTCCAAAAAGCCCTG‑C‑GCGGCGGCTGTCGGCCAGCTATTTTCCCGCCA‑AAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                     <  2:164420/139‑1 (MQ=255)                                                                                 ccAACACCCGCGTATTTACCCACCGCCA‑TGTGAAAAATTGCCTT‑C‑GGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                       <  1:862678/149‑1 (MQ=255)                                                                                   aacaccCGCGCTTTTATCCACCGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGTTATTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                   <  1:54560/151‑1 (MQ=38)                                                                                    acaacCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAACTTTAA‑CCCCTCATGTTTTGTGGGaa                                                                     >  1:1272931/1‑148 (MQ=25)                                                                                     cacccccgattttccccAACGCCA‑GGTGAAAAAATGCTTT‑C‑GGGGGGGTGGTCTGCCAGTTATTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                                  <  2:891505/135‑1 (MQ=255)                                                                                              aTTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACACAACCGCCGTTCTCTCCAGTAGTGGCGGCGCACGC‑G‑AAATATCATTTAAACCGTCTGGTGGTGTGGGCGTCACTGCTgg                                                           >  2:770576/1‑148 (MQ=21)                                                                                               ttcccccAACGCCA‑TGTGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCGGCTTTTTTCCCGCCA‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                        <  2:2457/145‑1 (MQ=255)                                                                                                       acccccttgacaaaaaaTGCCTC‑G‑GCGGCGCCTGTCGGCCAGTTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTgg                                               <  2:443190/139‑1 (MQ=255)                                                                                                          ccA‑TGTAAAAAAACCCCTT‑C‑GCGGCGGTTGTCTGCCATCTTTTTCCCCGCCG‑AAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑aa                                                    <  2:1022275/142‑1 (MQ=255)                                                                                                           ccttgtgAAAAAATGCCTT‑G‑GCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGaaa                                            <  1:1144117/147‑1 (MQ=255)                                                                                                                acaaaaagTGCTTT‑C‑GCGGTGTCTGCTTCCCAGTTTTTTTCCCGCCA‑AAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGaaa                                            <  2:550665/142‑1 (MQ=255)                                                                                                                   aaaaaTGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTg                                    <  1:1215201/150‑1 (MQ=255)                                                                                                                          ccTG‑C‑GGGGCGGCTGTCTCCCAGTTTTTTCCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAAt                            <  1:225901/151‑1 (MQ=255)                                                                                                                             tcc‑gcgGGGCCTGCCGGCCAGCTTTTTTCCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGaa                             <  2:176129/146‑1 (MQ=255)                                                                                                                                   ggcggtTTTCTCCCGGTTTTTTTCCCGCCA‑CAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:156785/144‑1 (MQ=255)                                                                                                                                   ggcggtTGTCTGCCAGCTTTTTTCCCCCCC‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGaa                             <  1:73981/143‑1 (MQ=255)                                                                                                                                               ccactttttttcCCGCCC‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa         <  1:266082/144‑1 (MQ=255)                                                                                                                                               ccATCTTTTTTCCCGCCG‑AAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  1:1165358/145‑1 (MQ=255)                                                                                                                                                ccgtcttttttCCCCCGAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                <  2:532071/138‑1 (MQ=255)                                                                                                                                                  gttATTTTCCCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa        <  1:579935/147‑1 (MQ=255)                                                                                                                                                  gCTTTTTCTCCCCGG‑CAACAAAAGCGGCCTTCTCTCCAGGAGTGGCGGGGCCAGA‑G‑GAATGGGTTGAAGGCCGTCATGTTGGGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:889798/144‑1 (MQ=255)                                                                                                                                                      ttttcCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg     <  1:987768/148‑1 (MQ=255)                                                                                                                                                      ttttcCCGCCG‑AAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg     <  2:360540/148‑1 (MQ=255)                                                                                                                                                      ttttcCCGCCA‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa    <  1:1171027/149‑1 (MQ=255)                                                                                                                                                       ttttCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAACTTTTAACCCCCCCTGTTGTGTGGGGATCGCTGCTGGGGACAAATGGGGAAGGGGAGGGGGGGGGAA‑CAACTCCCGAAACCCGCCAAtt      >  1:300546/1‑146 (MQ=255)                                                                                                                                                       tttcCCCCCG‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  1:701640/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                       ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCA‑TGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAC‑AATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026

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