Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,701,987	G→A	6.8%	L139L (CTG→CTA)	PP_3267 →	Clp protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,701,987	0	G	A	6.8%	172.4 / 2.8	59	L139L (CTG→CTA)	PP_3267	Clp protease
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base A (2/2);  total (33/26)
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.78e-01

ATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGT  >  NC_002947/3701842‑3702135                                                                                                                                                  |                                                                                                                                                      atatCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTAt                                                                                                                                                    <  3:290558/149‑1 (MQ=255)        ctTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGCCCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATgcg                                                                                                                                                 >  3:211583/1‑145 (MQ=255)            cTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGt                                                                                                                                         >  5:3407/1‑149 (MQ=255)             ttGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGTACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGAGCCTGTATGCGCAGATCGTg                                                                                                                                        >  4:145463/1‑149 (MQ=255)                ccGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCg                                                                                                                                     >  7:194941/1‑149 (MQ=255)                 cGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCgg                                                                                                                                    <  6:3407/149‑1 (MQ=255)                          aCGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑Cg                                                                                                                           >  6:120297/1‑148 (MQ=255)                            gCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGa                                                                                                                          >  1:59400/1‑147 (MQ=255)                                  cTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGAGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATAGTGGCGGAGCGCA‑CGAAAGGcgcg                                                                                                                  >  5:339718/1‑149 (MQ=255)                                         tttCAGCAAGTTTCTTTTCCCCCCGTTTCATTCGACCATGCATGGTTCGGTGGACCAGTGGCGGAAGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCAGTATGCGCAGATCGAGCCGGAGCGCA‑CGAAAGGCGCGGCCGAgg                                                                                                           >  8:245685/1‑149 (MQ=255)                                         tttCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAgg                                                                                                           >  8:191754/1‑149 (MQ=255)                                             aGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGttt                                                                                                        <  2:4340/148‑1 (MQ=255)                                              gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGTCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTg                                                                                                      >  1:98825/1‑149 (MQ=255)                                              gCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTg                                                                                                      >  1:194272/1‑149 (MQ=255)                                              gCAAGTTTCTTTTCCACCCGTTGCTTTGGACCTTGGATGCTTCGGTGGACCTTTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGCCCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTg                                                                                                      <  6:240750/149‑1 (MQ=255)                                              gCAAGTTTCTTTTCCACCCATTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTg                                                                                                      >  6:88652/1‑149 (MQ=255)                                               cAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTgg                                                                                                     >  5:179743/1‑149 (MQ=255)                                                aattttCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                    <  5:148621/146‑1 (MQ=255)                                                aaGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGa                                                                                                    <  5:148617/149‑1 (MQ=255)                                                   tttCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                                  >  8:278092/1‑148 (MQ=255)                                                   tttCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATg                                                                                                  >  2:32271/1‑148 (MQ=255)                                                     tCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGtt                                                                                                <  2:194272/148‑2 (MQ=255)                                                      cTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTgg                                                                                                     <  8:314816/142‑1 (MQ=255)                                                      cTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTgg                                                                                                     >  7:314816/1‑142 (MQ=255)                                                         ttCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTAcc                                                                                              <  2:98825/146‑1 (MQ=255)                                                                   ttCATTGGACATTTCATGGTTCGGTGTACCATTCGCTGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAAt                                                                                  <  6:339718/148‑1 (MQ=255)                                                                      aTTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGcc                                                                              >  4:49992/1‑149 (MQ=255)                                                                         ggACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTa                                                                            <  5:88652/148‑1 (MQ=255)                                                                               tGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAg                                                                      >  7:39317/1‑148 (MQ=255)                                                                                  aTGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGtttt                                                                                                       <  2:216502/112‑1 (MQ=255)                                                                                  aTGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATa                                                                  >  8:125772/1‑149 (MQ=255)                                                                                  aTGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATa                                                                  >  8:242978/1‑149 (MQ=255)                                                                                                 aTTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAg                                                   >  3:168959/1‑149 (MQ=255)                                                                                                      cgGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                              <  2:59400/149‑1 (MQ=255)                                                                                                        gATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTgg                                            <  8:39317/149‑1 (MQ=255)                                                                                                               gAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACgg                                     <  7:242978/149‑1 (MQ=255)                                                                                                                 gTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGga                                   <  7:94349/149‑1 (MQ=255)                                                                                                                     gtgATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                               >  8:89387/3‑149 (MQ=255)                                                                                                                       gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                               <  6:179792/147‑1 (MQ=255)                                                                                                                       gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATcc                             <  3:229802/149‑1 (MQ=255)                                                                                                                       gATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTCGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGCATACTGGCGCCGCGGGAGGCGCTGGAAAACGGgatgat                               >  5:179792/1‑147 (MQ=255)                                                                                                                        aTGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCa                            >  7:36502/1‑149 (MQ=255)                                                                                                                         tGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCAt                           <  6:179743/149‑1 (MQ=255)                                                                                                                           aGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATg                          <  7:89387/148‑1 (MQ=255)                                                                                                                             gcttgATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCTTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGaaa                       <  7:278092/145‑1 (MQ=255)                                                                                                                             ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAg                                                   <  3:243379/121‑1 (MQ=255)                                                                                                                             ccTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAg                                                   >  4:243379/1‑121 (MQ=255)                                                                                                                              cTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAAt                      >  6:31285/1‑149 (MQ=255)                                                                                                                               tccgccaCGATCTGCGCATACAGGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                              <  6:289636/118‑1 (MQ=38)                                                                                                                               tccgccaCGATCTGCGCATACAGGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCt                                              >  5:289636/7‑124 (MQ=38)                                                                                                                                 gATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcga                   <  3:145463/149‑1 (MQ=255)                                                                                                                                    caCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                <  4:35671/148‑1 (MQ=255)                                                                                                                                    caCGATCTGCGCATACAGGCGCAGGTCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                  >  3:35671/2‑147 (MQ=255)                                                                                                                                     aCGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgacga                <  7:125772/148‑1 (MQ=255)                                                                                                                                       gACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCATCGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATcgac                  >  1:61260/1‑145 (MQ=255)                                                                                                                                       gACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGAt             >  2:265264/1‑149 (MQ=255)                                                                                                                                       gACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGAt             >  6:281293/1‑149 (MQ=255)                                                                                                                                             cgcCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGAt       >  7:25170/1‑149 (MQ=255)                                                                                                                                                cTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATcgc    >  5:232226/1‑149 (MQ=255)                                                                                                                                                  gTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGt  >  3:229901/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      ATATCCTCTTCCTTGCCGGAGAGCGACGCACGGCTTGCAGTTTCAGCAAGTTTCTTTTCCACCCGTTTCATTGGACCTTGCATGGTTCGGTGGACCATTCGCGGATGTCCGAGTATGCGATGAGCCTGGATTACGACCTGCGCCTGTATGCGCAGATCGTGGCGGAGCGCA‑CGAAAGGCGCGGCCGAGGTTTTGGATGTACCGCGTTATCTAATGGCCTATCCAAGGATACTGGCGCCGCGGGAGGCGCTGGAAAACGGGATGATCCATGAAATCGACGAGATGCCGATCGCGT  >  NC_002947/3701842‑3702135

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