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

GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628486‑2628764                                                                                                                                             |                                                                                                                                            gCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGa                                                                                                                                             <  1:1297611/140‑1 (MQ=255)  cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                            >  1:1052697/1‑140 (MQ=255)          tAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGTGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGGAAATTAAAC‑GGGTATCCCCATGTGAAACGTTGTTCCGCCGGGGAGTAAAGCCCCTGTAAc                                                                                                                                    <  2:857826/140‑1 (MQ=255)          tAACCCTAAAAATGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCGGGGAAATTAACC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAAc                                                                                                                                    <  2:69974/140‑1 (MQ=255)          aaaTCCTAAAAATGCAGCAAACGGGAAAAGTATCAGAAACAGCGTTGCAGAAAACAAAATCTGCCAGAATATACGCGGGGAAATAAAAG‑GGTTATCCCAATGTGAAACTTTTTCACGCGGGGGAGAAAAGCCCTTGTagc                                                                                                                                    <  2:1295261/139‑3 (MQ=255)               cTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                               >  1:1163186/1‑140 (MQ=255)               cTAAAAATGCAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATAAAAC‑GGGTACCACCATGTAAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGccc                                                                                                                               <  2:1501894/140‑1 (MQ=255)                   aaaTGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAAc                                                                                                                           >  1:1501894/1‑140 (MQ=255)                    aaTGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGAAGGAAACACAATCTGCAAGAATATACGCGGGAAAATTAACC‑GGGTATCCCCAGGTAAAAGTTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:319538/140‑1 (MQ=255)                      tGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACaga                                                                                                                        <  1:454079/140‑1 (MQ=255)                      tGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCGGGAAAATAAAAC‑GGGTATCCCCATGTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACaga                                                                                                                        <  1:1427208/140‑1 (MQ=255)                      ggcagaaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGGAAATTAAAC‑GGTTATCCCCATGTAAAACGTTTTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACaga                                                                                                                        <  2:526233/139‑1 (MQ=255)                       gcagcaAACTGAAAAAGTATCAGAAACAGCGTGGAAGGAAACACAATCTGCCAAAATATCCGCTGGGAAATTAACC‑GGTTACCACCATGTGAAACGTTGTTGCGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                       <  2:1194921/140‑1 (MQ=255)                       gaagaaaaCAGGAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAAAATATACGCTGGGAAATTAAAC‑GGGTATCCCCATGTGAAGCGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                       <  2:1148277/135‑1 (MQ=255)                        cagcaAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGC‑GGGAAAATAAACGGGGTATCACCATGTGAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAgaga                                                                                                                      <  2:166365/140‑1 (MQ=255)                        aagcaAACTGAAAAAGTATCAGAAACAGCGTGGCAGAAACCCCAATCTGCCAGAATAAACGCGGGGAAATTAAAC‑GGGTATCACCATGTAAAGCTTTGTTCCGCCGGGAAAAAAAGCCCCTGTAACCGCCCCAACAgaga                                                                                                                      <  2:1059350/139‑1 (MQ=255)                          gcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaat                                                                                                                    >  1:331874/1‑140 (MQ=255)                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatg                                                                                                                   >  1:166365/1‑140 (MQ=255)                           caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaatg                                                                                                                   >  1:1148277/1‑140 (MQ=255)                                ggaaaaaGTATCAAAAACGGCGTTGCAGAAAACACAATCTGCCAAAATATCCGCGGGAAAATAAAAC‑GGGTATCCCAATGTAAAACGTTGTTACGCGGGGGAAAAAAGCCCCTGAAACCGCCCCAACAGAGaatgataat                                                                                                              <  2:742247/137‑1 (MQ=255)                                ggaaaaaGTATCAAAAACAGCGTGGCAGAAAACACAATCTGCCAGAATATACGCTGGAAAATAAAAC‑GGGTATCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                              <  2:553328/137‑1 (MQ=255)                                   aaaaTTATAAAAAACGGCGTGGCAGAAAACACAATCTGCAAGAATATACGCGGGAAAATAAACC‑GGTTATCCCCATGTAAAACGTTGTTACGCCGGGGAATAAACCCCCTGAAACCGCCCCAACAGAGaatgataatgat                                                                                                           <  2:791266/140‑1 (MQ=255)                                    aaaGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCGGGAAAATTAAAG‑GGGTATCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAatgataatgata                                                                                                          <  2:986404/140‑1 (MQ=255)                                     aaGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAAtgataatgataa                                                                                                         >  1:1409228/1‑140 (MQ=255)                                          tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:770151/1‑140 (MQ=255)                                          tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:816588/1‑140 (MQ=255)                                          tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  1:64756/1‑140 (MQ=255)                                            aaaaaCAGCGTTGCAGAAAACACAATCTGCCAAAATATGCGGGGGAAAATTAAAC‑GGGTATCCCCATGTGAACCGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                                  <  2:1083003/138‑1 (MQ=255)                                                aCAGCGTGGCAGAAAACACAATCTGCAAGAATATACGGTGGAAAATAAAAC‑GGGTATCCCCATGTAAAACGTTTTTCCGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                              <  2:475175/140‑1 (MQ=255)                                                  aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:475175/1‑140 (MQ=255)                                                        gCAGGAAACACAATCTGACAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAA‑GATAATAAAATTAGCCGTTAATGAt                                                                                     >  2:23758/1‑140 (MQ=255)                                                              aaCACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt                                                                                >  2:884380/1‑140 (MQ=255)                                                                aaaaaTCTGCGAGAATATCCGCTGGGAATTAAAGC‑GGTTAACCCCATGTAAAACTTTGTTACGCGGGAGAGAAAAAACCCGGTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                                              <  2:331874/137‑1 (MQ=255)                                                                   aaTCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                                           >  2:1222730/1‑140 (MQ=255)                                                                        gcaaGAATATCCGCGGGAAAATAAAAG‑GGGTATCACCTTGTAAAACTTTTTTCCGCGGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaaa                                                                      <  2:1052697/137‑1 (MQ=255)                                                                             aaTATACGCGGGGAAATAAAAC‑GGGTATCACCATGTAAAACGTTGTTCCGCCGGGGAAAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 <  2:1161891/140‑1 (MQ=255)                                                                                   cGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                           >  2:1313210/1‑140 (MQ=255)                                                                                    gCTGGAAAATTAAAC‑GGGTATCCCCATGTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCa                                                          <  2:1243521/140‑1 (MQ=255)                                                                                      tGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCata                                                        >  2:533957/1‑140 (MQ=255)                                                                                       gggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  2:1170340/1‑140 (MQ=255)                                                                                       gggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:893556/1‑140 (MQ=255)                                                                                       gggAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:1243521/1‑140 (MQ=255)                                                                                         gAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt                                                     >  1:768430/1‑140 (MQ=255)                                                                                           aaTTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                                   >  2:1418128/1‑140 (MQ=255)                                                                                                   ggggTACCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                            <  2:1226918/139‑1 (MQ=255)                                                                                                     ggTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa                                          >  2:1364774/1‑140 (MQ=255)                                                                                                        aTCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                       >  1:1221374/1‑140 (MQ=255)                                                                                                        aTCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                       >  1:703719/1‑140 (MQ=255)                                                                                                          cccaaTGTAAAACGTTGTTCCGCGGGGAAGAAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                     <  2:816588/136‑1 (MQ=255)                                                                                                           aCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa                                    >  1:429794/1‑140 (MQ=255)                                                                                                              aTGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                 >  1:1283263/1‑140 (MQ=255)                                                                                                               tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                >  1:689938/1‑140 (MQ=255)                                                                                                               tgtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc                                >  1:1226918/1‑140 (MQ=255)                                                                                                                gtgAAACTTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               <  1:320344/140‑1 (MQ=255)                                                                                                                gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  1:144708/1‑140 (MQ=255)                                                                                                                gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:651317/1‑140 (MQ=255)                                                                                                                ggaaaaTGTTTTTCCGCGGGGAAGAAAGGCCCCTGTAGCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               <  2:689938/137‑1 (MQ=255)                                                                                                                   aagcGTTGTGACGCCGGGGAGAAAAGCCCCTGCAGCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            <  2:837260/137‑1 (MQ=255)                                                                                                                   aaaCTTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            <  1:39902/140‑1 (MQ=255)                                                                                                                    aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc                           >  2:1300925/1‑140 (MQ=255)                                                                                                                       gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                        >  1:1397584/1‑140 (MQ=255)                                                                                                                       gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg                        >  1:1296370/1‑140 (MQ=255)                                                                                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    >  2:427636/1‑140 (MQ=255)                                                                                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCaa                    >  2:1258347/1‑140 (MQ=255)                                                                                                                                ggggggAGAAAAGCCCGTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccg               <  2:1397584/138‑1 (MQ=255)                                                                                                                                ccGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccg               <  1:427636/140‑1 (MQ=255)                                                                                                                                  ggggAAAAAAGCCCCTGTACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccgcc             <  2:1283263/140‑1 (MQ=255)                                                                                                                                   gggaaaaaaGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  2:422538/134‑1 (MQ=255)                                                                                                                                     gAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg          >  2:1195597/1‑140 (MQ=255)                                                                                                                                             aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  1:805745/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                            GCTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628486‑2628764

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