Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	51.6%	36.0 / 90.1	65	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (19/12);  new base A (0/33);  total (20/45)
Fisher's exact test for biased strand distribution p-value = 1.62e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-06

ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033                                                                                                                                                      |                                                                                                                                                  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)                                                                                                                                 ccagttttttccccgccaCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa                <  1:31447/139‑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)                                                                                                                                          ttcccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        <  2:76105/147‑1 (MQ=255)                                                                                                                                          tt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCACGAG‑AAAAAGATTTAAACCCGCCCCGTTGTGTGGGCATCGCCGGTGGGGAAAAAG‑TGGGAAGGGGGGGGGGGGGGATCAGCGCCCCAACCGCGCTAATTTAg          >  1:657167/1‑147 (MQ=255)                                                                                                                                           tcccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc       <  2:981555/148‑1 (MQ=255)                                                                                                                                             ccgccgCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg      <  2:228428/149‑1 (MQ=255)                                                                                                                                             ccgccaAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc    <  1:939078/151‑1 (MQ=255)                                                                                                                                                ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc  <  2:659672/150‑1 (MQ=255)                                                                                                                                                  gcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAG‑AAAAAACTTTAAAACCGGCATGTTTTGTGGGGAACGCTGCTGGGAACAAAG‑TGGGAAGGGGAGGGGGGCGGAACACCTCCCAAAACGCGCCCATTTAGTTTGCCTc  >  2:67624/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAT‑TTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357743‑358033

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