Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,591	0	A	.	35.8%	157.7 / 90.8	67	intergenic (+18/‑16)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (43/0);  new base . (22/2);  total (65/2)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947482‑2947740                                                                                                                |                                                                                                                                                      gtcatcatctcctGACTCTGGTTGGAGCAGCCTGGTAGCTCTGCGGTCTCATAACCCTAAGGTCGTCGGGTCTAATCATTCCCGCTCAACCAA‑TTATAATTTGATGAAGT‑AAAGCATTACTGTGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                                              <  4:218229/139‑1 (MQ=2)                                                                         aaaTCCGGCCCCCGCAACCAA‑TCA‑AATTTGAAGAAGTAAAAGCAGTAAGGTGAAGGGGGGTGGAGAAGACTGGAAGCTCGTCGCGCTCATATCACGAAGGTCGTAGGGTCAAATCCGGCCCGCGCAACCAATTATTGAACACCCTAACggg                                      >  4:742257/1‑151 (MQ=17)                                                                            tCCGGCCCCCGCAACCAA‑TGA‑ATTATGAAGATGTAAAAGCAGTACGGTGACGAGGTGTGCACAAGAATAGATGCACGTCGGGCACATAAACCGAACGTCGGAGGGACTATTCCGGCCCACGCAAAAATATATTGAACTCCCTaaatggttt                                   >  1:295761/1‑144 (MQ=255)                                                                                 cccccGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                              >  8:601623/1‑151 (MQ=255)                                                                                 cccccGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGtttttt                              >  2:1264330/1‑151 (MQ=255)                                                                                  ccccGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  5:927208/1‑151 (MQ=255)                                                                                  ccccGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCACCCGCACCCAATTATTGAACACCCTAACGGGTGttttttt                             >  1:1317374/1‑151 (MQ=255)                                                                                   cccGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTg                            >  8:1386464/1‑151 (MQ=255)                                                                                     cGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGtt                          >  3:436763/1‑151 (MQ=255)                                                                                     cGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGGTTTTTTGtt                          >  5:1252001/1‑151 (MQ=255)                                                                                      gCAACCAATTAA‑AATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTTAACACCCTCACgcctgagtatgtgttt                         >  7:562997/1‑136 (MQ=11)                                                                                      gCAACCAATTAA‑AATTTGATGAAGT‑AAAGCAGTAAGGTGAAGCGGGGTGGAGAAGCCAGGTAGCTCGTCGGGCTCATAACCCGAACGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACACCCTAACCCGTCtattttttta                          >  6:904506/1‑147 (MQ=11)                                                                                      gCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGttt                         >  6:857447/1‑151 (MQ=255)                                                                                      gCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGttt                         >  5:1483085/1‑151 (MQ=255)                                                                                       cAACCAATTAA‑AATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCAAATACGGTGTTTTATTGtttg                        >  1:1116716/1‑150 (MQ=11)                                                                                       cAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGtttg                        >  7:1390554/1‑150 (MQ=255)                                                                                       cAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTc                        >  7:541160/1‑151 (MQ=255)                                                                                        aaCCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  2:973270/1‑151 (MQ=255)                                                                                        aaCCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCAGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  4:449993/1‑151 (MQ=255)                                                                                         accaattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCTAACGGGTGTGTATTTGTTTCTg                      >  2:1023106/9‑151 (MQ=11)                                                                                          ccaattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGATCACCCACACGGCGGATTTTTGGTATCTgg                     >  1:1198720/8‑151 (MQ=1)                                                                                          ccaattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTATCGGGTGTTTTTTTGTATCggg                     >  5:254780/8‑151 (MQ=14)                                                                                           caattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTCACGGGGGTGTTTTTGTTTCTggg                    >  5:89032/7‑150 (MQ=14)                                                                                           caattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTggg                    >  5:980645/7‑150 (MQ=255)                                                                                           caattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGAAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTAACAACCTAAAAGGTCTTTTTTTTTTTCAGGt                    >  4:1452323/7‑151 (MQ=11)                                                                                           cAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  6:455343/1‑151 (MQ=255)                                                                                           cAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  6:400498/1‑151 (MQ=255)                                                                                           cAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  8:1342713/1‑151 (MQ=255)                                                                                           cAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  8:622188/1‑151 (MQ=255)                                                                                            aattaa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCAGGCTGTCGCGTGATTTGTTGTTTCTGGtg                   >  3:1351365/6‑150 (MQ=11)                                                                                            aa‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCACCGCAACCAATTATTGAACACCCTAACGGCTGTTTTTTTGTTTCTGGtc                   >  2:1186131/1‑151 (MQ=255)                                                                                            aa‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtg                   >  8:738214/1‑150 (MQ=255)                                                                                            aa‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  4:605059/1‑151 (MQ=255)                                                                                            aa‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  6:1140682/1‑151 (MQ=255)                                                                                             a‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtct                  >  5:1052025/1‑151 (MQ=255)                                                                                               taa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTCTCGGGGGTGTTGTTGTTTGTGGTCTcc                >  5:1103100/3‑151 (MQ=14)                                                                                               tCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc                 >  7:29503/1‑151 (MQ=255)                                                                                               tCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc                 >  6:723257/1‑151 (MQ=255)                                                                                                cA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTcc                >  7:1487213/1‑151 (MQ=255)                                                                                                aa‑aaTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAGCCTCACCGGGGTGTTTTTGGTTCTGGGCTccc               >  5:1043049/2‑151 (MQ=14)                                                                                                 a‑aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  3:1462943/1‑151 (MQ=255)                                                                                                   aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  2:593890/1‑150 (MQ=255)                                                                                                   aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  3:529008/1‑150 (MQ=255)                                                                                                    aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGCAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTATGTTTCTGGTCTCCCAt             >  2:249284/1‑151 (MQ=255)                                                                                                    aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa            <  2:978784/151‑1 (MQ=255)                                                                                                    aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAAAGGGTGTTTTTTTTTTTCTGGTCTCCCAt             >  3:900015/1‑150 (MQ=255)                                                                                                    aTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCAAAACCGTCTTTAATTTTATCGGGTCTCCCAt             >  8:229420/1‑150 (MQ=11)                                                                                                     tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  7:701398/1‑150 (MQ=255)                                                                                                     tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  7:42593/1‑150 (MQ=255)                                                                                                     tttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCTATAGGGGGATTTGTTGTATGTGGGCTCCCATa            >  7:255892/1‑150 (MQ=14)                                                                                                      ttGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATAGAACAACCAAACTCGGGTGTATTTTGTTCTGGTCTGCCATaa           >  5:1478025/1‑150 (MQ=11)                                                                                                       tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  6:1386412/1‑150 (MQ=255)                                                                                                       tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  2:1369496/1‑150 (MQ=255)                                                                                                       tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  5:112771/1‑150 (MQ=255)                                                                                                        gatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTTTTTGTTTCTGGTCTCCCATaaa          >  1:277624/1‑150 (MQ=255)                                                                                                        gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACCCTAACGCGTGTTTTTTTTTTAGTGGTCTCCCATaaaa         >  6:1295085/1‑150 (MQ=17)                                                                                                        gatgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAAAAACCTAACAAGTCTATAATTGTAACTGGTCTCCCaaaaaa         >  8:794585/1‑150 (MQ=11)                                                                                                         atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  1:879057/1‑150 (MQ=255)                                                                                                         atgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACAACCTAAAGGGTGTTTATTTGTATCTGGTCTCCCATaaaaac       >  6:1319574/1‑150 (MQ=17)                                                                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCGCCCGCAACCAATTATTGAACACCCAAAAAGGTGTTTTATTGGGGCGGGTCTCCCCTaaactcg      >  4:559818/1‑147 (MQ=14)                                                                                                          tgaAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACCCGTCTTTTATTGTTTCTGGTCTCCCATaaaaacg      >  6:1224423/1‑149 (MQ=17)                                                                                                           gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  5:1195370/1‑149 (MQ=255)                                                                                                             aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  5:111359/1‑151 (MQ=255)                                                                                                             aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  8:570893/1‑151 (MQ=255)                                                                                                             aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  3:416062/1‑151 (MQ=255)                                                                                                              gTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcgc   >  8:476483/1‑151 (MQ=255)                                                                                                               tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  8:825393/1‑151 (MQ=255)                                                                                                                aaaaGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  6:129031/1‑150 (MQ=255)                                                                                                                |                                                                                                                                                      AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCA‑AATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947482‑2947740

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