Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	41.9%	94.4 / 25.2	75	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (31/12);  new base C (0/31);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 2.90e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357730‑358029                                                                                                                                                            |                                                                                                                                                           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)                                                                                                                                             tctgcctctttttttc‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                    <  1:646660/136‑1 (MQ=255)                                                                                                                                             tCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGA‑G‑AAAAAAATTTTAAACCCCCCT‑GTTGTGTGGGGATAGCCGCTGGGGAAAATATGGGAAAGGGAGGGGGGGGGATCAGCTCCCGaa                    >  1:972720/1‑147 (MQ=255)                                                                                                                                                   gttttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                  <  1:1103892/138‑1 (MQ=255)                                                                                                                                                      tATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt          <  2:466298/148‑1 (MQ=255)                                                                                                                                                            t‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGAGG‑AAAAAAATTTAA‑CCCGGCAT‑GTTGTGGGGGAAAAGCTGGTGGGAAAAAAATGGGAAAGGGAGAGGGGGGGAACAACCGCACAAACCCGCTAATTGAGtttg  >  1:601971/1‑148 (MQ=255)                                                                                                                                                            |                                                                                                                                                           ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357730‑358029

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