Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,162,510	G→A	6.4%	P255L (CCA→CTA)	PP_5404 ←	TnsC‑like transposition protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,162,510	0	G	A	6.4%	180.7 / 2.2	62	P255L (CCA→CTA)	PP_5404	TnsC‑like transposition protein
Reads supporting (aligned to +/- strand):  ref base G (24/34);  new base A (2/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

GCAAGCTTCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAACAGTATGTCGGTCA  >  NC_002947/6162363‑6162657                                                                                                                                                      |                                                                                                                                                    gCAAGCTTCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTgg                                                                                                                                                    <  8:307504/149‑1 (MQ=255) gCAAGCTTCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTg                                                                                                                                                     >  1:232850/1‑148 (MQ=255) gCAAGCTTCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTg                                                                                                                                                     >  2:270909/1‑148 (MQ=255)    aGCTTCAGCAGTTGTCCGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATc                                                                                                                                                 <  2:342406/149‑1 (MQ=255)       ttCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATcc                                                                                                                                                >  4:348550/1‑147 (MQ=255)            cagTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCTCCAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCggg                                                                                                                                           >  1:88551/1‑148 (MQ=255)              gTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGgcagc                                                                                                                                       <  3:133586/149‑1 (MQ=255)              gTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGgcag                                                                                                                                        >  7:165345/1‑148 (MQ=255)                tGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGgcagca                                                                                                                                      <  3:67162/148‑1 (MQ=255)                    cGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATc                                                                                                                                 <  5:190918/149‑1 (MQ=255)                    cGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATc                                                                                                                                 <  5:6163/149‑1 (MQ=255)                      aCAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCt                                                                                                                               <  3:71493/149‑1 (MQ=255)                           ccGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGtt                                                                                                                          <  5:214733/149‑1 (MQ=255)                             gCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTgg                                                                                                                        <  7:36034/149‑1 (MQ=255)                              cTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTggg                                                                                                                       <  3:349561/149‑1 (MQ=255)                                    aCCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTCGCTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAag                                                                                                                   >  1:268028/1‑148 (MQ=255)                                    aCCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAag                                                                                                                   >  3:50315/1‑147 (MQ=255)                                     ccAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAagca                                                                                                                 <  2:88551/148‑1 (MQ=255)                                      cAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTGCTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAagcagc                                                                                                               >  3:158440/1‑149 (MQ=255)                                      cAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAagcagc                                                                                                               >  3:290732/1‑149 (MQ=255)                                       aaaaGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGcc                                                                                                              <  2:268028/149‑1 (MQ=255)                                        aaaGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCct                                                                                                             >  3:75462/1‑149 (MQ=255)                                              gAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCtt                                                                                                       >  5:154016/1‑149 (MQ=255)                                                   gTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCt                                                                                                  >  5:229698/1‑149 (MQ=255)                                                    tAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCt                                                                                                  >  1:206217/1‑148 (MQ=255)                                                      aaGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCAc                                                                                               <  4:290732/149‑1 (MQ=255)                                                        gCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTg                                                                                             <  1:270909/149‑1 (MQ=255)                                                             aGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGTCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTg                                                                                          <  1:173045/148‑1 (MQ=255)                                                                cGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGcgac                                                                                      <  3:348550/149‑1 (MQ=255)                                                                  cAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGcgacg                                                                                     >  5:243/1‑148 (MQ=255)                                                                         gTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCt                                                                             <  6:351647/149‑1 (MQ=255)                                                                            cTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTtgt                                                                          >  7:200060/1‑149 (MQ=255)                                                                                   tGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGccc                                                                      >  3:354536/1‑147 (MQ=255)                                                                                     aaCTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTGCTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGa                                                                    <  4:158440/147‑1 (MQ=255)                                                                                      aCTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAgg                                                                  <  2:206217/148‑1 (MQ=255)                                                                                      aCTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAgg                                                                  <  6:243/148‑1 (MQ=255)                                                                                       cTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTATTCCTCACTCTTGCGACGAAACAGCTTGTGCCCGAGGc                                                                 >  8:271277/1‑148 (MQ=255)                                                                                        tGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCta                                                               >  4:10663/1‑149 (MQ=255)                                                                                         gATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCtat                                                              <  8:165345/149‑1 (MQ=255)                                                                                         gATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTAGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAAc                                                                                 <  4:202363/130‑1 (MQ=255)                                                                                         gATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTAGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAAc                                                                                 <  4:197559/130‑1 (MQ=255)                                                                                         gATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTAGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAAc                                                                                 >  3:197559/1‑130 (MQ=255)                                                                                         gATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTAGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAAc                                                                                 >  3:202363/1‑130 (MQ=255)                                                                                            aCGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGt                                                           >  3:315289/1‑149 (MQ=255)                                                                                              gACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTg                                                         >  5:128530/1‑149 (MQ=255)                                                                                                    cACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGAcaca                                                   <  8:200060/149‑1 (MQ=255)                                                                                                    cACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCGCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGAcaca                                                   <  6:229698/149‑1 (MQ=255)                                                                                                     aCTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACaa                                                  <  4:359221/149‑1 (MQ=255)                                                                                                            cTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAAtcgt                                           <  2:16160/149‑1 (MQ=255)                                                                                                            cTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAAtcgt                                           <  6:23700/149‑1 (MQ=255)                                                                                                               aTGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTCTCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAAtcgtc                                          <  2:104410/147‑1 (MQ=255)                                                                                                               aTGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAAtcgtc                                          <  4:315289/147‑1 (MQ=255)                                                                                                               aTGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAAtcgtc                                          <  6:154016/147‑1 (MQ=255)                                                                                                                   aCTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATACTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGAAGAAACAGCTTGTGCCCGAGGATATAGTTGGACACAAAAATCGTCGGCAc                                     >  6:331218/1‑148 (MQ=255)                                                                                                                        ccTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGac                               <  4:354536/149‑2 (MQ=255)                                                                                                                                     cccAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATc                     <  4:50315/146‑1 (MQ=255)                                                                                                                                     cccAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCg                    <  3:10663/147‑1 (MQ=255)                                                                                                                                     cccAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGt                   >  4:144369/1‑148 (MQ=255)                                                                                                                                       cAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCt                                                                                                  <  7:79119/67‑1 (MQ=255)                                                                                                                                       cAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCt                                                                                                  >  8:79119/1‑67 (MQ=255)                                                                                                                                        aaTCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCACTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAAc               <  2:232850/149‑1 (MQ=255)                                                                                                                                        aaTCCTGGCTCCTTGGATCCTCGGGCAGCATTATACTGGGTTGGGAAAGCACCCTATGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCTAGGCTATAGTTTGACACAAAAATCGTCTGCACTCCTATGGCTACCATCGTCaa                >  8:74753/1‑148 (MQ=255)                                                                                                                                                cTCCTTGTATCCTCGGGCAGCATTATCCTGGGTTTGGAAAGCAGCCTCTGCTTTTCCTCACTGTTGCGACGAAACAGCTTGTTCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAACAGTATGTc       <  5:331218/149‑1 (MQ=255)                                                                                                                                                    ttGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAACAGTATGTCGGt    >  3:289318/1‑148 (MQ=255)                                                                                                                                                     tGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAACAGTATGTCGGTCa  >  8:35274/1‑149 (MQ=255)                                                                                                                                                      ggATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGacac                                                    >  4:344330/1‑98 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCAAGCTTCAGCAGTTGTACGACAAGCCGCTTGATACCAAAAGTGGAACGGTAAAGCT‑CAGCCGCAAATTCGTCCTTTC‑CTGAACTGATACGACCGCCACTGAGCCTCATGCACTCGCCTACGTACTCTGCCCAATCCTGGCTCCTTGGATCCTCGGGCAGCATTATCCTGGGTTGGGAAAGCAGCCTCTGCTTGTCCTCACTGTTGCGACGAAACAGCTTGTGCCCGAGGCTATAGTTGGACACAAAAATCGTCGGCACTCCGATGGCTACCATCGTCAACAGTATGTCGGTCA  >  NC_002947/6162363‑6162657

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