Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,651,377	C→A	9.4%	intergenic (‑185/+28)	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,377	0	C	A	9.4%	193.8 / 7.1	68	intergenic (‑185/+28)	PP_4959/fba	diguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base A (5/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

GGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGA  >  NC_002947/5651232‑5651520                                                                                                                                                  |                                                                                                                                                  ggATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTg                                                                                                                                                 <  3:58162/148‑1 (MQ=255) ggATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATAAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCAGGAATTATGTACAGGTTTCCATATCCCATAAAAAAACCCGCATCGCTGc                                                                                                                                                >  7:138057/1‑149 (MQ=255)  gATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGGAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGGCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGc                                                                                                                                                <  3:264456/148‑1 (MQ=255)     ccATTAGTAATCTAGACCTGAGGCTCCAACCCGCAATGCAATCACAACGCTAACAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGCAATTATGTACAGGTATCGATCTCCCAAAAAAAAACCCTCATCGATGCGGG‑t                                                                                                                                           >  2:258235/1‑149 (MQ=255)       aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑ttt                                                                                                                                         <  6:25188/149‑1 (MQ=255)            tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑g                                                                                                                                   >  3:70000/1‑149 (MQ=255)             aCTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑g                                                                                                                                   <  3:26926/148‑1 (MQ=255)              ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑gg                                                                                                                                  <  8:187981/148‑1 (MQ=255)              ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGG‑tttttttatg                                                                                                                                  <  7:89047/149‑6 (MQ=255)                  aGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTAc                                                                                                                             >  3:306057/1‑149 (MQ=255)                   gACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTAcc                                                                                                                            >  1:15467/1‑149 (MQ=255)                   gACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTAc                                                                                                                             <  7:180777/148‑1 (MQ=255)                          ggCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGa                                                                                                                      <  4:197546/148‑1 (MQ=255)                            cTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCa                                                                                                                   <  6:46328/149‑1 (MQ=255)                               cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAg                                                                                                                  <  2:278015/147‑1 (MQ=255)                               cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGtt                                                                                                                >  4:266661/1‑149 (MQ=255)                                aaCACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTATCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGtt                                                                                                                >  5:105856/1‑148 (MQ=255)                                  cacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGAc                                                                                                             <  4:31159/149‑1 (MQ=255)                                  cacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGAc                                                                                                             <  7:168878/149‑1 (MQ=255)                                   acGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGAcc                                                                                                            >  6:113205/1‑149 (MQ=255)                                            aatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTggcggc                                                                                                    >  5:110723/1‑148 (MQ=255)                                                   aCGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCAc                                                                                            <  8:140603/149‑1 (MQ=255)                                                    cGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCAcc                                                                                           >  2:57923/1‑149 (MQ=255)                                                     gCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCt                                                                                          >  8:183478/1‑149 (MQ=255)                                                         tCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTc                                                                                       <  6:105856/148‑1 (MQ=255)                                                            aTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                    >  5:259980/1‑148 (MQ=255)                                                             tGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCttt                                                                                        >  6:172014/1‑143 (MQ=255)                                                             tGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCttt                                                                                        <  5:172014/143‑1 (MQ=255)                                                             tGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTa                                                                                  <  5:113205/149‑1 (MQ=255)                                                              gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAg                                                                                 >  3:32793/1‑149 (MQ=255)                                                                 cAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTgg                                                                                                        <  8:153729/123‑1 (MQ=255)                                                                 cAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTgg                                                                                                        >  7:153729/1‑123 (MQ=255)                                                                   aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTAc                                                                                                                             <  4:157310/100‑1 (MQ=255)                                                                   aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTAc                                                                                                                             >  3:157310/1‑100 (MQ=255)                                                                   aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTg                                                                            <  6:299172/149‑1 (MQ=255)                                                                        tACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTg                                                                                                               <  1:134276/109‑1 (MQ=255)                                                                        tACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTg                                                                                                               >  2:134276/1‑109 (MQ=255)                                                                        tACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTa                                                                         <  4:70000/147‑1 (MQ=255)                                                                           ggTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGc                                                                    >  3:217889/1‑149 (MQ=255)                                                                                ccGACTAACGGTCGGTCTGCCTGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCTCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTcc                                                               <  1:160165/149‑1 (MQ=255)                                                                                  gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                                >  2:295827/1‑146 (MQ=255)                                                                                          gtcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGAc                                                                                                             >  2:282371/1‑93 (MQ=255)                                                                                              gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGAt                                                                                                                     >  6:237380/1‑81 (MQ=255)                                                                                              gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGAt                                                                                                                     <  5:237380/81‑1 (MQ=255)                                                                                              gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                                 <  3:233603/149‑1 (MQ=255)                                                                                               tcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                                <  2:15467/149‑1 (MQ=255)                                                                                                   ccGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGa                                               <  6:110723/146‑1 (MQ=255)                                                                                                     ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                          >  8:34497/1‑149 (MQ=255)                                                                                                       aaTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                        <  7:183478/149‑1 (MQ=255)                                                                                                        aTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGa                                        <  1:295827/148‑1 (MQ=255)                                                                                                        aTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                       <  4:217889/149‑1 (MQ=255)                                                                                                               aCAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                >  5:122542/1‑149 (MQ=255)                                                                                                                   gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                       <  2:318732/138‑1 (MQ=255)                                                                                                                   gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                       >  1:318732/1‑138 (MQ=255)                                                                                                                             cccATAAAAAAACCCGCATCGCTGGGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGAACTTGGCGGCCAGTTCACCTTTCAGCTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATATTCGAGGCATTGCCAGCGGTGCCGAAg                  >  7:154171/1‑149 (MQ=255)                                                                                                                              ccATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                       >  1:270664/1‑143 (MQ=255)                                                                                                                              ccATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                       <  2:270664/143‑1 (MQ=255)                                                                                                                              ccATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg                 <  7:34497/149‑1 (MQ=255)                                                                                                                              accctgaaaaaCCCGCAGCGATGCGGGTTTTTTTATGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                       <  7:15058/139‑1 (MQ=38)                                                                                                                              accctgaaaaaCCCGCAGCGATGCGGGTTTTTTTATGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                       >  8:15058/7‑145 (MQ=38)                                                                                                                                  tgaaaaaCCCGCAGCGATGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc             >  3:6488/3‑149 (MQ=255)                                                                                                                                    aaaaaCCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                >  5:175059/1‑128 (MQ=255)                                                                                                                                    aaaaaCCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                                <  6:175059/128‑1 (MQ=255)                                                                                                                                    aaaaaCCCGCAGCGATGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg            >  2:224866/1‑148 (MQ=255)                                                                                                                                      aaaCCCTCATCGATGCGGG‑TTTTTCA‑GGGTTACCTGTCGATAAGTTGACCTTGGCGGCCAGTTAACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGATTCgga          >  4:211397/1‑146 (MQ=255)                                                                                                                                      aaaCCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTaa         >  4:107247/1‑149 (MQ=255)                                                                                                                                        aCCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGAGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcg       >  5:54763/1‑149 (MQ=255)                                                                                                                                             cATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGa  <  6:259980/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGG‑TTTTTCA‑GGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGCGA  >  NC_002947/5651232‑5651520

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