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

CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764                                                                                                                                            |                                                                                                                                            cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                            >  1:992116/1‑140 (MQ=255) cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                            >  1:2617904/1‑140 (MQ=255)         tAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCCCCATGTAAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAAc                                                                                                                                    <  2:2390238/140‑1 (MQ=255)         aaacccTAAAACTGCACAAAACTGGAAAAGTATCAGAAACAGCGTGGCAGGAAACACAATCTGCCAAAATATACGCTGGAAAATTAAAC‑GGGTTTCACCATGTGAACCGTTGTTCCGCGGGGGAGAAAAGACCCTTTAAc                                                                                                                                    <  2:893409/136‑1 (MQ=255)                aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGACACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  1:2406816/1‑140 (MQ=255)                aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  1:905087/1‑140 (MQ=255)                 aaaaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTGGCGGAAAACACAATCGGCCAAAATATACGCTGGAAAATAAAAG‑GGGTATCCCGATGTGAAAGGTTGTTACGCGGGGGAGAAAACCCCCTGAAACCGCCCCaa                                                                                                                            <  2:208266/139‑1 (MQ=255)                   aaggcagcaAACGGAAAAAGTATCAAAAACAGCGTTGAAGAAAACAAAATCTGCCAAAATATACGCGGGAAAATTAAAC‑GGGTATCCCCATGTGAAACTTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:1185282/137‑1 (MQ=255)                   aaTGCAGCAAATGGGAAAAGTATCAAAAACGGCGTTGCAGAAAACACAATCTGCAAGAATATACGCTGGAAAATTAAAC‑GGGTATCACCATGTAAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:905087/140‑1 (MQ=255)                       cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                      >  1:2418634/1‑140 (MQ=255)                           aaaCTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatga                                                                                                                  >  2:679909/1‑140 (MQ=255)                               tGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                              >  1:990681/1‑140 (MQ=255)                               tGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                              >  1:1515912/1‑140 (MQ=255)                                 gAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataactg                                                                                                            >  1:1185282/1‑137 (MQ=255)                                         acaaaaaCAGCGTTGAAGAAAACACAATCTGCAAGAATATCCGCGGGGAAATAAAAC‑GGGTATCCCCATGTAAACCGTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    <  2:1650200/136‑1 (MQ=255)                                          aaaaaaCGGCGTGGAAGAAAACAAAACCTGCAAGAATATACGCGGGAAAATAAACC‑GGTTACCACCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAt                                                                                                   <  2:2608960/137‑1 (MQ=255)                                           aaaaaCAGCGTTGCAGAAAACACAATCTGCCAAAATATACGCTGGAAAATTAAAG‑GGGTATCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                                  <  2:2575345/138‑1 (MQ=255)                                           aGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                                  >  2:1092958/1‑140 (MQ=255)                                                 aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:2575345/1‑140 (MQ=255)                                                 aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:576073/1‑140 (MQ=255)                                                        cAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATGATGATAATAAAATTAGCCGTTAATGAt                                                                                     >  1:1141722/1‑140 (MQ=255)                                                        cAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                                     >  2:257935/1‑140 (MQ=255)                                                              acacAATCTGCAAAAATATGCGCTGGAAATTTAAAC‑GGGTATCACCATGTAAAACTTTTTTCCGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                                                                               <  2:1409685/140‑1 (MQ=255)                                                                  aaTCTGCCAAAATATACGCGGGAAAATAAAGC‑GGGTACCACCATGTGAAACGTTTTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                           <  2:483040/140‑1 (MQ=255)                                                                   aTCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa                                                                          >  1:1703656/1‑140 (MQ=255)                                                                    tCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                         >  1:483040/1‑140 (MQ=255)                                                                     cTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa                                                                        >  1:2174023/1‑140 (MQ=255)                                                                         aaaaaTATCCGC‑GGGAAAATAAACGGGGTACCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATGATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    <  2:1141722/137‑1 (MQ=255)                                                                              tataCGCTGGAAAATTAACC‑GGGTATCACCATGTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                               <  1:679909/140‑1 (MQ=255)                                                                              tataCGCTGGAAAATAAAAC‑GGTTCTCACCGGGTGAAACGTTTTTCCGCGGGGGAGTAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                               <  2:1515912/140‑1 (MQ=255)                                                                                      gggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:1310532/1‑140 (MQ=255)                                                                                       ggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  1:84089/129‑1 (MQ=255)                                                                                       ggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 >  2:84089/1‑129 (MQ=255)                                                                                       ggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:1463020/1‑140 (MQ=255)                                                                                       ggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:1578752/1‑140 (MQ=255)                                                                                           aTTAAAC‑GGGTATCACCAGGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                  <  1:1946168/140‑1 (MQ=255)                                                                                                 c‑GGGTATCCCCATGTGAAACGTTGTTCGGCGGGGGAGTAAAGTCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  2:1703656/140‑1 (MQ=255)                                                                                                 c‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            >  1:1080039/1‑140 (MQ=255)                                                                                                  ggggTATCCCCATGTGAAAGGTTTTTACGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  2:576073/139‑1 (MQ=255)                                                                                                  ggggAATCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  2:986788/139‑1 (MQ=255)                                                                                                     gTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                                         >  1:986788/1‑140 (MQ=255)                                                                                                       aTCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                       >  1:2520210/1‑140 (MQ=255)                                                                                                         ccccttgtaAAACGTTTTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                     <  2:369031/135‑1 (MQ=255)                                                                                                         ccccATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCAAATTTGCGGACAAAAATCCCt                                     <  2:990681/138‑1 (MQ=255)                                                                                                         agccATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                     <  2:1946666/138‑1 (MQ=255)                                                                                                             aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 <  1:949387/140‑1 (MQ=255)                                                                                                              tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                >  2:272948/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGATTGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:967968/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  1:85548/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:1817870/1‑140 (MQ=255)                                                                                                                tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa                              >  1:2606/1‑140 (MQ=255)                                                                                                                  aaaCTTTGTTCCGCGGGGGAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            <  2:1080039/140‑1 (MQ=255)                                                                                                                  aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            >  1:1872106/1‑140 (MQ=255)                                                                                                                   aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                           >  1:369031/1‑140 (MQ=255)                                                                                                                     cGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCt                         >  1:408782/1‑140 (MQ=255)                                                                                                                      gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcccct                        >  1:1946666/1‑138 (MQ=255)                                                                                                                          ttccgcgGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    <  2:2418634/137‑1 (MQ=255)                                                                                                                          ttACGCGGGGGAAAAAAGACCCTATAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    <  2:408782/140‑1 (MQ=255)                                                                                                                            aCGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATg                  >  2:1736381/1‑140 (MQ=255)                                                                                                                                cGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAGACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgc              >  1:1063406/1‑140 (MQ=255)                                                                                                                                 ggggAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc             <  2:85548/140‑1 (MQ=255)                                                                                                                                 ggggAGAAAACCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc             <  2:2606/140‑1 (MQ=255)                                                                                                                                  gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  2:2461658/1‑140 (MQ=255)                                                                                                                                  gggAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  1:2461658/140‑1 (MQ=255)                                                                                                                                   ggaaaaaaGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCa           <  2:1872106/135‑1 (MQ=255)                                                                                                                                    gaaaaaaGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg          <  2:1310532/136‑1 (MQ=255)                                                                                                                                            aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  2:1449350/1‑140 (MQ=255)                                                                                                                                            aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  2:1779206/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764

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