Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,167,904	A→T	8.0%	P328P (CCA→CCT) ‡	PP_5407 →	TnsD‑like transposition protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,167,904	0	A	T	8.0%	179.0 / 4.7	62	P328P (CCA→CCT) ‡	PP_5407	TnsD‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base A (35/22);  new base T (3/2);  total (38/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CCCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167759‑6168050                                                                                                                                                   |                                                                                                                                                   cccAAATGACCCGCAAGCCCCGTGGTTACTACCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGc                                                                                                                                                  <  4:54292/148‑1 (MQ=255)    aaaTGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAgg                                                                                                                                              >  7:91102/1‑149 (MQ=255)          cccGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAAc                                                                                                                                        >  8:62602/1‑149 (MQ=255)             gCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACtt                                                                                                                                      >  2:55295/1‑148 (MQ=255)             gCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTta                                                                                                                                     <  7:62602/149‑1 (MQ=255)               aatccccGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTTACTAGCCTCCCAGCGAGGCTAAACTtaac                                                                                                                                   <  1:235215/146‑1 (MQ=255)                             cTGCCATCCGTTGAAGCATCTGAGTCTTATTTGGGGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCATCAGCAACTTCCCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACActttctt                                                                                                                       <  1:294006/148‑1 (MQ=255)                                   tCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTa                                                                                                               >  3:333260/1‑149 (MQ=255)                                      gTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTAttt                                                                                                            <  8:30205/149‑1 (MQ=255)                                          aaGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                        <  6:118635/149‑1 (MQ=255)                                           aGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGc                                                                                                        >  3:138159/1‑148 (MQ=255)                                            gCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGcc                                                                                                       >  4:55823/1‑148 (MQ=255)                                                 tGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaagg                                                                                                  >  8:267587/1‑148 (MQ=255)                                                     tCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCaaggaag                                                                                               >  6:71947/1‑147 (MQ=255)                                                       tGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACGTTCTTTAAGCCTATTTCAGCCCAAGGAAGgag                                                                                            >  8:223625/1‑148 (MQ=255)                                                        gATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGaa                                                                                          >  8:216431/1‑149 (MQ=255)                                                         aTTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGgag                                                                                            <  4:138159/146‑1 (MQ=255)                                                                       cGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGaaa                                                                            <  5:10017/148‑1 (MQ=255)                                                                        gggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAAcc                                                                          >  8:32694/1‑149 (MQ=255)                                                                         ggAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAGcc                                                                          >  8:164607/1‑148 (MQ=255)                                                                            gCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                    >  5:304625/1‑103 (MQ=255)                                                                            gCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTaa                                                                                                                    <  6:304625/103‑1 (MQ=255)                                                                              ttAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGCCTAAACTTAACTAACACTTTCTTTAAGAATATTTCAGCCCAAGGATGGAGAACTGGTGCCCCGAAACCCaaaa                                                                     >  2:139661/1‑148 (MQ=255)                                                                               tAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTCTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCaaaaa                                                                    >  8:59426/1‑148 (MQ=255)                                                                                 gAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAAc                                                                  <  8:91102/148‑1 (MQ=255)                                                                                   aTCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTa                                                                <  4:333260/148‑1 (MQ=255)                                                                                       ttCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTaa                                                           >  4:178948/1‑149 (MQ=255)                                                                                       atCGTTAAGGCCTATCAACATGCCGAGAACCAGCAGCAACTT‑CCAGTGACAAGCCTCCCAGCGATGCTAAACTTAAATAACACTTTCTTTAAGCCTATTTCAGACCAAGGAACGAGAACTGTTGAACCCAACCCCAAAAAACTATTTaa                                                           >  1:393625/2‑149 (MQ=255)                                                                                         cGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAAt                                                          <  7:267587/148‑1 (MQ=255)                                                                                          gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGa                                                        >  2:34287/1‑149 (MQ=255)                                                                                          gTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTGCTTTAAGCCTATTTCATCCCAAGGAAGGAGAACTGGTGCCCCGACCCCCAAAAAACTATTTAATg                                                         >  8:370125/1‑148 (MQ=255)                                                                                           ttAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGAc                                                       <  5:71947/149‑1 (MQ=255)                                                                                           ttAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACCTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGAc                                                       <  7:223625/149‑1 (MQ=255)                                                                                             aaGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                     >  7:137815/1‑149 (MQ=255)                                                                                              aGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                     >  5:225809/1‑148 (MQ=255)                                                                                              aGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACtt                                                     >  5:220191/1‑148 (MQ=255)                                                                                                gCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGaa                                                  <  3:55823/149‑1 (MQ=255)                                                                                                  cTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCACCGAGGCTAAACTTAACTAACACTTTCTTTAACCCTATTTCAGCCCAAGGAACGAGAACTCGTCCCCCGAAACCCAAAAAACAATTTAATCACATCAAg                                                 >  4:126397/1‑148 (MQ=255)                                                                                                     tCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGgaga                                             <  4:47498/149‑1 (MQ=255)                                                                                                             ccGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAAtcctc                                     <  7:32694/149‑1 (MQ=255)                                                                                                                cgAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                 <  5:402551/134‑1 (MQ=255)                                                                                                                cgAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAg                                                 >  6:402551/1‑134 (MQ=255)                                                                                                                   aCCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGaa                               <  8:137815/149‑1 (MQ=255)                                                                                                                    acatcatcaACTT‑CCATTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAAc                              <  3:12866/142‑1 (MQ=255)                                                                                                                     cagcagCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACt                             >  8:317594/1‑149 (MQ=255)                                                                                                                         agcaACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACtgat                          >  7:227671/1‑148 (MQ=255)                                                                                                                               tt‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAAc                              >  1:408945/1‑138 (MQ=255)                                                                                                                               tt‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAAc                              <  2:408945/138‑1 (MQ=255)                                                                                                                               tt‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGgc                   >  2:311977/1‑149 (MQ=255)                                                                                                                                    aGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCa               >  6:277226/1‑149 (MQ=255)                                                                                                                                     aagtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAg              >  4:376038/6‑149 (MQ=255)                                                                                                                                       gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAAc                              >  2:134918/4‑131 (MQ=255)                                                                                                                                       gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAAc                              <  1:134918/128‑1 (MQ=255)                                                                                                                                       gtttAGCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGc             >  7:143804/4‑148 (MQ=255)                                                                                                                                        aCTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCa            >  6:97693/1‑148 (MQ=255)                                                                                                                                          tAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaa          >  4:56738/1‑148 (MQ=255)                                                                                                                                            gCCTCGCTGGGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCaaaaaa       <  8:143804/149‑1 (MQ=255)                                                                                                                                                cccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg   >  2:337697/1‑149 (MQ=255)                                                                                                                                                cccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg   >  2:74092/1‑149 (MQ=255)                                                                                                                                                cccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGAGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg   >  2:216523/1‑149 (MQ=255)                                                                                                                                                 ccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACg   >  2:134319/1‑148 (MQ=255)                                                                                                                                                 ccAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGt  <  5:97693/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCCAAATGACCCGCAAGCCCCGTGGTTACTGCCATCCGTTGAAGCATCTGAGTCTGATTTGGTGGCTATTCGGGAGCTTAGAATCATTCGTTAAGGCCTATCAACATGCCGCGAACCAGCAGCAACTT‑CCAGTGACTAGCCTCCCAGCGAGGCTAAACTTAACTAACACTTTCTTTAAGCCTATTTCAGCCCAAGGAAGGAGAACTGGTGCCCCGAAACCCAAAAAACTATTTAATGACTTGAAGGAGAAAATCCTCCGCGAACTGATGACTGGCGCCAGCAAAAAAGACGT  >  NC_002947/6167759‑6168050

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