Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,761	C→G	5.8%	A266G (GCC→GGC)	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,761	0	C	G	5.8%	222.8 / 3.6	69	A266G (GCC→GGC)	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base G (2/2);  total (33/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 = 7.40e-01

CCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGG  >  NC_002947/3593615‑3593903                                                                                                                                                   |                                                                                                                                               ccTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCg                                                                                                                                              >  7:210966/1‑149 (MQ=255)   tGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTg                                                                                                                                            >  7:319115/1‑149 (MQ=255)    gTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGc                                                                                                                                           >  1:222100/1‑149 (MQ=255)     ttCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCt                                                                                                                                          <  3:282941/149‑1 (MQ=255)     ttCGCCGCTGCCAGCGATGCTGAGGCGGACTACATGGCCAGCTCCCACCGCAAGTGGATGCTGGAGCTGCATGTCGGGCGCTTTGGCTTACTGCCCAAGCCGCAAGAGGGGTACCTCGAACGCTTGCCGGACCACGCGCGCGCCGTGCt                                                                                                                                          >  3:204098/1‑149 (MQ=255)      tCGCCGCTGCCATGGTTGCTGAGGGGGTCTATCTGGCCATCTCTCACCGCAATTGGATTCTGTACCTGCATGTCGGGCGCTTTTGCTTGCTGCCCAATCCGCAAGTGGGTTACCTCGAACGTTTGCCGGACCACGAGCGCGCCGTGCt                                                                                                                                          <  2:310496/148‑1 (MQ=255)        gccgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGa                                                                                                                                       <  7:137087/149‑1 (MQ=255)          cgcTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGAt                                                                                                                                      >  7:265154/1‑148 (MQ=255)            cTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgtggtc                                                                                                                                   <  8:210966/149‑6 (MQ=255)                 aGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGAt                                                                                                                              <  5:346731/149‑1 (MQ=255)                  gCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATg                                                                                                                             >  5:140087/1‑149 (MQ=255)                    gATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGc                                                                                                                           <  7:86262/149‑1 (MQ=255)                     aTGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGcc                                                                                                                          >  5:82580/1‑149 (MQ=255)                            ggCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGTCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCg                                                                                                                   >  6:150541/1‑149 (MQ=255)                               ggACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCAGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGACCGCTTGCCGGCCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAg                                                                                                                >  5:316504/1‑149 (MQ=255)                                     ccTGGCCCGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGa                                                                                                           <  8:265154/148‑1 (MQ=255)                                     ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGc                                                                                                                           <  7:161818/132‑1 (MQ=255)                                     ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGc                                                                                                                           >  8:161818/1‑132 (MQ=255)                                     ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAg                                                                                                          >  4:226675/1‑149 (MQ=255)                                     ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCgt                                                                                                                                       <  6:308492/120‑2 (MQ=255)                                     ccTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCACGGCGCGCTCGa                                                                                                                                       >  5:308492/1‑120 (MQ=255)                                      cTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAgg                                                                                                         <  3:101056/149‑1 (MQ=255)                                       tGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTATGGCTTGCTGCCCAAGCCGCAAAAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAggg                                                                                                        >  1:254203/1‑149 (MQ=255)                                         gCCACCTCCCACCGCAAGTGTATGCTGGACCTGCATGTCGGGCGCTTTGGCTTACTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCTCGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCTGAggg                                                                                                        <  4:204098/147‑1 (MQ=255)                                                  cACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGt                                                                                             >  4:133810/1‑149 (MQ=255)                                                     cGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAg                                                                                          <  8:319115/149‑1 (MQ=255)                                                         ggtgggagCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGaa                                                                                      <  6:225706/142‑1 (MQ=255)                                                             gATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGc                                                                                  >  5:107746/1‑149 (MQ=255)                                                                  tGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTgcgc                                                                             >  5:48020/1‑149 (MQ=255)                                                                       cTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGTCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGcgctcgct                                                                        <  5:150541/149‑1 (MQ=255)                                                                         gCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCCGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGa                                                                      >  4:309293/1‑149 (MQ=255)                                                                           aTGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATc                                                                    <  3:133810/149‑1 (MQ=255)                                                                             gTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATc                                                                    <  5:304654/147‑1 (MQ=255)                                                                             gTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATc                                                                    >  6:304654/1‑147 (MQ=255)                                                                                    gcTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGAcc                                                           <  6:48020/149‑1 (MQ=255)                                                                                      tttGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACgg                                                                                               <  1:9431/111‑1 (MQ=255)                                                                                      tttGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACgg                                                                                               >  2:9431/1‑111 (MQ=255)                                                                                       ttGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCggg                                                        >  3:309447/1‑149 (MQ=255)                                                                                          gCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCGTGCCGGACCACGAGCGCGCGGGGCTCGATCAGGGGTTGGCCGGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGc                                                      <  6:316504/148‑1 (MQ=255)                                                                                            ttGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGACCGCGCCCTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGTGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGa                                                   >  1:307795/1‑149 (MQ=255)                                                                                              gctgcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGt                                                                                             <  8:237320/105‑1 (MQ=255)                                                                                              gctgcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGAc                                                  <  6:82580/148‑1 (MQ=255)                                                                                                 gcCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGc                                              <  3:226675/149‑1 (MQ=255)                                                                                                  cccAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCGCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAg                                               <  3:307761/147‑1 (MQ=255)                                                                                                    cAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtga                                           <  6:140087/149‑1 (MQ=255)                                                                                                     aaGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtgat                                          >  7:85505/1‑149 (MQ=255)                                                                                                       gCCTCAAGAGTGTTACCTCGTACGCTTGCCGGACCACGAGCGCTCCTTGGTCGATCAGGTTATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGATTGAAGTGCTGCGCTCGCTTATGGAGCAGACCGGGGCGGACGAGGtgatga                                        <  2:307795/149‑1 (MQ=255)                                                                                                         cGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATc                                      <  4:309447/149‑1 (MQ=255)                                                                                                                       cTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCa                        >  6:292007/1‑149 (MQ=255)                                                                                                                       cTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCa                        >  5:281387/1‑149 (MQ=255)                                                                                                                        tCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCa                                                               >  2:210372/1‑109 (MQ=255)                                                                                                                        tCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCa                                                               <  1:210372/109‑1 (MQ=255)                                                                                                                        tCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCAt                       >  7:319278/1‑149 (MQ=255)                                                                                                                          gAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCTGGTGATGGCCTGCTCCGTCGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCAt                       <  2:222100/147‑1 (MQ=255)                                                                                                                           aaCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtgatg                                         <  5:329934/128‑1 (MQ=255)                                                                                                                           aaCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtgatg                                         <  5:325009/128‑1 (MQ=255)                                                                                                                           aaCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtgatg                                         >  6:329934/1‑128 (MQ=255)                                                                                                                           aaCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCtgatg                                         >  6:325009/1‑128 (MQ=255)                                                                                                                            aCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGAcc                   >  7:299459/1‑149 (MQ=255)                                                                                                                            aCGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGAcc                   <  5:292007/149‑1 (MQ=255)                                                                                                                              gCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACg                                                                                                <  6:106211/70‑1 (MQ=255)                                                                                                                              gCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACg                                                                                                >  5:106211/1‑70 (MQ=255)                                                                                                                                     ggACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgc          <  5:353482/149‑1 (MQ=255)                                                                                                                                      gACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGgcgtg         >  8:69194/1‑147 (MQ=255)                                                                                                                                       aCCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCgg        <  6:281387/149‑1 (MQ=255)                                                                                                                                         cACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCgggg                                                       >  7:243167/1‑100 (MQ=255)                                                                                                                                         cACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCgggg                                                       <  8:243167/100‑1 (MQ=255)                                                                                                                                            gAGCGCGCCGTGCTCGATCAGGTGATGGCCCGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCg   <  3:309293/149‑1 (MQ=255)                                                                                                                                             aGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCgg  >  6:24780/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCTGTTCGCCGCTGCCAGCGATGCTGAGGCGGACTACCTGGCCAGCTCCCACCGCAAGTGGATGCTGGACCTGCATGTCGGGCGCTTTGGCTTGCTGCCCAAGCCGCAAGAGGGTTACCTCGAACGCTTGCCGGACCACGAGCGCGCCGTGCTCGATCAGGTGATGGCCTGCTCCGTAGCCGGAGGGCCCGAACGGGTGAGGGAAGGGCTGCGCTCGCTGATCGAGCAGACCGGGGCGGACGAGCTGATGATCGACTGCCGCATCCATGACCCGCTGGCGCGGCAGCGG  >  NC_002947/3593615‑3593903

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