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

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357730‑358025                                                                                                                                                        |                                                                                                                                                     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt‑tttt                                                                                                                                                      <  1:1296992/151‑1 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt‑tttt                                                                                                                                                      >  1:349978/1‑151 (MQ=255)     gTAGGCCTGATAAGCGTTGCGCATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑ccg                                                                                                                                                  <  2:1078182/150‑1 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGcaaca                                                                                                                                          <  2:592377/149‑1 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGcaacaa                                                                                                                                         >  1:734887/1‑150 (MQ=255)              aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACaaa                                                                                                                                        <  1:1281186/151‑1 (MQ=255)                aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAAc                                                                                                                                      >  2:943476/1‑151 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCATACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAccc                                                                                                                                     >  1:805933/1‑151 (MQ=255)                 aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAAcc                                                                                                                                     <  2:1419831/151‑1 (MQ=255)                   cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAAccgc                                                                                                                                   >  1:568882/1‑151 (MQ=255)                   cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAAccg                                                                                                                                    <  2:1100193/150‑1 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                  >  2:206636/1‑151 (MQ=255)                         gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTtc                                                                                                                              >  2:1297746/1‑150 (MQ=255)                                   aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                    >  1:1034086/1‑150 (MQ=255)                                      ccACCCGCAGCCTTGTGTCAGCGGCTACCGTAATCCATTCCCCCCAACACCCGCCCTTTTACCCAACGCCATGTGCAAAAATCCCTTCGCGGCGGCTGTCGGCCAGCTTT‑TTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                <  2:963330/151‑1 (MQ=255)                                          ccGCAGCCCGGTGTCGGGGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATCCCTTCGCGGCGGCTGCCTGCCAGCTTT‑TTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                             <  2:9727/150‑1 (MQ=255)                                          ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggcg                                                                                                            >  1:441115/1‑149 (MQ=255)                                              aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                           >  1:691357/1‑147 (MQ=255)                                                 ccTGTGTCGGCGGCTCCCGTAATTCATTCCCCCCAACACCCGCGCATTTATCCACCGCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTT‑TTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                           <  2:761798/145‑1 (MQ=255)                                                   tgtgTCGGCGGCTACCGTTTTTCTTTCCCGCCAACACCCGCTATTTTACCCACCCCCATGTGAAAAAATGCTTTCGCGGCGGCTGTCTGCCAGTTTT‑TTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                          <  2:671532/144‑1 (MQ=255)                                                   tgggtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTATTTT‑CTGCCG‑GATACAACCTCCGCTCCCTTCAGCAGCGGCGTGGcccccgga                                                                                                    >  1:50112/4‑143 (MQ=14)                                                     tgtCGGCGGCCCCCGTTTTTCATTCCCGCCAACACCCCCGCATTTACCCACCGCCATGTGCAAAAATCCCTCCGCGGGGGCTGTCTCCCAGTTTT‑TTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                       <  2:774533/145‑1 (MQ=255)                                                      gtCAGCCGCTCCCTTGATTCTTTCCCCCCAACACCCCCCCTTTCACCCACCGCGATTCCAAAAAATCCCTCGGGGGCGCCTGTCGCCCTGTTTT‑TTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATa                                                                                                <  1:1379383/150‑1 (MQ=255)                                                         aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAG‑aaaaaaatt                                                                                            >  1:923980/1‑145 (MQ=255)                                                          ggggttcccttGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCCCAATTTGAAAAAATGCCTTGGCGGCGGCTTTCCGCCAGCTTT‑TTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                       <  2:30248/130‑1 (MQ=255)                                                          gcgcccccctgttttCTTTCCCCCCAACACCCGCCCTTTTACCCACCCCCTTTTCAAAAAGTGCCTCCGCGGGGGCGGTCTCCCTGTTTT‑TTCC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATa                                                                                                <  1:644287/133‑1 (MQ=255)                                                          gCGGCTCCCGTGATTCATTCCCCCCAACACCCGCCCTTTTACCCAACCCCTTGTCAAAAAATCCCTCCGCGGCGGCTGTCTCCCACTTAT‑TTTT‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAAt                                                                                                 <  2:1065892/145‑1 (MQ=255)                                                          gCGGCTACCGAGTTTCATTCCCCCCAACACCCCCGCTTTTACCCACCCCCATGAGCAAAAAGGCCTCGGCGGCGGTTTTCTGCCAGTTTT‑TTTC‑CCGCCAAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCttt                                                                                           <  1:153701/151‑1 (MQ=255)                                                           cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCACAAATGCCTTCGCGGCGGCTGTCTGCCAGCTCT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGCG‑GGATATCttt                                                                                           >  2:591839/1‑150 (MQ=255)                                                             gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcggg                                                                                                       >  1:1377573/1‑132 (MQ=255)                                                                 cggtttttCTTTCCCCCCAACACCCGCGTTTTTACCCACCGCCTTGTGAAAAAAGCCTTTGGGGGCGGTTGTCGGCCGGTTTT‑TTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                      <  2:564581/143‑1 (MQ=255)                                                                 ccGTGATTCTTTCCCCCCAACACCCGCGCATTTACCCACCCCCTTGTGAAAAAATCCTTCCGCGGCGGCTTTCTCCCAGTTTT‑TTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                      <  2:885229/149‑1 (MQ=255)                                                                 ccGTGATTCATTCCCGCAAACACCCGCGCATTTACCAACCCCCATTTAAAAAAAGGCCTGCGCGGGGGTTGTCTGCCAGTTTT‑TTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAg                                                                                        <  2:948807/147‑1 (MQ=255)                                                                   ttgttcccttcccGCCAACCCCCGCGATTTTCCCCACCCCCTTGTAAAAAATTGCCTCCGCGGGGGCTGCCGCCCGGTTTT‑TTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                      <  1:659285/136‑1 (MQ=255)                                                                   ttGATTCATTCCCCCCAACCCCCCCCCTTTTCCCCACCCCCATGTGAAAAAATCCCTTCGCGGCGGTTTCCTCCCATCTAT‑TTTC‑CCCCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                      <  1:354610/146‑1 (MQ=255)                                                                       ttCTTTCCCCCCAACACCCGCCTTTTTATCCACCGCCTTGTGAAAAAATGCCTTCGCGGCGGTTGTC‑GCCAGGTTT‑TTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                      <  1:726199/143‑1 (MQ=255)                                                                           ttccccccAACACCCGCCCATTTACCCACCCCCATGTGAAAAAAGGCCTCCGCGGCGGCTTTCTGCCAGTTAT‑TTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt                                                                              <  1:715639/141‑1 (MQ=255)                                                                             ccccccaccaCCCGCGCATTTACCCACCGCCGTGGGCAAAAATCCCTCCGCGGCGCCTGTCGCCCACTTAT‑TTTC‑CCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTtgtgt                                                                         <  1:934235/146‑1 (MQ=255)                                                                               ccccaccaCCCGCCCTTTTCCCCACCGCCGTGGGAAAAAATGCCTCGGGGGCGGCTGTCTGCCAGTTTT‑TTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                      <  1:572821/149‑1 (MQ=255)                                                                                          gcgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAACTTTTAACCCGCCCTGTTTTGTggggatag                                                                 >  2:189260/1‑139 (MQ=37)                                                                                                ttACCCACCCCCATGTAAAAAAATGCCTCCGCGGCGGCTGTCTCCCACCTTT‑TTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                      <  2:648829/150‑1 (MQ=255)                                                                                                     cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGAG‑AAAAAAAATTTAACCCGCCCTATTGTGTggggggcgcggct                                                            >  1:294597/1‑128 (MQ=21)                                                                                                      acccccATGTGAAAAAATCCCTCCGCGGCGGTTGTCTCCCATTTTT‑TTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑Tgg                                              <  2:35312/147‑1 (MQ=255)                                                                                                               taaaaaaaTCCTTGCGGGGGGGCTGTCTGCCAGTTTT‑TTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGaa                                      <  2:1191629/147‑1 (MQ=255)                                                                                                                  aaaagagCCTCCGCGGGGGCTTTCGGCCAGTTTT‑TTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGt                                    <  1:643199/143‑1 (MQ=255)                                                                                                                  aaaaaGGCCTCCGCGCCGGTTGCCTGCCATTTTT‑TTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTg                                   <  2:471829/150‑1 (MQ=255)                                                                                                                   aaaacccttccggggcgccTTTCGGCCAGTAAT‑TTTG‑CCCCAGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTgcg                                 <  1:985711/138‑1 (MQ=255)                                                                                                                   aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAGATTTAA‑CCCCCCCTGTTGTGTGGGGATCGCCGCCGGGTACAAAG‑TGgggaggggaggggg                                 >  2:131680/1‑147 (MQ=255)                                                                                                                        gCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGCG‑AAATAACTTTAAACCCGTCCTGTTTTGTGGGGATCGCTGCTGGTAACAATG‑TGGAAAGGGGAGGggggcg                              >  1:1041862/1‑149 (MQ=255)                                                                                                                         ccTTCGCGGCTGCTGTCTCCCAGCTAT‑TTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAAt                           <  2:39694/151‑1 (MQ=255)                                                                                                                          cTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑AAAAAAAATTTAAACCGCCCTGTTGTGTGGGGAACGCCGCTGGGAAAAATGTTGGAAAGGGagggggggg                               >  2:1412534/1‑139 (MQ=255)                                                                                                                          cTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAACCCCCCCTTGTT‑TGTGGGGATCGCCCCTGGTAACAAAA‑TGGGAAGGGGAgagggggg                              >  2:177824/1‑140 (MQ=255)                                                                                                                              ggggcggtTGTCTGCCATTTTT‑TTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCa                         <  2:994444/146‑1 (MQ=255)                                                                                                                                         cggcaactttt‑tttc‑CCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg           <  2:924138/141‑1 (MQ=255)                                                                                                                                         cTGCCAGTTTT‑TTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg             <  2:793042/149‑1 (MQ=255)                                                                                                                                            ccAGCTTT‑TTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                 <  2:735011/142‑1 (MQ=255)                                                                                                                                             cagtttt‑tttc‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa        <  2:350031/144‑1 (MQ=255)                                                                                                                                                tttt‑tttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                <  1:830751/136‑1 (MQ=255)                                                                                                                                                   t‑tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg    <  1:277514/148‑1 (MQ=255)                                                                                                                                                   t‑tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg  <  2:1149495/150‑1 (MQ=255)                                                                                                                                                   t‑tttc‑CCGCCCCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg    <  2:1140985/148‑1 (MQ=255)                                                                                                                                                   t‑tttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc          <  2:1185432/142‑1 (MQ=255)                                                                                                                                                       tt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCACGAG‑AAAAAGATTTAAACCCGCCCCGTTGTGTGGGCATCGCCGGTGGGGAAAAAG‑TGGGAAGGGGGGGGGGGGGGATCAGCGCCCCAACCGCGCTAATTTAg  >  1:657167/1‑147 (MQ=255)                                                                                                                                                        |                                                                                                                                                     ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357730‑358025

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