Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,953	0	A	C	41.8%	89.5 / 19.0	67	V661G (GTG→GGG)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (29/10);  new base C (0/28);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.15e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298805‑299103                                                                                                                                                          |                                                                                                                                                       ggTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CAcc                                                                                                                                                      >  1:669635/1‑151 (MQ=255) ggTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CAcc                                                                                                                                                      >  1:417793/1‑151 (MQ=255)          ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATAcccc                                                                                                                                             <  2:644099/151‑1 (MQ=255)            acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCg                                                                                                                                            <  1:1345842/150‑1 (MQ=255)            acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGc                                                                                                                                           <  1:377873/151‑1 (MQ=255)              acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGcc                                                                                                                                          >  1:1294220/1‑150 (MQ=255)                  cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATc                                                                                                                                     >  1:43259/1‑151 (MQ=255)                   gTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCg                                                                                                                                    >  1:402257/1‑151 (MQ=255)                    ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCa                                                                                                                                         >  2:1223203/1‑145 (MQ=255)                    ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTCTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCa                                                                                                                                         <  1:1222361/145‑1 (MQ=255)                          atggCCGCGCTCACGCCGCACAGGCCCAGCCCACCGACCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGTTGGGTATCACCCAGGAAAATCCCCTCCTGTTT‑GGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATCCCCCGCCATATCGTGAtt                                                                                                                               <  2:79895/146‑1 (MQ=255)                             gCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAAc                                                                                                                            <  2:988219/149‑1 (MQ=255)                             gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGCGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACg                                                                                                                           >  1:1069645/1‑150 (MQ=255)                              ccGCCCTCCCGCCGCACAGCCCCAGCTCCCCGCCCCCTTGGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTT‑TGGA‑ATATCCGC‑AGGAACCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                           <  1:882702/149‑1 (MQ=255)                               cccgcaggcccgccacgGGCCCACCTGCCCGACACCTTTGC‑CCTTCATCGGGAAGGTAAGGGGTTCGGTACAAACCGGAAAACTCCCCTCCTTTGG‑GGGG‑TTGCCCGT‑ATGCCCCGC‑CACCTAATACCCCGCCATATCGTGATTAACGAAg                                                                                                                        <  1:1016620/138‑1 (MQ=255)                                   cTCCCGCCGCACGGGCCCACCTCCCGGACCCTTTTGG‑CCTTCACCGGGGAGGATAGGGGGTCGGTTTCATCCAGAAAAATCCCCTCCTGTTT‑TGGG‑ATATCCGA‑AAGAACCGG‑CCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAg                                                                                                                     <  2:585918/150‑1 (MQ=255)                                     cACGCCGCACAGGCCCAGCTCACCGACCCCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTT‑TGGG‑ATATCCGC‑AGGAACCGG‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  <  2:984854/151‑1 (MQ=255)                                     cACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  >  1:636096/1‑151 (MQ=255)                                        cccgcATAAGCCCCGGCCACCGACACTTTCGG‑CCTTCGTGGGGGAGGCTGAGGTTGTGGAAACCCCCAGGAAAATCCCCTCCTTTTG‑TGAG‑AT‑TCCGCTAGGTACGGC‑CACAGAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  1:19106/148‑1 (MQ=255)                                          cgcACAGGCCCAGCCCCCCGACCCTTTGGG‑CCTTCATCGGGGAGGATAGGGGGTGGTTATCATCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGACCCGC‑CCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              <  1:1306045/150‑1 (MQ=255)                                             acaGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          >  1:928543/1‑151 (MQ=255)                                                ggCCCAGCTCCCCGACCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTT‑TGGG‑ATACCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                            <  2:1086778/146‑1 (MQ=255)                                                ggCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          <  2:420866/148‑1 (MQ=255)                                                  cccAGCCCCCCGCCCCCTTTGG‑CCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      <  2:916002/150‑1 (MQ=255)                                                  cccAGCCCCCCGCCACCTTTGG‑CCTTAAGGGGGGGGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTT‑TGGA‑AAACCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                     <  1:27753/151‑1 (MQ=255)                                                      gCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                                 >  1:751242/1‑151 (MQ=255)                                                         cACCGCCCCTTTTGC‑CCTTCATGGGGGAGGATAGGGGGTCGTTACCACCAAGAAAACTCCCCTCTTTTTT‑GGGT‑ATATCCGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                               <  2:265051/150‑1 (MQ=255)                                                           ccgcccccTTTGC‑CCTCCATGGGGGAGGATAGGGGGCATGTACCACCCAGGAAAATCCCCTCTTGTTT‑TGGG‑TTATCCGC‑ATGACCGGC‑CCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct                                                                                            <  1:651930/145‑1 (MQ=255)                                                            cGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                           >  1:121353/1‑151 (MQ=255)                                                                cccTTTGGCCCTTCAGGGGGGGGGATGGGGGGGCGTTATCACCCA‑GAAAAACCCCCCCTTTTT‑TGGGTATACCCGA‑AGGACCGGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                          <  1:73765/148‑1 (MQ=255)                                                                      gg‑CCTTCATGGGGGAGGAATGGGGGTAGATAGCACCCAGAAAAACCCCCTCCTGTTT‑GGGG‑ATACCCGA‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg                                                                                  <  1:1269810/150‑1 (MQ=255)                                                                       gcccTTCATGGGGGAGGATAGGGGGTCGTTATCACCCAGAAAAATCCCCCCCTGTTT‑GGGG‑ATACCCGC‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt                                                                                 <  2:838740/149‑1 (MQ=255)                                                                       gcccTTCAGCGGGGAGGATAGGGGTTCGTTATCCCCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATATCCGC‑ATGACCCCC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                   <  2:477981/147‑1 (MQ=255)                                                                           ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                             >  1:1273360/1‑151 (MQ=255)                                                                                 ggggAGGATAGGGGGTGGTTATCACCAAGGAAAATCCCCTCCTGTTT‑TGGG‑ATACCCGC‑ATAACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:409509/146‑1 (MQ=255)                                                                                 ggggAGAATTGGGGGTGGGTACCA‑CCAAGAAAACCCCCCCCTTTTTGTGGT‑ATATCCGG‑TTAACCCGG‑CCCCTAACCCCCCGCCATATCGTGATTAGCGAAGTAGCCCAAGCGGTCATCCTCCGCCAGCTCCTCCATCAGTGc                                                                               <  1:213657/143‑1 (MQ=255)                                                                                 ggggAGAATTGGGGGCCGGTTTCACCCGGAAAAACCCCCCCCTTTTG‑AGGG‑ATA‑CCGC‑AAGAACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:121372/146‑1 (MQ=255)                                                                                  gggAGGATAGGGGGTGGGTATCACCCAGAAAAATCCCCTCTTGTTT‑TGGG‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:468769/151‑1 (MQ=255)                                                                                  gggAGGATAGGGGGTCGGTATCACCCAGAAAAATCCCCCCCTGTTT‑TGGG‑ATATCCGC‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:1233507/151‑1 (MQ=255)                                                                                            ggggTGGGTATCACCCAGAAAAATCCCCTCCTTTTT‑TGGG‑ATACCCGC‑ATAACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                             <  2:1320131/150‑1 (MQ=255)                                                                                            ggggTCGTTATCACCGAGGAAAACCCCCCCTTTTTT‑TGGG‑ATATCCAC‑AGGACCGGC‑CCCCTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                <  2:94029/147‑1 (MQ=255)                                                                                            ggggTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            >  1:1368278/1‑151 (MQ=255)                                                                                                   gTACCACCCAGGAAACTCCCCTCCTGTTT‑GGGT‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  1:696477/151‑1 (MQ=255)                                                                                                    tATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                    >  2:1229410/1‑151 (MQ=255)                                                                                                          ccAGAAAAATCCCCTCTTTTTT‑TGGG‑ATACCCGC‑ATGACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  2:1065382/144‑1 (MQ=255)                                                                                                          ccAGAAAAATCCCCTCCTTTTT‑GGGG‑ATCTCCGC‑ATGACCCGC‑CACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                              <  2:87355/151‑1 (MQ=255)                                                                                                            aGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                            >  1:1178795/1‑151 (MQ=255)                                                                                                               aaaaTCCCCTCCTGTTT‑GGGG‑ATCCCCCC‑ATGACCCGC‑CCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:102854/144‑1 (MQ=255)                                                                                                                 aaccccctcctGTTT‑TGGG‑ATATCCGA‑ATGACCCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:528718/146‑1 (MQ=255)                                                                                                                  accccctcttGTTT‑TGGG‑ATACCCGC‑AGAACCCGG‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:585478/146‑1 (MQ=255)                                                                                                                   tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:1065351/1‑151 (MQ=255)                                                                                                                   tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  2:658472/1‑151 (MQ=255)                                                                                                                    cACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:581665/1‑150 (MQ=255)                                                                                                                     aCCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                   >  1:395882/1‑151 (MQ=255)                                                                                                                       ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:726682/1‑150 (MQ=255)                                                                                                                       ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  1:1240883/1‑150 (MQ=255)                                                                                                                       ctcctGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAATCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:690327/1‑150 (MQ=255)                                                                                                                           tGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  1:1381215/1‑151 (MQ=255)                                                                                                                             ttt‑tGTT‑ATCTCCGC‑TTGTACCGG‑CACCTTTCTCCCCGCCATATCGTGCGTGACGAGGTAACCCAGCCGGTCATCTACCGGCAGCTCCTGCAGCAGGGCGGCGGCCTTGCCCATAGTCATTGCGCCAAGGTCCTGGCGGCGCGCGGttt                             >  1:56818/1‑149 (MQ=255)                                                                                                                              tt‑tGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  1:159232/1‑151 (MQ=255)                                                                                                                                 tGGG‑ATATCCGC‑ATAAACCGG‑CACCTCATACCCCGCCATATCGTGCTTAACGAAGTAGACCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCAGGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGGTTTCgg                         >  1:492882/1‑150 (MQ=255)                                                                                                                                       tatCCGC‑ATGAACCGG‑CAACTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCACTGACCTGGCCGCGCGCAGTTTTCGGATTc                     >  1:376573/1‑149 (MQ=255)                                                                                                                                          ccGC‑ATGAACCGG‑CCCCTCATCCCCCGCCATATCGTGATCAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa                    <  2:392565/147‑1 (MQ=255)                                                                                                                                          ccGC‑ATGAACCGC‑CCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  <  2:395803/149‑1 (MQ=255)                                                                                                                                            gC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc              >  2:1044251/1‑151 (MQ=255)                                                                                                                                                   aCCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg        >  1:450130/1‑151 (MQ=255)                                                                                                                                                    ccGG‑CCCGTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt         <  1:1132032/149‑1 (MQ=255)                                                                                                                                                          aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  2:757136/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       GGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298805‑299103

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