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

CCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCG  >  NC_000913/298771‑299070                                                                                                                                                       |                                                                                                                                                           ccAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCa                                                                                                                                                            <  1:701638/150‑1 (MQ=255) ccAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAc                                                                                                                                                           >  1:1360595/1‑151 (MQ=255)  cAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcc                                                                                                                                                          >  2:73778/1‑151 (MQ=255)   agagTGAGGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcct                                                                                                                                                         >  2:306106/1‑151 (MQ=255)    gagTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAcctc                                                                                                                                                        >  1:392449/1‑151 (MQ=255)     agTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATAGGGTCGGTATCATCCAGGAAAATCAcctcc                                                                                                                                                       <  1:1353566/151‑1 (MQ=255)        gATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGt                                                                                                                                                    <  1:14057/151‑1 (MQ=255)            ggATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑Tg                                                                                                                                               <  1:105771/151‑1 (MQ=255)             gATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑Tgg                                                                                                                                              <  2:320688/151‑1 (MQ=255)               tAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑Tggg                                                                                                                                             >  1:838080/1‑150 (MQ=255)                     cGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCg                                                                                                                                     <  2:306088/151‑1 (MQ=255)                        aCCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑At                                                                                                                                 <  2:190426/151‑1 (MQ=255)                         cccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATg                                                                                                                                <  1:906376/151‑1 (MQ=255)                            gAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAAc                                                                                                                             <  2:454456/151‑1 (MQ=255)                             aaTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAAcc                                                                                                                            <  2:29362/151‑1 (MQ=255)                             aaTACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAAcc                                                                                                                            <  1:979554/151‑1 (MQ=255)                                 ccGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑c                                                                                                                        >  1:395687/1‑150 (MQ=255)                                   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)                                                                                                                                                     tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg    >  2:658472/1‑151 (MQ=255)                                                                                                                                                     tCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg    >  1:1065351/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)                                                                                                                                                       |                                                                                                                                                           CCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGG‑CCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTT‑TGGG‑ATATCCGC‑ATGAACCGG‑CACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCG  >  NC_000913/298771‑299070

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