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

CTCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAA  >  NC_000913/2628487‑2628761                                                                                                                                            |                                                                                                                                         cTCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                         >  2:2284242/1‑140 (MQ=255) cTCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGCATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                         >  2:2092850/1‑140 (MQ=255)  tCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                                                                                        <  1:74280/140‑1 (MQ=255)  tCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                                                                                        <  1:2243119/140‑1 (MQ=255)   cGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAccc                                                                                                                                       <  2:2058594/140‑1 (MQ=255)         tAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAAc                                                                                                                                 <  2:2235835/140‑1 (MQ=255)         aaaTCCTACAAATGCAGAAAAC‑TGAAAAAGTATAAAAAACAGCGTTGCAGGAAACACAATCTGCCAAAATATCCGCTGGAAAATTAACCGGGTATCACCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAAc                                                                                                                                 <  2:2052081/139‑1 (MQ=255)            cccTACAAATGCAGCAAAC‑TGGAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCTGGAAAATAAAACGGGTATCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGc                                                                                                                              <  2:915094/139‑1 (MQ=255)             ccTAAAAATGCAGCAAAC‑GGAAAAAGTATAAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGGAAATAAAGCGGGTATCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGcc                                                                                                                             <  2:228256/140‑1 (MQ=255)              caaaaaaTGCAGCAAAC‑TGAAAAAGTATCAGAAACAGCGTGGCAGGAAACACAATCTGCCAGAATATCCGCTGGGAAATTAACCGGGTACCCCCATGTGAAAGGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGccc                                                                                                                            <  2:816686/136‑1 (MQ=255)              cTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                            >  1:106188/1‑140 (MQ=255)                aCAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                          >  1:1516088/1‑140 (MQ=255)                aCAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                          >  1:211903/1‑140 (MQ=255)                  aaaTGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAAc                                                                                                                        >  1:2336981/1‑140 (MQ=255)                   aaTGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACa                                                                                                                       >  2:518854/1‑140 (MQ=255)                   aaTGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATAGGCTGGAAAATTAAACGGGTATCCCCATGGGAACCGTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACa                                                                                                                       <  2:1545096/140‑1 (MQ=255)                   aaTGCAGCAAAC‑TGGAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                       <  2:2155149/140‑1 (MQ=255)                   aaTGCAGAAAAC‑TGAAAAAGTATCAAAAACAGCGTGGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATAAACCGGGTATCACCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                       <  2:211903/140‑1 (MQ=255)                   aaTGAAGCAAAC‑TGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                       <  2:809798/140‑1 (MQ=255)                     atggacaAACTTGGAAAAGTATCAGAAACGGCGTTGCAGGAAAAACACACTGCCAAAAAATCCGCTGGAAAATTAAAGGGGTTCCACCTGGTGAAAGTTTGTTACGCCGGGGAGAAAAGACCCTGTAACCGCCCCAACag                                                                                                                      <  2:625705/135‑1 (MQ=255)                      gaagaaaaC‑GGAAAAAGTATAATAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATAAACCGGTTATCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                    <  2:1425932/135‑1 (MQ=255)                       cagcaAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                   >  1:1363750/1‑140 (MQ=255)                       cagcaAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                   >  2:1103177/1‑140 (MQ=255)                              c‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataa                                                                                                            >  1:783435/1‑140 (MQ=255)                               cggaaaaaGAATCAGAAGCAGCGTCGAAGGAAGCACAATCGCCAAGAAAATGCGGTGGTAATTTAAAGGGATACCACTAGATGAACCGTTGTTACCCGGGGGAGTAATGCCCTTGTAGCCGCCCCAACAGAGaatgataa                                                                                                            <  2:1516088/136‑1 (MQ=255)                                       gTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGataata                                                                                                    >  1:734932/1‑140 (MQ=255)                                            aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                               >  1:247942/1‑140 (MQ=255)                                              aaaCAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                             <  2:734932/140‑1 (MQ=255)                                                 cAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                          >  1:785271/1‑140 (MQ=255)                                                    cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                                       >  1:314551/1‑140 (MQ=255)                                                              aaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt                                                                             >  1:254314/1‑140 (MQ=255)                                                                aaaaaTCTGCAAGAATATACGAGGGGAAATAAAAGGGATACCCCCATGTAAACCTTTTTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                           <  2:1705039/137‑1 (MQ=255)                                                                   aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                        >  2:848003/1‑140 (MQ=255)                                                                     tCTGCCAGAATATCCGCGGGAAATTTAAACGGGTATCCCCATGTAAAAGGTTGTTCCGCGGGGGAAAAAACCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                      <  2:1097207/140‑1 (MQ=255)                                                                     tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                      >  1:1531458/1‑140 (MQ=255)                                                                     tCTGCAAAAATATACGCTGGAAAATAAAGCGGTTATCCCCATGTGAAATGTTTTTCCGCCGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                      <  2:1648248/140‑1 (MQ=255)                                                                           aGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                >  1:1273375/1‑140 (MQ=255)                                                                             aaTATCCGCGGGAAAATTAAACGGGTATCACCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                              <  2:144895/140‑1 (MQ=255)                                                                             aaTATCCGCGGGAAAATAAAAGGGTTACCCCCATGTAAAACTTTGTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                              <  2:2560047/140‑1 (MQ=255)                                                                             aaTATCCGCGGGAAAATAAAACGGGTACCACCATGTAAAAGGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                              <  2:696769/140‑1 (MQ=255)                                                                             aaTATACGCGGGAAAATAAAACGGGTATCCCCATGTAAAACTTTTTTCCGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                              <  2:1273375/140‑1 (MQ=255)                                                                               tatcCGCGGGGAAATTAAAGGGGTATCCCCATGTGAAACGTTTTTACGCGGGGGAATAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                            <  2:189039/140‑1 (MQ=255)                                                                               tataCGCTGGAAAATTAAACGGGTATCCCCATGTAAAACTTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                            <  1:1469394/140‑1 (MQ=255)                                                                                   cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                        >  1:2188004/1‑140 (MQ=255)                                                                                   cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                        >  2:1584143/1‑140 (MQ=255)                                                                                     cTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCat                                                      >  2:14357/1‑140 (MQ=255)                                                                                        ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                   >  1:189039/1‑140 (MQ=255)                                                                                         gAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAACTCTGCATAttt                                                  >  1:1528681/1‑140 (MQ=255)                                                                                           aaTTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                >  1:148682/1‑140 (MQ=255)                                                                                                  cGGGTACCACCATGTAAACCTTTTTTACGCGGGGAAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                         <  2:785271/140‑1 (MQ=255)                                                                                                             ttgtgAATCGTTTTTCCGCCGGGGAGTAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                              <  2:783435/139‑1 (MQ=255)                                                                                                             aTGTGAACCGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                              <  1:864522/140‑1 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                            >  2:188856/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                            >  2:1325122/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                            >  1:2084122/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                            >  2:2413083/1‑140 (MQ=255)                                                                                                               gtaaaccGTTTTTCCGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                            <  2:148682/134‑1 (MQ=255)                                                                                                                  aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                         >  2:168968/1‑140 (MQ=255)                                                                                                                   aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                        >  1:1428459/1‑140 (MQ=255)                                                                                                                      gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                     >  2:835503/1‑140 (MQ=255)                                                                                                                      gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                     >  2:949364/1‑140 (MQ=255)                                                                                                                        tgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                   >  2:283320/1‑140 (MQ=255)                                                                                                                        tgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGc                   >  1:150292/1‑140 (MQ=255)                                                                                                                         gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCa                  >  2:1475938/1‑140 (MQ=255)                                                                                                                         gttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCa                  >  1:561880/1‑140 (MQ=255)                                                                                                                          ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                 >  2:111518/1‑140 (MQ=255)                                                                                                                                  gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc         >  2:1470/1‑140 (MQ=255)                                                                                                                                  gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc         >  1:1159398/1‑140 (MQ=255)                                                                                                                                   ggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCa        >  1:1247988/1‑140 (MQ=255)                                                                                                                                     aGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGc      <  1:835503/140‑1 (MQ=255)                                                                                                                                         aaGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa  <  1:2284242/140‑1 (MQ=255)                                                                                                                                         aaGCCCCTATAGCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa  <  1:2092850/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                         CTCGCCAGTAATCCTACAAATGCAGCAAAC‑TGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAA  >  NC_000913/2628487‑2628761

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