Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,651,374	T→G	7.1%	intergenic (‑182/+31)	PP_4959 ← / ← fba	diguanylate cyclase/fructose‑bisphosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,651,374	0	T	G	7.1%	250.4 / 6.3	81	intergenic (‑182/+31)	PP_4959/fba	diguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base G (3/3);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGA  >  NC_002947/5651230‑5651520                                                                                                                                                 |                                                                                                                                                    ttGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCt                                                                                                                                                 <  2:54295/149‑1 (MQ=255)    gATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCg                                                                                                                                              >  4:324235/1‑149 (MQ=255)     aTCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                             <  5:315245/149‑1 (MQ=255)     aTCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                             <  1:179816/149‑1 (MQ=255)      tacATTAGTACTAAAGAACTGAGGCTGCAACACGCAATGCAATCAGAACGATATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTCCCGCAATTATGTACAGGTTTCAATCTCCCATAAAAAAACCCGAAGCGCTGCggc                                                                                                                                            >  8:71133/3‑148 (MQ=255)         aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGtt‑t                                                                                                                                        >  6:268471/1‑149 (MQ=255)         aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGtt‑t                                                                                                                                        >  1:41622/1‑149 (MQ=255)         aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGtt‑t                                                                                                                                        <  2:117242/149‑1 (MQ=255)          ttAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGtt‑tt                                                                                                                                       >  7:280197/1‑149 (MQ=255)          ttAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGtt‑tt                                                                                                                                       >  8:4503/1‑149 (MQ=255)             gTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtt‑ttttt                                                                                                                                    <  8:45404/149‑3 (MQ=255)               aCTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAg                                                                                                                                   <  2:27467/148‑1 (MQ=255)                ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAggg                                                                                                                                 <  5:168410/149‑1 (MQ=255)                ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCAATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAggg                                                                                                                                 >  7:3437/1‑149 (MQ=255)                 tctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGt                                                                                                                                >  4:38723/1‑149 (MQ=255)                 tctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGt                                                                                                                                >  3:79997/1‑149 (MQ=255)                       ccTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCt                                                                                                                           <  4:79997/148‑1 (MQ=255)                         tGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTc                                                                                                                        >  5:30050/1‑149 (MQ=255)                                 ccgttctcaatgcaatCAGTACGCTTTCAATGATCAAACTGTACGGTATCCGACTTACGGTGGGTCTTCCGGAATTATGTACAGTTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGt                                                                                                                 <  2:295646/141‑1 (MQ=255)                                 cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGtt                                                                                                                >  6:128705/1‑149 (MQ=255)                                 cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGtt                                                                                                                >  5:234847/1‑149 (MQ=255)                                     acGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGAcc                                                                                                            >  5:136474/1‑149 (MQ=255)                                      cGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtt‑tttttat                                                                                                                                  <  6:24835/126‑5 (MQ=255)                                      cGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtt‑tttttat                                                                                                                                  >  5:24835/1‑122 (MQ=255)                                       gcaatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCtt                                                                                                          <  7:4503/149‑1 (MQ=255)                                         aatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTgg                                                                                                        >  2:263448/1‑149 (MQ=255)                                             caatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTggcggc                                                                                                    <  5:215150/149‑1 (MQ=255)                                                tCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCAATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAg                                                                                                 >  6:53662/1‑149 (MQ=255)                                                   gAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAgg                                                                                                                                  <  8:191059/113‑1 (MQ=255)                                                   gAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAgg                                                                                                                                  >  7:191059/1‑113 (MQ=255)                                                    aaCGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCa                                                                                             <  1:193575/149‑1 (MQ=255)                                                            cAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAg                                                                                     >  3:189704/1‑149 (MQ=255)                                                             aaTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    >  4:34137/1‑149 (MQ=255)                                                             aaTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    >  6:161482/1‑149 (MQ=255)                                                             aaTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAg                                                                                     <  5:128705/148‑1 (MQ=255)                                                                gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAg                                                                                 >  2:182743/1‑149 (MQ=255)                                                                 aTCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTCTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTa                                                                                  <  1:263448/148‑1 (MQ=255)                                                                   cAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    >  1:126007/1‑143 (MQ=255)                                                                   cAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    <  2:126007/143‑1 (MQ=255)                                                                    aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTc                                                                                       <  5:149291/139‑1 (MQ=255)                                                                    aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTc                                                                                       >  6:149291/1‑139 (MQ=255)                                                                    aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgcgc                                                                              <  2:263926/148‑1 (MQ=255)                                                                     aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    <  4:47561/141‑1 (MQ=255)                                                                     aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    >  3:47561/1‑141 (MQ=255)                                                                          tACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                             >  6:3961/1‑80 (MQ=255)                                                                            cGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACAt                                                                      <  6:30050/148‑1 (MQ=255)                                                                                    gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                                >  2:159090/1‑146 (MQ=255)                                                                                    gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                                >  4:202458/1‑146 (MQ=255)                                                                                     aCTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                            <  1:13031/149‑1 (MQ=255)                                                                                       tAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                            <  3:91747/147‑1 (MQ=255)                                                                                       tAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                              >  4:91747/1‑145 (MQ=255)                                                                                         aCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagaga                                                          <  3:324235/147‑1 (MQ=255)                                                                                         aCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                            >  5:141266/1‑145 (MQ=255)                                                                                               ggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCt                                                  <  8:280197/149‑1 (MQ=255)                                                                                                gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                                 <  5:268471/149‑1 (MQ=255)                                                                                                gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                                 >  5:4624/1‑149 (MQ=255)                                                                                                gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAACGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                                 >  1:169291/1‑149 (MQ=255)                                                                                                 tcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                                 >  8:230893/1‑148 (MQ=255)                                                                                                       ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                          >  5:16913/1‑149 (MQ=255)                                                                                                         aaTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGTCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                        <  1:159090/149‑1 (MQ=255)                                                                                                         aaTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                          <  6:234847/147‑1 (MQ=255)                                                                                                           ttATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAgg                                      >  8:107851/1‑149 (MQ=255)                                                                                                                 aCAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTA‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                >  1:237401/1‑149 (MQ=255)                                                                                                                  cAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                <  7:230893/148‑1 (MQ=255)                                                                                                                  cAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                <  3:34137/148‑1 (MQ=255)                                                                                                                    ggTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcg                             <  6:208824/149‑1 (MQ=255)                                                                                                                          cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                       <  5:15474/149‑1 (MQ=255)                                                                                                                               cccATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAg                  >  1:277981/1‑149 (MQ=255)                                                                                                                                ccATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg                 >  1:210712/1‑149 (MQ=255)                                                                                                                                 cATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGc                >  1:320628/1‑149 (MQ=255)                                                                                                                                  cctgaaaaaCCCGCAGCGATGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCt               >  3:292512/5‑149 (MQ=255)                                                                                                                                  aTAAAAACACCCGCATCGCTGCGGTTT‑TTTCAGGGTTTCCTGTCGAGCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTTGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCt               <  7:71133/149‑1 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                                >  1:269145/1‑116 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                                <  2:269145/116‑1 (MQ=255)                                                                                                                                      aaaaaCCCGCAGCGATGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc             <  4:292512/147‑1 (MQ=255)                                                                                                                                           cccGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgc      <  2:41622/149‑1 (MQ=255)                                                                                                                                             cGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg                 <  6:37705/136‑1 (MQ=255)                                                                                                                                             cGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg                 >  5:37705/1‑136 (MQ=255)                                                                                                                                             cGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc    <  2:277981/149‑1 (MQ=255)                                                                                                                                               cATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGa  <  1:182743/149‑1 (MQ=255)                                                                                                                                               cATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGa  <  2:169291/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTT‑TTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGA  >  NC_002947/5651230‑5651520

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