Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	45.9%	56.6 / 39.2	63	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (26/7);  new base A (0/28);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 3.60e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357740‑358029                                                                                                                                                                |                                                                                                                                            cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAAAAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGc‑aac                                                                                                                                             >  1:636452/1‑151 (MQ=255)    aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACaaa                                                                                                                                          <  1:278950/151‑1 (MQ=255)     tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACaaa                                                                                                                                          >  1:319272/1‑150 (MQ=255)      aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAAc                                                                                                                                        >  2:922606/1‑151 (MQ=255)           tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAAccgccg                                                                                                                                   >  1:238777/1‑151 (MQ=255)            aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAAccgccg                                                                                                                                   >  2:1050806/1‑150 (MQ=255)            aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGt                                                                                                                                  >  2:305732/1‑151 (MQ=255)              cgcacaaGGCAATCCCAGCCGCAACCCTTTTTCACCGCCTACCGTGTTTCATTCCCG‑CCAACACCCGCCCATTT‑AC‑CCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTCCCG‑CCGC‑AACAAAACCGCCGTtc                                                                                                                                <  1:1002521/145‑1 (MQ=18)              cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAAccgcc                                                                                                                                    <  2:903654/147‑1 (MQ=255)              cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTtc                                                                                                                                >  2:615444/1‑151 (MQ=255)               gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTtct                                                                                                                               >  2:307164/1‑151 (MQ=255)                cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTtct                                                                                                                               >  1:1155451/1‑150 (MQ=255)                   cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTcc                                                                                                                           >  2:604462/1‑151 (MQ=255)                    aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AATAAAACCGCCGTTCTCTcc                                                                                                                           >  1:613551/1‑150 (MQ=255)                        aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCagt                                                                                                                        >  1:609362/1‑149 (MQ=255)                           tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                   >  1:1074371/1‑151 (MQ=255)                           tCCAGCCGCAGAACTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                   >  1:1195034/1‑151 (MQ=255)                            ccACCCGCAGCCCTGTGTCGGCGGCTCCCGTAATTCATTCCCG‑CCAACACCCGCGCTTTT‑AC‑CCACCGCCATGTGAAAAAATGCTTC‑CGCGGCGG‑C‑TGTCGGCC‑AGTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                     <  2:118233/148‑1 (MQ=255)                             cagccgcagccgTGTGTCAGCGCCTCCCTTGTTTCATTCCCG‑CAAACACCCGCCCTTTT‑AC‑CCACCGCCTTGTGAAAAAATCCCTC‑CGCGGGGG‑C‑TGTCGGCC‑AGTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCagta                                                                                                                       <  2:689218/145‑1 (MQ=255)                                ccGCAGCCCGGTGTCGGCGGCTACCTTGATTCTTTCCCC‑CCACCACCCGCGCATTT‑AT‑CCACCGCCATGTGAAAAAATGCCTC‑CGCGGCGG‑C‑TGCCTGCC‑AGCTATTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                               <  2:238777/150‑1 (MQ=255)                                ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  1:1179454/1‑151 (MQ=255)                                   gcttccTGTGGCAGGGGCTCCGGTATTTCATTCCCC‑CCAACAACCG‑GCATTTCAG‑CCCCCCCAATG‑GCAAAAAAGCCCC‑CCGGGGAG‑C‑CTTCTTCCAAGTTATTTTCCCC‑CCGC‑AAAAACACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              <  1:645399/144‑1 (MQ=255)                                   cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  2:913766/1‑148 (MQ=255)                                   cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  2:30420/1‑148 (MQ=255)                                    aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                            >  2:1001280/1‑149 (MQ=255)                                       ccTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  2:376343/1‑144 (MQ=255)                                           tgccgcctgCTCCCGTTATTCTTTCCCC‑CCACCACCCCCTCTTTT‑AC‑CCACCCCCTTGTAAAAAAATCCTTC‑CGGGGGGG‑C‑TGTCCGCC‑AGTTATTTCCCCG‑CCAA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                         <  1:1056301/137‑1 (MQ=255)                                            gtCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggggaaaa                                                                                                    >  2:1189703/1‑139 (MQ=255)                                              tcggggtccccggtaTTCCTTTCCCCCCAAAACCCGCACTTTT‑CC‑CCAACCCCTGGTGAAAAAATCCTTC‑CGGGGGGG‑C‑TTCCTGCC‑AGTTATTTCCCCG‑CCGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                  <  1:1172012/136‑1 (MQ=255)                                              cggCGGCTCCCGTTTTTCTTTCCCC‑CAAACACCCGCGCATTT‑AC‑CCACCGCCTGGTGAAAAAATCCCTC‑CGGGGCGG‑C‑TGTCTGCC‑AGTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                  <  2:842211/147‑1 (MQ=255)                                               aGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGaaaaaaat                                                                                                >  1:645656/1‑145 (MQ=255)                                                 cGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaa                                                                                                   >  2:804429/1‑143 (MQ=255)                                                        cagtgatgaTTTCCG‑GCAACCAACCCGCGATT‑CTCCCAAGGCTATTTGCAAAAATG‑CTT‑CGGCGCGGAC‑TCTCTCCC‑AG‑TTTTTTTCCGACCTC‑ATAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                                <  2:725452/134‑1 (MQ=255)                                                         gtatttCATTCCCG‑CCACCACCCCCGCTTTT‑AC‑CAACCCCCATGT‑CAAAAATTCCCT‑CCGGGGGG‑CTTGTCTCCC‑AGTTTTTTTCCCC‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                      <  1:842044/146‑1 (MQ=255)                                                          tGATTCATTCCCG‑CCAACACCCGCATTTTT‑AT‑CCACCGCCTTGTCAAAAACTCCCTT‑CGCGGCGC‑T‑TTTCTGCC‑A‑CTTTTTTTCCC‑CCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                                     <  2:258494/150‑1 (MQ=255)                                                            aTTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACATTaaaacactcc                                                                                    >  2:704840/1‑141 (MQ=25)                                                                          accacccGCGCATTT‑AC‑CCACCCCCTTGTGAAAAAATCCCTG‑CGCGGCGG‑C‑TGTCTGCC‑ATTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                     <  2:504270/146‑1 (MQ=255)                                                                          aacaacCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAAAAATAACCTTAAGcaccgcat                                                                                   >  1:556694/1‑130 (MQ=25)                                                                           acaacCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAAATTAAAACCCTCCTGTTTTGTggggaac                                                                     >  2:700593/1‑145 (MQ=25)                                                                              cccGCGCTTTT‑AC‑CCACCCCCATGAAAAAAAATGCTTG‑CGGGGCGG‑C‑TGTCTCCC‑ATTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                 <  1:695074/150‑1 (MQ=255)                                                                               ccGCGCTTTT‑AC‑CCAACCCCTTGTAAAAAAATGCTTC‑CGGGGCGG‑C‑TGTCTGCA‑AGTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct                                                                  <  1:391260/149‑1 (MQ=255)                                                                                      ttt‑AC‑CCACCGCCTTGGGAAAAAATCCCTC‑GGGGGGGG‑C‑TGTCTCCC‑TGTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                          <  1:307163/150‑1 (MQ=255)                                                                                        ttaccccagcgctttgtgAAAAAATGCCTT‑CGGGGGGG‑C‑TGTCGGCA‑GGTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                           <  2:651532/135‑1 (MQ=255)                                                                                                    ctttgcaaaaaaTTCCTTGCGGGGCGG‑C‑TGTCTGCC‑ATTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                               <  1:815299/143‑1 (MQ=255)                                                                                                       gtgAAAAATTCCTTT‑CGGGGCGG‑C‑TGTCGGCA‑AGTTATTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                             <  2:782439/148‑1 (MQ=255)                                                                                                        tgaaaaaaTGCCTC‑CGCGGCGG‑C‑TGTCTGCC‑AGTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                          <  2:1102306/147‑1 (MQ=255)                                                                                                           aaaaaTGCCTT‑CGGGGCGG‑C‑TGTCTGCC‑AGCTATTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                      <  2:1099713/151‑1 (MQ=255)                                                                                                                 gCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTAAGCCCGCCCCGGTTTGTGGGAGCAGCCGCTGGGTACAAAGTTGGAAGGGGAGggcggc                                   >  2:66383/1‑148 (MQ=255)                                                                                                                       cgcgGGGG‑C‑TGTCTGCC‑ACTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                             <  2:720403/149‑1 (MQ=255)                                                                                                                          ggcgg‑c‑TTCCGGCC‑ATTTTTTTTCCCG‑CCGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg                            <  1:517257/147‑1 (MQ=255)                                                                                                                          ggcgg‑c‑TGTCTGCC‑AGTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                              <  1:519928/145‑1 (MQ=255)                                                                                                                          ggcgg‑c‑TGTCTCCC‑AGTTTTTTTCCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                        <  2:760038/151‑1 (MQ=255)                                                                                                                             ggggttGTCGGCC‑AGTTTTTTTCCCG‑CCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                       <  2:434461/146‑1 (MQ=255)                                                                                                                                   gTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGAAATAACTTTTAACCCCTCCTTTTTTTTGGGGATCGCTGCTGGGTAAAATATGGGAAGGGGAGGGGGGCGGATCAACTccccaaac                  >  2:962767/1‑144 (MQ=255)                                                                                                                                     tcccccggtttttttcCCC‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                  <  2:1131785/138‑1 (MQ=255)                                                                                                                                     cTGCC‑AGTTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                 <  1:1066648/149‑1 (MQ=255)                                                                                                                                      cgccctcattttgttCCG‑CCAA‑AAGAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                    <  2:336047/133‑1 (MQ=255)                                                                                                                                          cattttttttcCCG‑CCCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc              <  1:828011/142‑1 (MQ=255)                                                                                                                                           agtttttttcCCG‑CCAC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                  <  2:199134/138‑1 (MQ=255)                                                                                                                                               ttttttCCC‑CCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt         <  1:598987/148‑1 (MQ=255)                                                                                                                                                 tttcCCG‑CCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt     <  1:50477/150‑1 (MQ=255)                                                                                                                                                    tCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAATTTTAACCCCCCCTGTTTTGTGGCGATCGCTGCTGGGGACAATGTGGGAAGGGGAGTGGGGGGGAACAACTCCCCAAACGCGCTCATTTAg      >  2:670888/1‑146 (MQ=255)                                                                                                                                                    cccg‑ccgC‑AAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg  <  1:376342/149‑1 (MQ=255)                                                                                                                                                    cccg‑ccgA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt    <  1:177827/147‑1 (MQ=255)                                                                                                                                                                |                                                                                                                                            CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCG‑CCAACAACCGCGCATTT‑AT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGG‑C‑TGTCTGCC‑AGCTATTTTTCCG‑CCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357740‑358029

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