Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,625	0	A	C	27.0%	157.5 / 28.6	74	V105G (GTC→GGC)	pdeF	cyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (47/7);  new base C (1/19);  total (48/26)
Fisher's exact test for biased strand distribution p-value = 5.40e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAA  >  NC_000913/2628490‑2628762                                                                                                                                         |                                                                                                                                           gCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                        <  2:2053603/140‑1 (MQ=255)    aGTAATCCTAAAAAGGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGAAAATAAACCGGGTATCCCCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCTTGTa                                                                                                                                     <  2:2339233/140‑1 (MQ=255)      tAATCCTACAAATGCAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑Ac                                                                                                                                  <  2:1985352/140‑1 (MQ=255)      tAATCCTAAAAATGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGAAAATTAAACGGTTACCCCCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑Ac                                                                                                                                  <  2:1920675/140‑1 (MQ=255)         tCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGc                                                                                                                               >  2:51693/1‑140 (MQ=255)             aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCa                                                                                                                           >  1:2563033/1‑140 (MQ=255)              aaaaTGCAGAAAACTGAAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCaa                                                                                                                          <  2:1671161/139‑1 (MQ=255)                aaggcagcaAACTGGAAAAGTATAAAAAACAGCGTTGCAGGAAACACAATCTGCAAAAATATCCGCGGGGAAATAAACCGGGTATCCCCATGTGAAAC‑GTTGTTACGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACa                                                                                                                        <  2:177205/137‑1 (MQ=255)                aaTGCAGCAAATTGGAAAAGTATCAAAAACAGCGTTGAAGGAAACACAGTCTGCAAGAATATACGCGGGGAAATTAAAGGGGTATCCCCATGTAAAAG‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACa                                                                                                                        <  2:1489698/140‑1 (MQ=255)                aaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACa                                                                                                                        >  2:2834106/1‑140 (MQ=255)                aaTGCAGCAAACTGAAAAATTATCAAAAACAGCGTTGCAGGAAACAAAATCTGCCAGAATATACGCGGGAAAATTAAAGGGGTACCCCCATGTAAAAC‑TTTGTTACGCGGGGGAGAAAAGCCCTTGTA‑ACCGCCCCAACa                                                                                                                        <  2:2372377/140‑1 (MQ=255)                aaTCCAGAAAACTGAAAAAGTACCAAAAACAGGGTGGCAGAAAACAAAATCGGCAAGAATATCCGGGGGAAATTTAAGCGGGAACCCCCATGTGAAAC‑GTTTTTCCGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACa                                                                                                                        <  2:775308/140‑1 (MQ=255)                  tGCAGCAAACGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAAAATATACGCGGGGAAATTAAACGGGTATCACCATGTGAAAC‑TTTGTTCCGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACaga                                                                                                                      <  2:2617822/140‑1 (MQ=255)                  tGCAACAAACTGGAAAAGTATCAAAAACAGGTTGGCAGAAAACACAATCTGCCAGAATATACGCGGGGAAATTAAAGGGGTACCACCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACaga                                                                                                                      <  2:177540/140‑1 (MQ=255)                   gcagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACagag                                                                                                                     >  1:2372377/1‑140 (MQ=255)                   gcagcaAACTGAAAAATTATAAAAAACAGCGTGGAAGAAAACAAAATCTGCAAGAATATACGCGGGGAAATAAAGCGGGTATCCCCTTGGAAAAC‑GTTGTTCCGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACagag                                                                                                                     <  2:2199114/140‑1 (MQ=255)                   gcagcaAACGGGAAAATTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTACCACCATGTGAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACagag                                                                                                                     <  2:223416/140‑1 (MQ=255)                    cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAgaga                                                                                                                    >  1:775308/1‑140 (MQ=255)                    cagaaAACTGAAAAAGTATCAAAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATCCGCGGGAAAATAAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAgaga                                                                                                                    <  2:1466268/140‑1 (MQ=255)                    caacaAACTGAAAAATTATAAGAAACGGCGTTGAAGGAAACACAATCGGCAAGAATATCCGGTGGGAAATTAAGGGGGTCTCACCTTGAGAAAT‑GTTGTTCCGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAgaga                                                                                                                    <  2:2545147/137‑1 (MQ=255)                     agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGaa                                                                                                                   >  2:577758/1‑140 (MQ=255)                     agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGaa                                                                                                                   >  2:1905957/1‑140 (MQ=255)                        aaaCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGaatga                                                                                                                >  1:1927223/1‑140 (MQ=255)                         aaCTGAAAAAGTAAAAGAAACAGCGTTGCAGGAAACACAACCTGCCAAAATATACGCGGGAAAATAAAGGGGGTACCCCCATGTGAAAC‑TTTGTTCCGCGGGGGAGAAAAGCCCCTGAA‑ACCGCCCCAACAGAGaatgat                                                                                                               <  2:2563033/140‑1 (MQ=255)                          aCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGaatgata                                                                                                              >  1:2367441/1‑140 (MQ=255)                            tGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGaatgataat                                                                                                            >  2:2191973/1‑140 (MQ=255)                                 aaGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCTGGGAAATTAAAGGGGTATCACCATGTGAAAC‑GTTGTTACGCGGGGGAGTAAAGCCCCTGAA‑ACCGCCCCAACAGAGAAtgataatgataa                                                                                                       <  2:2325885/140‑1 (MQ=255)                                   gTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGataata                                                                                                     >  1:1158412/1‑140 (MQ=255)                                      tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                  >  1:2822078/1‑140 (MQ=255)                                      tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGCAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                  >  1:239579/1‑140 (MQ=255)                                          aaagatcGTTGCAGAAAACACAATCTGCAAGAATCGCCGCTGGGAAATTAAACGGCTAAC‑CCAGTTGGAACTTTTGTAACGCCGGGGAGAAAAGACCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                              <  2:1592596/134‑1 (MQ=255)                                          aaaaggCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGAAAATAAAGCGGGTATCCCCATGTGAAAC‑TTTGTTACGCGGGGGAGTAAAGCCCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                              <  2:1927223/135‑1 (MQ=255)                                          aaaCAGCGTTGCAGGAAACACATTCTGCAAGAAGATACGCTGGGAAATAAAACGGGTATCACCATGTAAATC‑GTTGTACCGCGGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                              <  2:2367441/140‑1 (MQ=255)                                          aaaCAGCGTTGCAGAAAACAAAATCTGCCAAAATATACGCGGGAAAATAAAAGGGGTACCCCCATGTAAAAC‑TTTTTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                              <  2:358690/140‑1 (MQ=255)                                          aaaCAGCGTGGCGGGAAACAAAATCTGCAAAAATATGCGCGGGAAAATTAAACGAGTACCCCCATGTAAAGC‑TTTTTAACGCCGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                              <  2:239579/140‑1 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                          >  1:358690/1‑140 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                          >  2:1525250/1‑140 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                          >  1:1123692/1‑140 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                          >  1:2352609/1‑140 (MQ=255)                                               gCGTTGCAGAAAACACAATCTCCAAAAATACCCGCGGGAAAATTAAACGGGTACCCCCATGTAAAAC‑GTTTTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                                                                                         <  2:1910763/140‑1 (MQ=255)                                               gCGTTGAAGAAAACACAACCGCCCAAAATATCCGCGGGAAAATAAAGCGGGTACCACCATGGGAAAG‑TTTTTTAGCCGGGGGAAAAAACACCC‑GTATACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                                                                                         <  2:1811983/140‑1 (MQ=255)                                               gCGTTGAAGAAAACAAAATCTGCAAGAATATCCGGGGGAAAATAAAAGGGGTATCCCCATGTAAACC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                                                                                         <  2:1158412/140‑1 (MQ=255)                                                cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                                        >  2:2427035/1‑140 (MQ=255)                                                          aaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt                                                                              >  1:1910763/1‑140 (MQ=255)                                                                 tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                       >  2:1921524/1‑140 (MQ=255)                                                                      cAGAATATCCGCGGGGAAATTAACCGGGTATCCCCATGTAAAAC‑GTTGTTACGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                  <  1:326112/140‑1 (MQ=255)                                                                      aaGAATATCCGCGGGAAAATAAAACGGGTATCACCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                  <  2:2624689/139‑1 (MQ=255)                                                                         aaTATGCGCTGGAAAATTAACCGGTTACCCCCATGTAAAAC‑GTTGTTACGCGGGGGAGAAAAGCCCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                               <  1:2191973/140‑1 (MQ=255)                                                                         aaTATACGCTGGAAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                               <  2:1025202/140‑1 (MQ=255)                                                                           tataCGCTGGAAAATTAAACGGGTATCCCCATGTAAAAC‑GTTGTTACGCCGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                             <  2:1123692/140‑1 (MQ=255)                                                                             taCGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                           >  2:1878976/1‑140 (MQ=255)                                                                             taCGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                           >  1:2337879/1‑140 (MQ=255)                                                                             taCGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATACCGTTTCGAAAAAACCCACctct                                                           >  1:974961/1‑140 (MQ=255)                                                                               cGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                         >  1:2624689/1‑140 (MQ=255)                                                                                   gggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                     >  1:1833384/1‑140 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAACGTTTCGTAAAAACCCAGCTCTGCATAtt                                                    >  1:1535452/1‑140 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                    >  1:114613/1‑140 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                    >  2:1647018/1‑140 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                    >  2:335479/1‑140 (MQ=255)                                                                                       aaTTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                 >  2:1152457/1‑140 (MQ=255)                                                                                          aaaaCGGGTATCCCCATGTAAAAC‑GTTTTTACGCGGGGGAGAAAAGCCCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                                              <  2:2352609/139‑1 (MQ=255)                                                                                             aCGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaa                                           <  1:1921524/140‑1 (MQ=255)                                                                                              cGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                          <  1:577758/140‑1 (MQ=255)                                                                                              cGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                          >  2:133078/1‑140 (MQ=255)                                                                                                   acccccATGTAAAAC‑GTTTTTCCGCGGGGGAGAAAAGCCCTTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                     <  2:2751569/136‑1 (MQ=255)                                                                                                    ccACCATGTAAAAC‑GTTTTTACGCGGGGGAAAAAAGCCCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATccc                                    <  2:2822078/139‑1 (MQ=255)                                                                                                     ccccATGTAAAAC‑TTTGTTACGCGGGGGAGAAAAGCCCCTGAA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                   <  2:766018/138‑1 (MQ=255)                                                                                                     cACCATGAAAAAT‑GTTATTAGGCCGGGGAGAAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                   <  2:1535452/140‑1 (MQ=255)                                                                                                         aTGTGAAAT‑TTTTTTACGCGGGGGAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                               <  2:1833384/140‑1 (MQ=255)                                                                                                         aTGTAAAAG‑GTTGTTACGCGGGGGAATAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                               <  2:1270065/140‑1 (MQ=255)                                                                                                          tgtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                              >  2:965322/1‑140 (MQ=255)                                                                                                           gtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                             >  1:1507142/1‑140 (MQ=255)                                                                                                           gtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                             >  1:644176/1‑140 (MQ=255)                                                                                                           gtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                             >  1:1518620/1‑140 (MQ=255)                                                                                                           gtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                             >  1:181696/1‑140 (MQ=255)                                                                                                           gtgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                             >  2:2665243/1‑140 (MQ=255)                                                                                                            tgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                              <  1:1368746/138‑1 (MQ=255)                                                                                                            tgAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                              >  2:1368746/1‑138 (MQ=255)                                                                                                                aC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATccc                        >  2:1456575/1‑140 (MQ=255)                                                                                                                    ttgttACGCGGGGGAGAAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                     <  2:114613/140‑1 (MQ=255)                                                                                                                     tgttACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                    >  2:1893957/1‑140 (MQ=255)                                                                                                                     tgttACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                    >  1:1270065/1‑140 (MQ=255)                                                                                                                       ttACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                  >  2:1914974/1‑140 (MQ=255)                                                                                                                            ccGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccg             >  2:594463/1‑140 (MQ=255)                                                                                                                             cGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgc            >  1:2380403/1‑140 (MQ=255)                                                                                                                              ggggAGTAAAGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc           <  1:1456575/140‑1 (MQ=255)                                                                                                                               gggAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc          <  1:1647018/140‑1 (MQ=255)                                                                                                                               gggAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc          >  1:144319/1‑140 (MQ=255)                                                                                                                               ggaataaaaGCCCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc          <  2:1507142/134‑1 (MQ=255)                                                                                                                                    tAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTg     >  1:2094171/1‑140 (MQ=255)                                                                                                                                     aaaGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGa    <  1:1525250/140‑1 (MQ=255)                                                                                                                                      aaGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa   <  1:2427035/140‑1 (MQ=255)                                                                                                                                       aGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaa  >  2:1337308/1‑140 (MQ=255)                                                                                                                                       aGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaa  >  2:95379/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                           GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAAC‑GTTGTTACGCCGGGGAGTAAAGACCCTGTA‑ACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAA  >  NC_000913/2628490‑2628762

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