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

GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628490‑2628764                                                                                                                                        |                                                                                                                                            gCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                         >  1:2142440/1‑140 (MQ=255)   cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGt                                                                                                                                       >  1:775984/1‑140 (MQ=255)    aGTAATCCTAAAAATGCAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTa                                                                                                                                      <  2:2794796/140‑1 (MQ=255)      tAACCCTAAAAATGCAGCAAACTGGAAAAGTATCAAAAACGGCGTGGCAGAAAACACAATCTGCAAGAATATCCGCTGGAAAATAAAACGGGTACCCCCATGTGAAACGTTGTTACGCGGGGGAAAAAAGCCCCTGTAAc                                                                                                                                    <  2:2702997/140‑1 (MQ=255)        aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCg                                                                                                                                  >  1:448500/1‑140 (MQ=255)        aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCg                                                                                                                                  >  1:1576656/1‑140 (MQ=255)             aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  1:2196869/1‑140 (MQ=255)              aaaaTGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACAAAATCTGCCAGAATATACGCGGGGAAATAAAACGGGTATCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCaa                                                                                                                            <  2:2231075/139‑1 (MQ=255)                aaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACa                                                                                                                          >  1:2702997/1‑140 (MQ=255)                aaTGCAACAAACTGAAAAAGTAAAAGAAACAGCGTTGAAGAAAACAAAATCTGCCAGAATATACGTTGGAAAATAAAACGGTTACCACCATGTAAAACGTTGTTACGCGGGGGAGTAAAGCCCCTTTAACCGCCCCAACa                                                                                                                          <  2:2221989/140‑1 (MQ=255)                aaTCCAGCAAACGGAAAAAGTATAAAAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:2196869/140‑1 (MQ=255)                 atgaaaaaaaCTGAAAAAGTATAAAAAACAGCGTTGCAGAAAACAAATTCTGCCAAAATATACGCTGGGAAATTAACCGGGTATCACCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACag                                                                                                                         <  2:3065398/133‑1 (MQ=255)                   gcagcaAACTGAAAAAGTATAAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCGGGAAAATAAAAGGGGAATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                       <  2:306158/140‑1 (MQ=255)                    cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                      >  1:2648579/1‑140 (MQ=255)                    cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAgaga                                                                                                                      >  1:2642264/1‑140 (MQ=255)                    cagaaAACTGGAAAAGTATAAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGGAAATAAACCGGGTATCCCCATGTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGAAACCGCCCCAACAgaga                                                                                                                      <  2:1277043/140‑1 (MQ=255)                     agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaa                                                                                                                     >  1:1997557/1‑140 (MQ=255)                     agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaa                                                                                                                     >  1:1277043/1‑140 (MQ=255)                     agcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaa                                                                                                                     >  1:203069/1‑140 (MQ=255)                      gaaaaCGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAAAATATCCGCTGGGAAATTAAACGGGTATCACCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaat                                                                                                                    <  2:1948000/138‑1 (MQ=255)                            ggaaaaaGAATAAAAAACGGCGTTACAGGAAACACAATCTGCCAGAATATACCGTGGAAAATTAAACGGATATCACCATGTGAACCGTTGTTGCGCGGGGGAGAAAAGACCGTGTAACCGCCCCAACAGAGaatgataat                                                                                                              <  2:1281415/137‑1 (MQ=255)                              gAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatgataatga                                                                                                            >  1:1281415/1‑140 (MQ=255)                                      tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:1752860/1‑140 (MQ=255)                                      tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:1991019/1‑140 (MQ=255)                                      tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:334917/1‑140 (MQ=255)                                        aaaaaCAGCGTTGCAGAAAACACAATCTGCCAAAATATACGCTGGAAAATAAAACGGGTACCACCATGTAAAACGTTGTTACGCGGGGGATTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                                  <  2:2648579/138‑1 (MQ=255)                                          aaaCAGTTTTGGAGAAAACACAATCTGCCAGAATATCCGCGGGGAAATAAAGCGGGTATCGCCAGATGAAACGTTGTTACGCGGGGGATTAAGGGCCTTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:334917/140‑1 (MQ=255)                                          aaaCAGCGTTGAAGAAAACACAATCTGCAAAAATATACGCGGGAAAATTAAACGGGTACCACCATGTAAAACTTTGTTCCGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:1997557/140‑1 (MQ=255)                                          aaaCAGCGTCGCAGGAAACACACTCTGCGAGAATATAGGCGGGAAAATTAACCGGGTACCACCTTGTGAAGGTTTGTTACCCCGGGGAGGAAAGCCCCTGTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:2775230/140‑1 (MQ=255)                                            acggCGTTGCAGAAAACACAATCTGCAAAAATATACGCGGGAAAATAAAGGGGGTACCCCCATGAAAAACGTTGTTCCGCGGGGGAGAAAACCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                              <  2:786116/137‑1 (MQ=255)                                            aCAGCGTTGCAGAAAACAAAATCTGCCAGAATATACGCGGGAAAATTAAACGGGTATCCCCATGTAAAACTTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                              <  2:1701387/140‑1 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:626029/1‑140 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:786116/1‑140 (MQ=255)                                              aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:1420272/1‑140 (MQ=255)                                              aGCGTTGCAGAAAACACAATCTGCAAAAATATACGCGGGAAAATAAAACGGGTATCACCATGTGAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            <  2:626029/140‑1 (MQ=255)                                                 gTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTaa                                                                                         >  1:2775230/1‑140 (MQ=255)                                                               aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                           >  2:2171104/1‑140 (MQ=255)                                                               aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAAACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                           >  1:1637071/1‑140 (MQ=255)                                                                    gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaaa                                                                      >  1:2914675/1‑140 (MQ=255)                                                                      cAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    >  2:1230937/1‑140 (MQ=255)                                                                       aGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                   >  1:2013382/1‑140 (MQ=255)                                                                         aaaatccgcgGGAAAATAAAACGGTAACCCCAAAGTCAACCGTTGTTACGCCGGGGAGAAAAGCCCCTAAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  2:2790728/134‑1 (MQ=255)                                                                         aaTATCCGGGGGAAAATTAAACGGGTATCACCATGTAAAAGGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  2:856796/140‑1 (MQ=255)                                                                         aaTATACGCTGGAAAATAAAAGGGGTACCCCCATGTAAAACGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  1:2789095/140‑1 (MQ=255)                                                                         aaTATACGCGGGAAAATAAACCGGGTATCACCATGTAAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  2:2013382/140‑1 (MQ=255)                                                                           tataCGCGGGGAAATTAACCGGGTACCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAACCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                                               <  2:1637071/140‑1 (MQ=255)                                                                             taCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                             >  1:2107178/1‑140 (MQ=255)                                                                             taCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                             >  1:1444992/1‑140 (MQ=255)                                                                                gcgggaaaaTTAAACGGGTATCACCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                          <  2:2107178/134‑1 (MQ=255)                                                                                gcgggaaaaTTAAACGGGTACCCCCATGTAAAACGTTTTTACGCGGGGGAGTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                          <  2:1444992/134‑1 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:2741527/1‑140 (MQ=255)                                                                                    ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:856796/1‑140 (MQ=255)                                                                                           aaaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGAc                                               >  2:1364486/1‑140 (MQ=255)                                                                                              gggttaccACCATGTAAAACTTTGTTACGCGGGGGAGAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  2:1420272/133‑1 (MQ=255)                                                                                              cGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  1:1230937/140‑1 (MQ=255)                                                                                                   aTCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTACGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                       >  1:1917795/1‑140 (MQ=255)                                                                                                         aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 <  1:1364486/140‑1 (MQ=255)                                                                                                         aTGTGAAACGTTGTTACGCCGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 <  1:2171104/140‑1 (MQ=255)                                                                                                         aTGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 <  2:963970/140‑1 (MQ=255)                                                                                                          tgtgAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  1:1972056/128‑1 (MQ=255)                                                                                                          tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            >  2:1972056/1‑128 (MQ=255)                                                                                                          tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                >  1:2668193/1‑140 (MQ=255)                                                                                                           gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCTCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:1763968/1‑140 (MQ=255)                                                                                                           gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCTCTAAGGCg                               >  2:518840/1‑140 (MQ=255)                                                                                                              aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            >  1:963970/1‑140 (MQ=255)                                                                                                               aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                           >  2:173252/1‑140 (MQ=255)                                                                                                                      ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    >  1:1260660/1‑140 (MQ=255)                                                                                                                        aCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATg                  <  2:2741527/140‑1 (MQ=255)                                                                                                                            cGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgc              >  1:954314/1‑140 (MQ=255)                                                                                                                              gggAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  2:954314/140‑1 (MQ=255)                                                                                                                              gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAGTGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGCCAATGCCGccc            >  2:1392620/1‑140 (MQ=255)                                                                                                                              gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  2:2015384/1‑140 (MQ=255)                                                                                                                              gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  2:1213493/1‑140 (MQ=255)                                                                                                                              gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTACGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  2:1917795/140‑1 (MQ=255)                                                                                                                                gAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg          <  1:2015384/140‑1 (MQ=255)                                                                                                                                        aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  2:3166103/130‑1 (MQ=255)                                                                                                                                        aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  1:3166103/1‑130 (MQ=255)                                                                                                                                        aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  <  1:173252/140‑1 (MQ=255)                                                                                                                                        aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa     >  1:843156/1‑137 (MQ=255)                                                                                                                                        aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa     <  2:843156/137‑1 (MQ=255)                                                                                                                                        |                                                                                                                                            GCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628490‑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: ‑