Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	917,209	T→A	6.1%	N120Y (AAT→TAT)	PP_0798 ←	GGDEF domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	917,209	0	T	A	6.1%	198.3 / 2.7	65	N120Y (AAT→TAT)	PP_0798	GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (26/35);  new base A (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ATACACCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGATTGGTCG  >  NC_002947/917061‑917354                                                                                                                                                      |                                                                                                                                                  aTACACCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATAtt                                                                                                                                                   <  4:355773/149‑1 (MQ=255) aTACACCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATAtt                                                                                                                                                   <  4:350090/149‑1 (MQ=255)    cacCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGc                                                                                                                                                 <  2:26431/148‑1 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCct                                                                                                                                               >  5:13171/1‑149 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCct                                                                                                                                               >  6:115517/1‑149 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCct                                                                                                                                               >  3:184879/1‑149 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCct                                                                                                                                               >  7:330111/1‑149 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATAGCct                                                                                                                                               >  3:200939/1‑149 (MQ=255)     acCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATAGCct                                                                                                                                               >  3:195697/1‑149 (MQ=255)         gTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCcttgct                                                                                                                                           <  4:327001/149‑1 (MQ=255)            cAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGt                                                                                                                                        >  3:79472/1‑149 (MQ=255)            cAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTAGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTAGCCAACGATATTGCCTTGCTTGt                                                                                                                                        >  4:120850/1‑149 (MQ=255)              ggTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTcc                                                                                                                                      <  8:71178/149‑1 (MQ=255)               gTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCcttgctt                                                                                                                                          <  2:72161/144‑1 (MQ=255)               gTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCcttgctt                                                                                                                                          >  1:72161/1‑144 (MQ=255)                  gCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGc                                                                                                                                   >  4:164624/1‑148 (MQ=255)                  gCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGc                                                                                                                                   <  1:368345/148‑1 (MQ=255)                   ccAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCca                                                                                                                                 <  2:287810/149‑1 (MQ=255)                   ccAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCca                                                                                                                                 >  4:204450/1‑149 (MQ=255)                    cAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCca                                                                                                                                 >  6:50365/1‑148 (MQ=255)                           ggCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGTCCTTGCGCGGCACGCTACTGGTGGAGTCTCCAACGATATTGCCTTGCTTGTCCAGCCACAGCAc                                                                                                                           <  7:265871/147‑1 (MQ=255)                           ggCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCAc                                                                                                                           <  6:22070/147‑1 (MQ=255)                              cTTGAAAGGGGGGCCCACCAGCAAGCCGAGGGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTc                                                                                                                        >  1:29945/1‑148 (MQ=255)                                   aaGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCgcgc                                                                                                                 <  2:47572/149‑1 (MQ=255)                                    aGGGGGGCCCACCAGCAAGACGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCACGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCGACTGGTGGAGTCGCCCACGATATTGCCTTGATTGTCCAGCCAAAGCACGTCAACGCGCt                                                                                                                >  8:186151/1‑149 (MQ=255)                                      gggggCCCACCAGCAAGCCGA‑TGTTGGCATTGTTCTGGTGTGCCGTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAAAGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTg                                                                                                              <  7:186151/149‑1 (MQ=255)                                      gggggCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTg                                                                                                              >  8:237524/1‑149 (MQ=255)                                       ggggCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGc                                                                                                             <  7:307323/149‑1 (MQ=255)                                            ccaccaGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTc                                                                                                         >  3:70706/1‑148 (MQ=255)                                             caccaGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACCCCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCACCGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTcga                                                                                                       >  8:65956/1‑149 (MQ=255)                                                caGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTcgacg                                                                                                     >  2:199896/1‑148 (MQ=255)                                                   cAAGCCGA‑GGGTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAgg                                                                                                 <  3:204450/149‑1 (MQ=255)                                                     aGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCt                                                                                               >  4:104207/1‑149 (MQ=255)                                                        cGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCa                                                                                            <  4:184879/149‑1 (MQ=255)                                                        cGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCa                                                                                            <  2:335467/149‑1 (MQ=255)                                                        cGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCa                                                                                            <  3:265466/149‑1 (MQ=255)                                                             gTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGa                                                                                        <  7:51742/149‑1 (MQ=255)                                                             gTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGa                                                                                        <  4:79472/149‑1 (MQ=255)                                                                 gCATTGTGCTGGTTGGCCTTGAACACCCCTGTTTCGCCGAAGTTGGCCTTGCGCGGGACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCTCCGCGCTTGCGGTCGACGAAGGCTTCATTGaaaa                                                                                     <  6:2031/148‑1 (MQ=255)                                                                    ttGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAg                                                                                 <  6:298266/149‑1 (MQ=255)                                                                            gtggGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTgg                                                                          <  1:388198/148‑1 (MQ=255)                                                                                  cTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTCTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCg                                                                   <  7:178227/149‑1 (MQ=255)                                                                                     gAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCggg                                                                 <  4:70706/148‑1 (MQ=255)                                                                                       acacCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGGAAAGCAGGCGCTGGCGTATCGGGcc                                                               <  5:50365/148‑1 (MQ=255)                                                                                           cccGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCAt                                                                                           >  1:20613/1‑116 (MQ=255)                                                                                           cccGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCAt                                                                                           <  2:20613/116‑1 (MQ=255)                                                                                            ccGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGAc                                                         >  6:381393/1‑149 (MQ=255)                                                                                             cGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACt                                                        >  2:217536/1‑149 (MQ=255)                                                                                               gTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTc                                                      <  1:199896/149‑1 (MQ=255)                                                                                               gTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATAGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTc                                                      <  4:200939/149‑1 (MQ=255)                                                                                               gTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATAGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTc                                                      <  4:195697/149‑1 (MQ=255)                                                                                                   cGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGGCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCgg                                                   <  5:115517/148‑1 (MQ=255)                                                                                                      cGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCggg                                                            <  6:195974/136‑1 (MQ=255)                                                                                                      cGAAGTTGGCCTTGCGCGGCACGCTACTGGAGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCggg                                                            >  5:195974/1‑136 (MQ=255)                                                                                                               ccTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGCGCCGGGGACTTCCGGCAGTTCCTTGCgg                                      <  1:217536/149‑1 (MQ=255)                                                                                                                         cACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACGGccc                            <  8:330111/149‑1 (MQ=255)                                                                                                                            gCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGcc                             >  2:348172/1‑145 (MQ=255)                                                                                                                                cTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCaa                     <  3:164624/149‑1 (MQ=255)                                                                                                                                      gAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGCAATTCCGGCAGTTCCTTGCGGCCAACCCCCCAGAACAAGCCTTg               >  1:324966/1‑149 (MQ=255)                                                                                                                                        gTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCa             <  6:13171/149‑1 (MQ=255)                                                                                                                                           gCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAgg           >  4:27597/1‑148 (MQ=255)                                                                                                                                            ccAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGa         >  2:387458/1‑149 (MQ=255)                                                                                                                                                cGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGTGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGATTgg     <  3:120850/149‑1 (MQ=255)                                                                                                                                                  atatTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGATTGGTc   <  3:225394/149‑1 (MQ=255)                                                                                                                                                   tatTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGATTGGTCg  <  3:104207/149‑1 (MQ=255)                                                                                                                                                    atTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAgg           >  6:140333/1‑139 (MQ=255)                                                                                                                                                    atTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAgg           <  5:140333/139‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATACACCAGTCCAGGTCGCCAAGCTTGGCCTTGAAAGGGGGGCCCACCAGCAAGCCGA‑GGTTGGCATTGTGCTGGTGGGCCTTGAACACCCCGGTTTCGCCGAAGTTGGCCTTGCGCGGCACGCTACTGGTGGAGTCGCCAACGATATTGCCTTGCTTGTCCAGCCACAGCACGTCACCGCGCTTGCGGTCGACGAAGGCTTCATTGAAAAGCAGGCGCTGGCGTATCGGGCCGGGGACTTCCGGCAGTTCCTTGCGGCCAACCGCCCAGATCAAGCCTTGCAGGGATTGGTCG  >  NC_002947/917061‑917354

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