Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,294	0	T	C	37.9%	85.6 / 24.8	68	intergenic (+100/+129)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (12/29);  new base C (25/0);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 2.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423146‑3423443                                                                                                                                                       |                                                                                                                                                      tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCCG‑Ccccccccccccc                                                                                                                                                     >  1:118310/1‑151 (MQ=255)   cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑Ccccccccccccccg                                                                                                                                                   >  1:710199/1‑149 (MQ=255)    tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGCCTG‑CCCGCCCCCCCCCcgg                                                                                                                                                  >  2:891601/1‑151 (MQ=255)     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCggg                                                                                                                                                  <  1:1322423/150‑1 (MQ=255)     aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGc                                                                                                                                                 <  1:43154/151‑1 (MQ=255)        tCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGTCCGT‑CCCCCCCCCCCCCCGGcc                                                                                                                                                >  1:8006/1‑149 (MQ=255)        tCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCCG‑CCCCCCCCCCCCCCGGcc                                                                                                                                                >  1:486339/1‑149 (MQ=255)         cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCCCCCCCCCCCCGGcc                                                                                                                                                >  1:1274303/1‑148 (MQ=255)                tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCCT‑CCCCCCCCCCCCCCGGCCGTTTCtt                                                                                                                                         >  2:1233267/1‑148 (MQ=255)                  cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGCTCCCT‑CCCCCCACCCCCCcgggc                                                                                                                                                >  2:589192/1‑137 (MQ=255)                  cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCCG‑CCCCCCCCCCCCCCGGCCggt                                                                                                                                             >  1:530159/1‑140 (MQ=37)                   gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGTCCTG‑CCCCCCCCCCCCCCGGCCGGTGCtt                                                                                                                                         >  1:334738/1‑145 (MQ=255)                       gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGTTGTCCTG‑CCCGCCACCCTCCCGGGCGGTGCttttcaa                                                                                                                                    >  1:1161632/1‑141 (MQ=255)                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTT‑CCCGCCACCCTCCCGGCCCGTGCTTCGCCACGtt                                                                                                                                >  1:930686/1‑149 (MQ=255)                        cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCCCCCCCCCCCGGGCCGTTGGTTTGCaaactt                                                                                                                                >  2:1178032/1‑145 (MQ=25)                         aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCCCCACCCTCCCGGCCGGTTCTTTGCAACGtt                                                                                                                                >  2:37384/1‑148 (MQ=255)                           ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTGTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaaa                                                                                                                            <  2:1052345/150‑1 (MQ=255)                               aTCGCAGTTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                        <  2:318741/150‑1 (MQ=255)                               aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGGTTCTTTCCCACGTTCAAATTCGc                                                                                                                       >  1:142596/1‑151 (MQ=21)                               aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCGCCCCCCTCCGCCCCGCTGCTTCTCAACGTTCAAAAcccc                                                                                                                       >  2:840833/1‑149 (MQ=255)                                 cGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                     <  1:902262/151‑1 (MQ=255)                                    aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                  <  1:530172/151‑1 (MQ=255)                                      gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                <  1:785082/151‑1 (MQ=255)                                       tAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGCCTT‑TCCGCCCACCTCCGGGCCGTTGGTTTTCAACGTTCAAATACACTCCcg                                                                                                                  >  2:629987/1‑148 (MQ=255)                                           ccATTAATTGAATGTTAGTTCTAAAAGCAAAAGGGCCATCCTTTCGGAGGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                            <  2:410846/150‑1 (MQ=255)                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                           <  1:626153/151‑1 (MQ=255)                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                           <  2:1320703/151‑1 (MQ=255)                                           ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                            <  2:527340/150‑1 (MQ=255)                                            cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTT‑GCCGCCACCCTCCGGGCCGGTTGTTCGCAACCTTCCCATACCCCCCCGCCGGagt                                                                                                           >  2:1212537/1‑148 (MQ=18)                                                    gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCCGC‑CCCCCCCCCCCCCGGGCGGTTGC‑TCTCGAAGTTCAAATTCGCTCCCCGCGGGGTTGtgcct                                                                                                    >  2:147622/1‑145 (MQ=255)                                                     aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGCTTCTTCTCCACATTCAAAACCCCCCCCCGCGGGttt                                                                                                          >  1:870260/1‑142 (MQ=18)                                                      aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                <  1:1033429/151‑1 (MQ=255)                                                          tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                             <  1:779046/150‑1 (MQ=255)                                                                gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGCCTT‑CCCGCCCCCCCCCCGGGCGGTTTTTCTCAAAGTTCAAATCCCCCCCCCCCGGGATTGTGCTTCTCcgggg                                                                                            >  2:1217259/1‑143 (MQ=255)                                                                gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATCGC‑GGGGGTCCCG‑CCCGCCACCCCTCGCGCCGCTTCCTCCCTACGCTCAAATACACTccc                                                                                                                   >  2:1401013/1‑122 (MQ=255)                                                                 aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTC‑CCCCCCCCCCTCCCGGCCGGTGCCTTGCAACGCTCAAATCCGCCCCCCGCGGGGTTTTCCCTCTCTGGagag                                                                                          >  1:397833/1‑146 (MQ=17)                                                                  aatgCAAAAAGGCCATCCGTCAGGATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCCCCCCCCCGGCCCGTGGTTTGCCACGTTCAAATCCGCCCCCCGCGGGTTTTTTCTtctccagcgagcg                                                                                        >  1:602561/4‑138 (MQ=11)                                                                    aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCGCCACCCTCCCGGCCGGTGCTTCTCCACGTTTAAATCCCCTCCCCGCGGGTTTTTTCTACTCCGGGGAGTGTTcctc                                                                                  >  2:415761/1‑148 (MQ=255)                                                                       aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCCCCCCCCCCCCCGCCGGTGCTTTGCCACGTTTAAATCCCCTCCCCGCGGGTTTGTCCTTCTCCGGGGAGGGttgtccgcc                                                                               >  2:115769/1‑144 (MQ=255)                                                                        aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGCCTG‑GCCGCCCCCCTCCCGGGCGGTGCTTCCCCAAGTTCAAATACCCCCCCCCCGGGTTTGTTCTTCTTTGGGGAGGGTTTACCgcc                                                                               >  2:1289146/1‑148 (MQ=255)                                                                         aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAGGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                               <  2:1014975/101‑1 (MQ=255)                                                                           aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCCG‑CCCCCCCCCCCCCGGGCCGGTTCTTCGCAACGTTCAAATACCCCCCCCGCGGGTTTGTTCTACTCCGGAGAGAGTGTACCGccaaa                                                                            >  2:1199716/1‑150 (MQ=14)                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                          <  2:249231/151‑1 (MQ=255)                                                                            ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCCCCCCCCCCCCGGCCGTTGCTTCGCAACCATCAAAACCGCCCCCCGCGGCGTTTTTCTACTCTGGGGAGGGTTTCCCCacaaa                                                                            >  1:77380/1‑149 (MQ=255)                                                                                  ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGAt                                                                    <  1:1261766/151‑1 (MQ=255)                                                                                  ccTTTCGGATGGCCTTTCGCTGGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCGCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGAt                                                                    <  1:279717/151‑1 (MQ=255)                                                                                   cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGCCTG‑CCCCCCACCCCCCCGGCCGGTTCTTCGCCACGTTCCAATACCCCCCCCGCGGGTTTGTTCTACTTCGGGGAGTGTTTACCCACAAaaaacaaa                                                                     >  1:1230452/1‑147 (MQ=255)                                                                                    tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGTCCCG‑CCCCCCCCCCCCCCGGCCGTTGCTTTGCAACCCTTAAATCCCCCCCCCGCGGGTTTGTCCTACTCTGGGGAGTGTTCACCCCCAAACAAAAGATa                                                                   >  1:551034/1‑150 (MQ=255)                                                                                    tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCCCCCCGCCGTTGCTTCGCAAAGTTTTAATCCGCGCCCGGCGGATTTGTCCTACCCGGGAGAGTGTTCTCCCACAAACAACAGATa                                                                   >  1:127535/1‑150 (MQ=14)                                                                                       cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTT‑CCCGCCACCCCCCGCGCCCTTGCTTCTCAACGTTTCAATCCGCCCCCCGCGGGTTTTTTCTTCTCGGGGGAGGGTTCACCGACAAACAACAAATaaaa                                                                >  1:1208991/1‑150 (MQ=14)                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACa                                                              <  1:933536/151‑1 (MQ=255)                                                                                        ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGGTTTGTCCTACTCGGGGGAGGGTTCACCGACAAACAACAGATaaaa                                                                >  1:948545/1‑149 (MQ=255)                                                                                          aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                            <  2:869999/151‑1 (MQ=255)                                                                                              ccTTTCGCTGGATTGGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGc                                                        <  2:1394525/151‑1 (MQ=255)                                                                                                 ttCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGGGGGCTG‑CCCCCCCCCCACCCGGCCGCTGGTTTGCAACGTTCAAATCCCCCCCCCGCGGGGTTTTCCTTCTCGCGGGAGAGTTTCCCCACAAAAAACAGAGAAAACAAAAGGacc                                                      >  2:653364/1‑150 (MQ=255)                                                                                                   cGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                    <  2:791652/150‑1 (MQ=255)                                                                                                       tgatttgatGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                                 <  1:1172088/149‑1 (MQ=40)                                                                                                           ttgatGTCTGGCAGTTTATGGC‑GGGCGGCCCGCCCCCCCCCCCCCCGGGCGGTGGTTTCGAAAGTGTCACA‑CCGCCCCCCGCCGGTTTGTTCTACCCGGGAGAGGGGTCCCCCCCACACCACAaaaaaaaaaaaaagacc                                                      >  2:560847/1‑135 (MQ=255)                                                                                                            tgatGTCTGGCAGTTTATGGC‑GGGCGGCCTG‑CCCGCCCCCCCCCCGGGCGTTGCTTCGCAACGTTCAAATCCCCCCCCCGCGGGATTTTCCTACTCCGGGGAGGGTTTACCGCCAAACAAAAAATAAAAAAAAAAGGCCAGTCTTcc                                              >  1:1050004/1‑147 (MQ=2)                                                                                                                gTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAg                                       <  2:1102389/150‑1 (MQ=38)                                                                                                                gTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                      <  2:709992/151‑1 (MQ=37)                                                                                                                   tGGCAGTTTATGGC‑GGGCGTGCCG‑CCCGCCCCCCCCCCGGCCCTGGTTTTTCGACGTTCAAATACCCTCCCCGCGGGCTTGTTCTTCTCGGGGGAGTGTTTACCCACAAACAACAGAAAAAAAAAAAAACCCCGTCCTCCGACCGAGccct                                   >  2:727760/1‑149 (MQ=2)                                                                                                                    ggCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                   <  2:486171/150‑1 (MQ=35)                                                                                                                                c‑gggcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                       >  1:236274/1‑150 (MQ=38)                                                                                                                                  gggcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                     <  2:535462/151‑1 (MQ=38)                                                                                                                                      gTCCTG‑CCCGCCACCCTCCTGGCCGTTGCTTCGCATCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGACAGTGTTCACCGACAAACAACAGATGAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTTTTTTATTTGagg                    <  2:832582/148‑3 (MQ=14)                                                                                                                                      gTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                 >  1:974726/1‑151 (MQ=40)                                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc       <  1:559832/151‑1 (MQ=40)                                                                                                                                                     ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct   <  2:871150/151‑1 (MQ=40)                                                                                                                                                     ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct   <  2:877885/151‑1 (MQ=40)                                                                                                                                                       tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc  >  2:609893/1‑150 (MQ=40)                                                                                                                                                       |                                                                                                                                                      TCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGC‑GGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTC  >  NC_000913/3423146‑3423443

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