Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,873	0	C	T	36.2%	110.4 / 36.2	81	S273N (AGC→AAC)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (35/16);  new base T (0/29);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 1.28e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAA  >  NC_000913/357724‑358011                                                                                                                                                           |                                                                                                                                                cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCt                                                                                                                                                >  1:13369/1‑151 (MQ=255)  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)                                                                                                                                                           |                                                                                                                                                CACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑G‑CC‑TT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAT‑GTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAA  >  NC_000913/357724‑358011

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