Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,875	0	A	T	40.3%	103.6 / 33.8	79	N272K (AAT→AAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (33/13);  new base T (0/31);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT  >  NC_000913/357725‑358021                                                                                                                                                            |                                                                                                                                                        aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTa                                                                                                                                                         <  1:191766/151‑1 (MQ=255)  cTGAACGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTa                                                                                                                                                         <  2:15998/150‑1 (MQ=255)   tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                                       >  2:431460/1‑151 (MQ=255)    gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                      >  1:856842/1‑151 (MQ=255)    gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                                       <  1:564417/150‑1 (MQ=255)      aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                    >  1:1019839/1‑151 (MQ=255)         cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccg                                                                                                                                                <  2:1283020/151‑1 (MQ=255)         cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                 <  1:281065/150‑1 (MQ=255)          gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgc                                                                                                                                               <  2:984708/151‑1 (MQ=255)           tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgcc                                                                                                                                              >  2:38858/1‑151 (MQ=255)            aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgcc                                                                                                                                              <  1:1139329/150‑1 (MQ=255)             ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑ccgccg                                                                                                                                             <  1:473598/150‑1 (MQ=255)             ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                            >  2:1103237/1‑151 (MQ=255)               ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGca                                                                                                                                           >  2:1233317/1‑150 (MQ=255)                  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaa                                                                                                                                          <  2:38067/148‑1 (MQ=255)                  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACAGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                        >  2:694653/1‑150 (MQ=255)                     aaGCGTAGCGCATCATGAAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc                                                                                                                                    <  2:1141729/151‑1 (MQ=255)                      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc                                                                                                                                   >  2:1009346/1‑151 (MQ=255)                      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc                                                                                                                                   <  1:41165/151‑1 (MQ=255)                        cgcaGCCCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTCCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGCCTTCCAGCAATTTTT‑CCTCCCCACCAAAAccg                                                                                                                                  >  2:938793/4‑150 (MQ=17)                         gTAGCGCACCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGTCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                <  2:1054173/151‑1 (MQ=255)                              gcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCAGCGGCTGTCTGCCAGCTGTAGTT‑TCTCCGCAACACAAcggcaactc                                                                                                                            >  2:274509/1‑145 (MQ=14)                               cATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCAGCGGCTGTCTGCCAGCTATTTTT‑TCGCCGGAATACAAccagaacg                                                                                                                             >  2:41159/1‑142 (MQ=11)                                aTCAGGCAATTCCAGCCGCAGACCTGTGTTAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                                         >  1:632759/1‑151 (MQ=255)                                aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctc                                                                                                                          >  1:1075833/1‑150 (MQ=255)                                  cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                       >  1:100978/1‑151 (MQ=255)                                   aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                       >  1:305618/1‑150 (MQ=255)                                       aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                    >  1:478952/1‑149 (MQ=255)                                           ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                              >  1:718399/1‑151 (MQ=255)                                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                            >  1:335326/1‑151 (MQ=255)                                               ccGCAGACTTGTGTCAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCATTTGCAAAAAT‑C‑CC‑TT‑C‑GGGGCGGCTGTCTCCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggg                                                                                                          <  1:140607/150‑1 (MQ=255)                                                 gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                         >  1:965456/1‑149 (MQ=255)                                                     cccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCTTTTACCCACCGCCATGTGAAAAAAT‑G‑CC‑TC‑C‑GGGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga                                                                                                    <  1:804919/149‑1 (MQ=255)                                                       gtgtgTCAGCGCTTCCGGGGATTCTTTCCCCCCAACCCCCAAGCATTTACCAAACCCTATGACCAAAAAT‑‑‑CC‑TT‑CGGCGGCGTCTGTTTGCCGTTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                         <  2:349882/142‑1 (MQ=255)                                                        tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                         >  1:598246/1‑142 (MQ=255)                                                         gtgtCGGCGGCTACCGTGATTCTTTCCCGCAAACACCCGCGCATTTACCCAACGCCATGTGAAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑aaa                                                                                              <  1:680265/150‑1 (MQ=255)                                                          tgccggCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGATTTTACCCACCGCCTTGTAAAAAAAT‑C‑CC‑TT‑C‑GCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAAt                                                                                             <  1:702814/145‑1 (MQ=255)                                                          tgccAGCGCCTCCTGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATTTGAAAAAAT‑C‑CC‑TG‑C‑GCGGGGGTTGTCTGCCGGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAga‑g‑a                                                                                                <  1:399420/144‑1 (MQ=255)                                                              ggcgctaccccgttttCATTCCCCCAAACACCCGC‑CATTTCTCCCACGCCTTGAGAAAAAAA‑G‑CCTTT‑C‑GCGGGGGCTGCCTGAAAGTTTTTTTT‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAg                                                                                           <  2:576318/134‑1 (MQ=255)                                                               gCGGCTCCCGTGATTCATTCCCCCCAACACCCGCTCATTTACCCACCCCCATGTAAAAAAAT‑G‑CT‑TC‑C‑GGGGGGGCTTTCGCCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAAt                                                                                             <  1:431867/145‑1 (MQ=255)                                                                 gcctcctgtgtTTCATTTCCGCAAACACCCGCGCATTTACCCACCGCCATGTGAAAAAAT‑G‑CC‑TC‑C‑GGGGCGTTTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTaa                                                                                     <  1:381073/144‑1 (MQ=255)                                                                     aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGG‑AAAAAACTTTAAccc                                                                                  >  1:427304/1‑151 (MQ=25)                                                                      ccttgtttCTTTCCCCCCAACACCCGCGCATTCCCCCACCGCCATGTGCAAAATG‑C‑CC‑TC‑C‑GGGGCGGTGGTCTGCCGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGc                                                                                   <  2:1252020/142‑1 (MQ=255)                                                                      ccgtatttCATTCCCCCAAAAACCCGCGCTTTTATCCACCGCCATGTAAAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAg                                                                                    <  1:263594/141‑1 (MQ=255)                                                                      ccgtattcaaTTCCCGCCAACACCCGCCCTTTTACCCACCGCCTGGTG‑AAAAAA‑G‑CC‑TTCC‑GCGGCGGCGCTCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGccc                                                                                 <  2:174607/141‑1 (MQ=255)                                                                      ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑G‑aaaaaaactttaaacc                                                                                 >  1:461433/1‑137 (MQ=25)                                                                      ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGACCAGA‑a‑aaagaagattcaccgc                                                                                 >  1:1027173/1‑137 (MQ=25)                                                                           aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGG‑A‑AAAAAGCTTTa                                                                                      >  2:422978/1‑141 (MQ=25)                                                                              cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑G‑AAAAAACTTTaaaaccccccagtt                                                                         >  1:715327/1‑138 (MQ=25)                                                                                ttcccGCAAACACCCCCGCTTTTCCCCACCGCATTGTAAAAAAAT‑G‑CC‑TC‑C‑GCGGGGGCTTTCTCCCATTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAt‑gt                                                                         <  1:617610/144‑1 (MQ=255)                                                                                  cccccaaacaccCGCGCTTTCCCCCACCGCCATGGGAAAAAAG‑C‑CC‑TC‑C‑GGGGCGGCTGTCGGCAATTGTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTtgtgtg                                                                   <  2:333896/145‑1 (MQ=255)                                                                                  cccGCCAACACCCGCCTTTTTACCCACCCCCTTGTGAAAAAAT‑G‑CC‑TT‑C‑GCGGCGGTTGTCTGCCAGTTTTTTTC‑CCGCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTtgtg                                                                     <  2:457935/148‑1 (MQ=255)                                                                                   cccccAACAACCCCGTATTTACCCAACGCCATTTGCAAAAAT‑G‑CC‑TG‑C‑GCGGGGGCTTTCTGCCGGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTtgtgt                                                                    <  1:1254925/146‑1 (MQ=255)                                                                                         acaccccctcttttACCCACCCCCATGTAAAAAAAG‑C‑CC‑TC‑C‑GGGGGGTCTGTCTGCCTTTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATc                                                             <  1:1089607/139‑1 (MQ=255)                                                                                            acCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑G‑AAAAAACTTTTATCCCCCCCT‑GTTTTGTGGGGATCg                                                            >  1:521794/1‑148 (MQ=25)                                                                                            acCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAAAAACTTTAAAACCGTCAT‑GTTGTGTGGGAATc                                                             >  1:636108/1‑147 (MQ=255)                                                                                              cgcgCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGG‑G‑GAAAAAATTTATACCCGGCAT‑CTTGTGTggggagtggtgc                                                        >  2:146945/1‑140 (MQ=25)                                                                                                 gcttttACCCACCGCCATGTGAAAAAAT‑G‑CC‑TC‑C‑GGGGCGGTTTTCTGCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGctgctg                                                      <  1:975438/146‑1 (MQ=255)                                                                                                  cATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑G‑AAAAAACTTTAAAACCCCCCT‑GTTTTGTGGGAAACg                                                            >  1:1011957/1‑142 (MQ=37)                                                                                                    tttcccaacccccgtgtgAAAAATT‑C‑CC‑TG‑C‑GGGGCGGTGTTCGGCCAGTTATTTGC‑CCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGt                                                    <  2:1231295/134‑1 (MQ=255)                                                                                                    tttacccaccccattgtgAAAAAAT‑G‑CC‑TT‑C‑GGGGCGTCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTa                                                   <  1:880875/135‑1 (MQ=255)                                                                                                    tttATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGG‑G‑GAAAAACATTAAGCCCATCCTGGTTGTG‑GGGGATCCCTGCTgg                                                     >  2:241619/1‑147 (MQ=25)                                                                                                    tttACCCACCGCCTTGTAAAAAATT‑C‑CT‑TT‑G‑GGGGCGGCTGTCTGCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTgg                                                     <  1:458797/147‑1 (MQ=255)                                                                                                    tttACCCAACTCCTGGTGAAAAAAG‑C‑CC‑TT‑G‑GGGGCGCCTGCCTGCAAGTTTTTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTaa                                                  <  1:691768/150‑1 (MQ=255)                                                                                                         caagcgcttgggcaaaaaaT‑C‑CC‑TC‑C‑GCGGCGTCTTTCT‑CCGGGTGTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GTGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAAc                                                 <  2:20796/132‑1 (MQ=255)                                                                                                           accGCTATGTGAAAAAAT‑G‑CC‑TG‑C‑GCGGGGGCTTTCTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAAtg                                             <  2:13366/146‑1 (MQ=255)                                                                                                           accGCCATGTGAAAAAAT‑C‑CC‑TC‑C‑GCGGCGGCTGTCTCCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                            <  1:366642/147‑1 (MQ=255)                                                                                                             ccccATGT‑CAAAAATGG‑CC‑TG‑C‑GGGGCGGTTGTCTGCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                          <  1:682148/147‑1 (MQ=255)                                                                                                               cccttggaaaaaaT‑GCCC‑TT‑C‑GGGCCGGCTGTCTCCCTGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                          <  2:434569/140‑1 (MQ=255)                                                                                                                  tgtgAAAAAAT‑G‑CC‑TG‑C‑GGGGGGGTTGTCGCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                    <  1:725379/150‑1 (MQ=255)                                                                                                                  tggcaaaaaaT‑G‑CC‑TC‑C‑GGGGCGGCTTTCTCCCATTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                     <  1:274649/144‑1 (MQ=255)                                                                                                                   gtgCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑G‑AAAAAAAATTTAAACCGCCCT‑GTTGTGTGGGGATCCGTGCTggggaaaaagtg                                           >  1:388109/1‑132 (MQ=255)                                                                                                                      cAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑G‑AAAAAAATTTTAAACCCTCCC‑GTTTTGTGGGAAGCGCCTCTGGGTAAAAAGTGGGAAggg                                    >  1:339747/1‑146 (MQ=255)                                                                                                                      aaaaaaT‑C‑CT‑TG‑C‑GGGGGGGCTGTCGGCTTGTTTTTTCC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAAc                                                 <  2:536214/132‑1 (MQ=255)                                                                                                                       aaaaaT‑C‑CC‑TC‑C‑GCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                               <  1:186414/150‑1 (MQ=255)                                                                                                                       aaaaaG‑G‑CC‑TC‑G‑GGGGTGGCTGTTTGCGAGTTTTTTTG‑TCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                    <  2:1052209/145‑1 (MQ=255)                                                                                                                              g‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAACCTTTAAACCC‑CCCT‑GTTTTGTGGGGATCGCTGCTGGGGACAAAGTGGGAAGGGGAgggggggg                          >  1:540473/1‑142 (MQ=255)                                                                                                                                          ggcggcTGTCTGCCAGCTATTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                <  1:1239325/151‑1 (MQ=255)                                                                                                                                             ggttGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg              <  1:695138/147‑1 (MQ=255)                                                                                                                                                  tCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGA‑G‑AAAAAAATTTTAAACCCCCCT‑GTTGTGTGGGGATAGCCGCTGGGGAAAATATGGGAAAGGGAGGGGGGGGGATCAGCTCCCGaa            >  1:972720/1‑147 (MQ=255)                                                                                                                                                           tATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt  <  2:466298/148‑1 (MQ=255)                                                                                                                                                            |                                                                                                                                                        ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAT  >  NC_000913/357725‑358021

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