Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,968,542	G→A	6.8%	D224N (GAT→AAT)	PP_3498 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,968,542	0	G	A	6.8%	162.9 / 2.1	59	D224N (GAT→AAT)	PP_3498	transposase
Reads supporting (aligned to +/- strand):  ref base G (29/26);  new base A (2/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CGCGCTGGGTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATCGTGGC  >  NC_002947/3968400‑3968690                                                                                                                                               |                                                                                                                                                     cgcgCTAGGTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGt                                                                                                                                                >  5:28612/1‑149 (MQ=17)      tGGGTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCCTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCt                                                                                                                                            >  4:121457/1‑148 (MQ=17)       gggTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGa                                                                                                                                          >  7:218111/1‑149 (MQ=32)       gggTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGa                                                                                                                                          >  2:79443/1‑149 (MQ=32)               aGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTgg                                                                                                                                  >  4:153890/1‑149 (MQ=32)               aGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTg                                                                                                                                   <  8:218111/148‑1 (MQ=32)               aGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTg                                                                                                                                   <  2:120153/148‑1 (MQ=32)                  tCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGtt                                                                                                                                <  4:180697/148‑1 (MQ=32)                     gTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGtct                                                                                                                            <  2:127311/149‑1 (MQ=32)                      tGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGtctc                                                                                                                           <  8:316762/149‑1 (MQ=32)                        gTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTctct                                                                                                                          <  3:153890/148‑1 (MQ=32)                        gTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTa                                                                                                                         <  6:297395/149‑1 (MQ=32)                        gTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTa                                                                                                                         <  6:28612/149‑1 (MQ=32)                                         ggCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGgc                                                                                                        >  1:115505/1‑149 (MQ=32)                                          gCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTgg                                                                                                                                  <  5:139398/122‑1 (MQ=32)                                          gCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTgg                                                                                                                                  >  6:139398/1‑122 (MQ=32)                                          gCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATa                                                                                                             <  8:120010/143‑1 (MQ=32)                                          gCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATa                                                                                                             >  7:120010/1‑143 (MQ=32)                                                ccATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGAGTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAgc                                                                                                 >  4:237719/1‑149 (MQ=17)                                                   tGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATAGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAgccgc                                                                                              <  1:79443/149‑1 (MQ=17)                                                      cccAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCgg                                                                                           <  3:209532/149‑1 (MQ=32)                                                      cccAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCgg                                                                                           >  3:104779/1‑149 (MQ=32)                                                        cAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCt                                                                                         >  5:274253/1‑149 (MQ=32)                                                        cAGGGCGATGCCGAGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTAGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACTCATAACGGCGCTGAGCCGCCGGCt                                                                                         <  2:123188/149‑1 (MQ=12)                                                              gATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCCTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGc                                                                                   <  3:121457/149‑1 (MQ=17)                                                                   cgcgATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACg                                                                                                          <  6:118563/121‑1 (MQ=32)                                                                   cgcgATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACg                                                                                                          >  5:118563/1‑121 (MQ=32)                                                                          aaTAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGccc                                                                           >  6:363227/1‑145 (MQ=34)                                                                          aaTAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCaa                                                                         <  5:363227/147‑1 (MQ=34)                                                                          aaTAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAaga                                                                       <  2:221825/149‑1 (MQ=34)                                                                             aaGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGa                                                                    >  5:65485/1‑149 (MQ=34)                                                                              aGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGAt                                                                   <  2:181763/149‑1 (MQ=34)                                                                                      ccTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAgg                                                                     >  3:232635/1‑139 (MQ=34)                                                                                      ccTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAgg                                                                     <  4:232635/139‑1 (MQ=34)                                                                                           ggCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGc                                                      <  4:104779/149‑1 (MQ=32)                                                                                                ttGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCAATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGa                                                    <  5:313801/146‑1 (MQ=17)                                                                                                ttGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCAATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGa                                                    <  1:239621/146‑1 (MQ=17)                                                                                                ttGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCAATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGa                                                    >  2:239621/1‑146 (MQ=17)                                                                                                ttGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCAATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGa                                                    >  6:313801/1‑146 (MQ=17)                                                                                                       tGCTTCCTTGCCTCTAACCCACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACAGGTTGTCTCTAGGTGAACGCATAACGGCGCTGATCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCa                                          >  8:81201/1‑149 (MQ=255)                                                                                                        gCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGc                                                                             >  5:93740/1‑113 (MQ=32)                                                                                                        gCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGc                                                                             <  6:93740/113‑1 (MQ=32)                                                                                                        gCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAaga                                                                       >  4:326812/1‑119 (MQ=32)                                                                                                          ttCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTc                                                         <  5:191890/131‑1 (MQ=32)                                                                                                          ttCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTc                                                         >  6:191890/1‑131 (MQ=32)                                                                                                          ttCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGAt                                               <  1:27109/141‑1 (MQ=34)                                                                                                          ttCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGAt                                               >  2:27109/1‑141 (MQ=34)                                                                                                                  ctctACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTc                               >  1:275885/1‑149 (MQ=34)                                                                                                                    ctACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCgg                             >  4:172010/1‑149 (MQ=34)                                                                                                                      aCCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGcc                           >  1:156234/1‑149 (MQ=34)                                                                                                                       cccGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGccc                          >  5:406231/1‑149 (MQ=34)                                                                                                                          gACCTGTTGCGCGAGTGGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTAt                       <  6:233407/149‑1 (MQ=21)                                                                                                                             cTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACAGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCATATCAt                    >  5:43683/1‑149 (MQ=21)                                                                                                                               gTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTg                  >  1:305262/1‑149 (MQ=34)                                                                                                                                ttGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTg                  <  2:156234/148‑1 (MQ=34)                                                                                                                                      gAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACgg            >  3:322780/1‑148 (MQ=34)                                                                                                                                        gTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCaa         >  4:153053/1‑149 (MQ=35)                                                                                                                                            cTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATc       <  6:162363/147‑1 (MQ=37)                                                                                                                                            cTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATc       >  5:162363/1‑147 (MQ=37)                                                                                                                                               gATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATCGTg    <  4:322780/147‑1 (MQ=37)                                                                                                                                               gATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATCGTGGc  <  3:153053/149‑1 (MQ=37)                                                                                                                                               |                                                                                                                                                     CGCGCTGGGTAATCAGATCCGTGGTTTGCTACTTGAGCATGGCGTATCCATGCCCCAGGGCGATGCCGCGATCAATAAGCATGTGCCTCGGGCGCTTGAGGATGCTTCCTTGCCTCTACCCGACCTGTTGCGCGAGTTGCTCGATGAGTTGCTGATGGACTGGTTGTCTCTAGGTGAACGCATAACGGCGCTGAGCCGCCGGCTCGAAGCGACTGCCCAAGAGGATAAAGTCGCTCAGCGACTGATCACCATTCGCGGCGTCGGCCCTATCATTGCCACGGCAATCGTGGC  >  NC_002947/3968400‑3968690

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