Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,625	0	A	C	40.0%	178.0 / 17.9	95	V105G (GTC→GGC)	yfgF	cyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base A (48/9);  new base C (0/38);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 2.80e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATC  >  NC_000913/2628481‑2628771                                                                                                                                                    |                                                                                                                                                    tctcTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTg                                                                                                                                               <  3:705222/150‑1 (MQ=255)  ctctGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGt                                                                                                                                              >  1:1196499/1‑150 (MQ=255)     tGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAAcc                                                                                                                                          >  3:34065/1‑151 (MQ=255)        tCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGcc                                                                                                                                       >  2:410851/1‑151 (MQ=255)         cGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGccc                                                                                                                                      >  3:1418624/1‑151 (MQ=255)          gCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGcccc                                                                                                                                     >  3:1096862/1‑151 (MQ=255)            cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCa                                                                                                                                    >  2:23812/1‑150 (MQ=255)             aGTAATCCTACAAATGCAGCAACCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCGGGAAAATTAACCGGTTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGC‑CCCTGTAACCGCCCCaa                                                                                                                                   <  1:1267161/150‑1 (MQ=255)             aGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGACA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAAc                                                                                                                                  >  3:1338213/1‑151 (MQ=255)             aGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAAc                                                                                                                                  >  3:748590/1‑151 (MQ=255)               tAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAAc                                                                                                                                  >  1:1167397/1‑149 (MQ=255)                 aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACag                                                                                                                                >  1:1324549/1‑149 (MQ=255)                  tCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAAAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACag                                                                                                                                >  1:103792/1‑148 (MQ=255)                   ccTACAAATGCAGCAAACGGAAAAAGTATCAAAACCAGCGTTGCAGAAAA‑CACAATCTGCCAGAATATCCGCTGGGAAATTAAACGGGTATCCCCAT‑GGAAAA‑CGTTGTTCCGCGGGGAAGTAAAGC‑CCCTGTAACCGCCCCAACagag                                                                                                                              <  1:236314/149‑1 (MQ=255)                   ccTAAAAATGCAACAAACTGAAAAAGTATCAGAAACAGCGTTGCGGAAAC‑CAAAATCTGCCAAAATATCCGCGGGAAAATTAAACGGTTATCCCCAT‑GTAAAC‑CTTTGTTCCGCCGGGTAGTAAAGC‑CCCTGTACCCGCCCCAACagag                                                                                                                              <  1:1110316/149‑1 (MQ=255)                      aaaaaTGCAGAAAACTGGAAAAGTATCAGAAACAGCGTGGCGGGAAA‑CACAATCTGCAAGAATATACGCGGGGAAATTAAACGGGTATCACCTT‑GTGAAA‑GGTTGTTACGCCGGGGAGTAAAGC‑CCTTGTAACCGCCCCAACAGAGaat                                                                                                                           <  2:435360/147‑1 (MQ=255)                      aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGaatg                                                                                                                          >  3:140790/1‑150 (MQ=255)                      aCAAATGCAGCAAACTGGAAAAGTATCAAAAACAGCGTTGCAGAAAA‑CACAATCTGCCAGAATATCCGCGGGGAAATAAACCGGGTATCCCCAT‑GTGAAC‑CTTTGTTACGCGGGGGATTAAAGC‑CCCTGTAACCGCCCCAACAGAGaatg                                                                                                                          <  3:1416071/150‑1 (MQ=255)                        aaaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGaatga                                                                                                                         >  3:113043/1‑149 (MQ=255)                        aaaTGCAGCAAACTGAAAAAGTATCAAAAACAGCGTGGAAGGAAC‑CACAATCTGCCAGAATATACGCGGGGAAATAAAACGGGTATCACCAT‑GTAAAC‑CGTTGTTCCGCGGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGaatg                                                                                                                          <  4:408707/148‑1 (MQ=255)                        aaaTGCAGCAAACGGAAAAAGTATCAAAAACAGCTTTGCAGGAAA‑CACAACCTGCCAAAATATACGCGGGAAAATTAACCGGGTACCACCAT‑GTAAAA‑CGTTGTTACGCCGGGGAGAAAAGC‑CCCTGAAACCGCCCCAACAGAGaat                                                                                                                           <  4:282027/147‑1 (MQ=255)                        aaaTGCAACAAACGGAAAAAGTATAAGAAACAGGTTTGAAGAAAC‑CACAATCGGCCAGAATATACGCGGGGAAATTAACCGGGTATCCCCGT‑TTCAAA‑CGTTGTTCCGCGGGGGAAAAAAGC‑CCCGTAAACCGCCCCAACAGAGaatgat                                                                                                                        <  4:569414/150‑1 (MQ=255)                        aaaTCCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAA‑CACAATCTGCCAAAATAAACGCGGGAAAATAAACCGGGTATCACCAG‑GGAAAA‑CTTTGTTACGCCGGGGAGAAAAGC‑CCCTGAAACCGCCCCAACAGAGaatgata                                                                                                                       <  1:106776/151‑1 (MQ=255)                         aaTGCACAAAACTGGAAAAGAACAAGAAACAGCGTTCAAGAAAA‑CCCAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTACCACCAT‑GGAAAC‑CGTGGTTAGGAGGGGGAGTAAGGC‑CCCAGTAACCGCCCCAACAGAGaatgataa                                                                                                                      <  4:332065/151‑1 (MQ=255)                             cagcaaatgaaaaaaTTATCAAAAACAGTTTTGAAGAAAA‑CAAACTCTGCGAAAATGTCGGTTGGGAAATAAAACGGGTACCCCCAG‑TTAAAA‑CGTTGTTGCCCGGGGGAAAAAGGC‑CCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                                     <  4:22071/137‑1 (MQ=255)                                caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATgataatgataat                                                                                                               >  2:71172/1‑151 (MQ=255)                                caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGGATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                                     >  1:961007/1‑145 (MQ=255)                                 aaaCTGAAAAAGTATAAGAAACAGCGTTGCAGGAAA‑CACAATCTGCAAGAATATACGCGGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCGGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGAATgataatgataat                                                                                                               <  3:1306728/150‑1 (MQ=255)                                 aaaCGGGAAAAGTATCAGAAACAGCGTTGCAGGAAC‑CACAATCTGCCAAAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTAAAA‑CGTTGTTACGCGGGGGAGTAAAGC‑CCCTGAACCCGCCCCAACAGAGAAtgataatgataa                                                                                                                <  1:689807/149‑1 (MQ=255)                                 aaaCGGGAAAAGTATAAAAAACAGTGTTGCAGAAAA‑CAAAACCGGCCAAAATATCCGCTGGAAAATAAAACGGTTATCACCAT‑GTAAAA‑CGTTGTTACGCGGGGGAATAAAGC‑CCCTGTAACCGCCCCAACAGAGAatgataatgata                                                                                                                 <  1:486435/148‑1 (MQ=255)                                  aaCGGGAAAAGTATCAAAAACAGCGTTGCAGGAAA‑CACAATCTGCAAGAATATACGCTGGAAAATAAAACGGTTATCCCCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGAtaataa                                                                                                             <  3:928083/151‑1 (MQ=255)                                   aCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGAtaataa                                                                                                             >  3:366313/1‑150 (MQ=255)                                     tGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAt                                                                                                          >  3:1163801/1‑151 (MQ=255)                                         aaaGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                                                                                      >  1:912074/1‑151 (MQ=255)                                         aaaGTATCAAAAACAGCGTTGCAGAAAA‑CACAATCTGCCAAAATATACGCTGGGAAATAAAACGGGTATCACCAT‑GTAAAA‑CGTTTTTACGCGGGGGAAAAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGc                                                                                                      <  4:227975/151‑1 (MQ=255)                                              ataaaaaaCGGGGTGGCAGAAAA‑CAAAACCTGCAAAAATATGCGCTGGAAAATTAAAC‑GGCAACCCCAT‑TGGAAATTGTTGTTCGGCCGGGGAGAAAAGA‑CCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                                                                                                  <  1:215346/145‑1 (MQ=255)                                                   aaaCAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                                            >  1:1023686/1‑151 (MQ=255)                                                   aaaCAGCGTGGCAGAAAA‑CAAAATCGGCCAGAATATCCGCGGGAAAATAAAACGGTTACCCCCAT‑GTAAAA‑CGTTGTTCCGCCGGGAAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                                   <  1:857770/144‑1 (MQ=255)                                                    aaacaggtttgaGGAAACCACAATATTCAAGAATTCACGCTGGAAAATTAAAGGGTTAACCCAAT‑TGGAAG‑CTTTTTTCCCCGGGGGAGAAAAGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                                            <  3:1167152/140‑1 (MQ=255)                                                    aaCAGCGTTACAGGAAA‑CACATACGGCAAAAAAATACGC‑GGGAAAATAAACCGCAATACCCATGGTAAGA‑CTTTGTTACTCGGGGTAGAAAAGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATa                                                                                           <  1:1014642/151‑1 (MQ=255)                                                      cAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGc                                                                                         >  1:632218/1‑151 (MQ=255)                                                         cGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                                                                                      >  3:66689/1‑151 (MQ=255)                                                          gTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                     >  1:1411612/1‑151 (MQ=255)                                                           ttGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc                                                                                    >  3:620384/1‑151 (MQ=255)                                                           ttGCAGAAAA‑CACAATCTGCCAGAATATACGCTGGAAAATTAAACGGGTATCCCCAT‑GTAAAA‑CGTTTTTCCGCCGGGGATAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                                     <  3:1113058/150‑1 (MQ=255)                                                                gaaaa‑CACATCCTCCCAAAATATCCGCTGGAAAATAAACCGGTTACCCCCAT‑GTAAAA‑CGTGGTTACGCGGGGAAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                                <  1:300973/148‑1 (MQ=255)                                                                gaaaa‑CACAATCTGCCAGAATATACGCTGGAAAATTAAACGGGTATCACCAT‑GTAAAA‑CGTTGTTACGCGGGGGAGAAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa                                                                               <  3:1064744/149‑1 (MQ=255)                                                                  aaa‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa                                                                              >  3:1115836/1‑150 (MQ=255)                                                                  aaa‑CACAATCTGCCAAAATATACGCGGGAAAATTAAACGGTTATCCCCAT‑GTAAAA‑CGTTGTTACGCGGGGGAGTAAAGC‑CCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaaa                                                                             <  3:13804/151‑1 (MQ=255)                                                                    a‑cacAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc                                                                            >  1:31024/1‑150 (MQ=255)                                                                       acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                          >  1:259912/1‑150 (MQ=255)                                                                        caATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                          >  1:502431/1‑149 (MQ=255)                                                                           tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc                                                                     >  1:431940/1‑151 (MQ=255)                                                                               ccAGAATATACGCGGGAAAATTAAACGGGTCTCCCCAT‑GTAAAA‑CTTTGTTACCCGGGGGAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                                                   <  2:224993/149‑1 (MQ=255)                                                                               caaGAATATACGCTGGGAAATAAAACGGGTATCACCAT‑GTAAAC‑CGTTGTTCCGCGGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                                    <  3:87414/146‑1 (MQ=255)                                                                                 aaaaTATCCGCTGGGAAATTAAACGGGTACCCCTAT‑GAAAAA‑CGTTGTTACGCCGGGGAGTAAAGC‑CCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                                                   <  3:1268430/145‑1 (MQ=255)                                                                                 aGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                                 >  1:537425/1‑149 (MQ=255)                                                                                  gAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                              >  1:1344528/1‑151 (MQ=255)                                                                                   aaTATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                             >  3:68828/1‑151 (MQ=255)                                                                                   aaTATACGCGGGAAAATAAAACGGTTACCCCCAT‑GTGAAA‑CTTTGTTCCGCGGGGGAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctc                                                                     <  3:1395866/143‑1 (MQ=255)                                                                                           cTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGAc                                                      >  3:1023505/1‑150 (MQ=255)                                                                                           cTGGAAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCGGGGGAGAAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                         <  2:288119/147‑1 (MQ=255)                                                                                             gaaaaaTAAACCGGGTACCACCAG‑GTGAAA‑TGTTGTTCCGCGGGGAAGAAAACA‑CCCTTTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACa                                                     <  1:110686/146‑1 (MQ=255)                                                                                                 aaTTAAAAGGGGATCACCAT‑GGGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATc                                               >  1:396555/1‑151 (MQ=255)                                                                                                 aaTAAACCGGGTATCACCAT‑GTAAAA‑CGTTGTTACGCGGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATc                                               <  1:977933/151‑1 (MQ=255)                                                                                                   ttAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATccc                                             >  4:21718/1‑151 (MQ=255)                                                                                                    tAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                            >  3:714603/1‑151 (MQ=255)                                                                                                     aaaCGGTTATCACCAT‑GTGAAA‑CGTTGTTCCGCCGGGGAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                           <  3:745146/151‑1 (MQ=255)                                                                                                     aaaCGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                            >  1:467932/1‑150 (MQ=255)                                                                                                     aaaCGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                            >  1:960177/1‑150 (MQ=255)                                                                                                       aCGGATTTCGCCAT‑TAGAAA‑CGTTATTTCGCGGGGGATAAAACA‑CCTCATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTaa                                          <  3:1021631/150‑1 (MQ=255)                                                                                                        cGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                        >  3:1301037/1‑151 (MQ=255)                                                                                                         gggTATCCCCAT‑GTAAAA‑CTTTGTTCCGCCGGGGAGTAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                        <  3:226181/150‑1 (MQ=255)                                                                                                         gggTATCACCAT‑GAAAAA‑CGTTTTTACGCCGGGGAGAAAAGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                       <  1:433650/151‑1 (MQ=255)                                                                                                         gggTACCACCAT‑GTGAAA‑CGTTTTTACGCCGGGGAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                       <  3:764739/151‑1 (MQ=255)                                                                                                         gggTACCACCAT‑GAAAAA‑TTTTTTTACGCGGGGGAAAAAAGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                        <  4:370807/150‑1 (MQ=255)                                                                                                            tATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGAt                                    >  4:372841/1‑151 (MQ=255)                                                                                                               ccccAT‑GTAAAA‑CGTTGTTACGCCGGGGAGAAAAGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                                  <  3:749176/148‑1 (MQ=255)                                                                                                                    t‑gtgAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                             >  3:558603/1‑150 (MQ=255)                                                                                                                       tgAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAt                          >  1:302065/1‑151 (MQ=255)                                                                                                                                ttttACGCGGGGGATAAAGGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc                   <  3:1023814/148‑1 (MQ=255)                                                                                                                                   tACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCt               >  3:280516/1‑151 (MQ=255)                                                                                                                                   tACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCt               >  1:762651/1‑151 (MQ=255)                                                                                                                                     cGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTg              >  3:78225/1‑150 (MQ=255)                                                                                                                                      gCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCCGa             >  3:1183906/1‑150 (MQ=255)                                                                                                                                         ggggAGAAAAGC‑CCTTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGa             <  3:1058015/147‑1 (MQ=255)                                                                                                                                         ggggAGAAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaa          <  3:602780/150‑1 (MQ=255)                                                                                                                                          gggAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAc        >  4:73221/1‑151 (MQ=255)                                                                                                                                          gggAGTAAAACCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTg              <  1:1392232/146‑1 (MQ=255)                                                                                                                                             agaaaaGC‑CCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatc     <  1:581068/148‑1 (MQ=255)                                                                                                                                             aGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatc     >  1:28427/1‑151 (MQ=255)                                                                                                                                                aaaGC‑CCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatca    <  1:1115705/149‑1 (MQ=255)                                                                                                                                                aaaGC‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcat   <  3:862697/150‑1 (MQ=255)                                                                                                                                                aaaGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcatc  >  1:1257681/1‑151 (MQ=255)                                                                                                                                                aaaGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAAcatcata  <  1:366267/151‑2 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAA‑CACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAA‑CGTTGTTACGCCGGGGAGTAAAGA‑CCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAACATCATC  >  NC_000913/2628481‑2628771

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