Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,488,992	C→G	7.6%	intergenic (+345/+226)	PP_3981 → / ← PP_3982	transposase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,488,992	0	C	G	7.6%	213.1 / 6.4	66	intergenic (+345/+226)	PP_3981/PP_3982	transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base G (2/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCAC  >  NC_002947/4488846‑4489131                                                                                                                                                   |                                                                                                                                            gATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCt                                                                                                                                           >  4:33206/1‑149 (MQ=255) gATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCt                                                                                                                                           >  6:398715/1‑149 (MQ=255)    gCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTccc                                                                                                                                        >  2:148758/1‑149 (MQ=255)     cccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCg                                                                                                                                       <  8:10076/149‑1 (MQ=255)     cccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCg                                                                                                                                       >  6:151751/1‑149 (MQ=255)      ccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATTTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGc                                                                                                                                      >  7:44830/1‑149 (MQ=255)      ccAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGc                                                                                                                                      >  6:386657/1‑149 (MQ=255)       cAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCATATGGGCGGGAGGTTCCCGcc                                                                                                                                     >  7:406664/1‑149 (MQ=255)          ggCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTCCTTGCACACGGGCGGGAGGTTCCCGcccgt                                                                                                                                  <  1:37847/149‑3 (MQ=255)            cTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCATATGGGCGGGAGGTTCCCGcccgtgt                                                                                                                                <  8:406664/149‑5 (MQ=255)            cTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAGGTTCCCGcccgtgt                                                                                                                                <  5:345245/149‑5 (MQ=255)             tCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATg                                                                                                                               >  2:165603/1‑149 (MQ=255)                aCATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgcgc                                                                                                                            <  1:14878/149‑1 (MQ=255)                 cATTGCGCCATCACGACTGGAAAATGAAGGTTGTCGAAAATCCACCGATCGCCAAGATCAATCATCGAGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTCAGGCCTTGATTGCACACGGGAGGGACCCTCCCGCCCATATgcgc                                                                                                                            >  4:280045/1‑148 (MQ=255)                     gcgcCATGACGACTGGAAACTGCAGGAAGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTAc                                                                                                                        <  2:81226/148‑1 (MQ=255)                                 cTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgc                                                                                                                              <  7:217503/130‑1 (MQ=255)                                 cTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATgc                                                                                                                              >  8:217503/1‑130 (MQ=255)                                     aaaaTGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATTTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATc                                                                                                       <  8:44830/149‑1 (MQ=255)                                              gTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATa                                                                                                                  >  4:227880/1‑129 (MQ=255)                                              gTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATa                                                                                                                  <  3:227880/129‑1 (MQ=255)                                              gTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTcc                                                                                                 >  5:98878/1‑146 (MQ=255)                                                 gTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGAc                                                                                            <  3:174395/148‑1 (MQ=255)                                                    gAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTca                                                                                        >  4:403113/1‑149 (MQ=255)                                                    gAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTca                                                                                        <  5:339382/149‑1 (MQ=255)                                                      aaaTCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTcaca                                                                                      >  1:276555/1‑149 (MQ=255)                                                         tCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACAtct                                                                                   >  3:152835/1‑149 (MQ=255)                                                          ccACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACAtct                                                                                   >  7:207194/1‑148 (MQ=255)                                                              cGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACtt                                                                              <  3:33206/149‑1 (MQ=255)                                                                     cAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGaa                                                                       >  4:232282/1‑149 (MQ=255)                                                                         aTCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTgg                                                                   <  5:386657/149‑1 (MQ=255)                                                                               cATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTa                                                              >  6:62427/1‑148 (MQ=255)                                                                                  cgcgCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCAc                                                                                >  8:301569/1‑127 (MQ=255)                                                                                  cgcgCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCAc                                                                                <  7:301569/127‑1 (MQ=255)                                                                                     gcCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAActct                                                       <  1:148758/149‑1 (MQ=255)                                                                                       ccAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTcc                                                     <  5:398715/149‑1 (MQ=255)                                                                                           cGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAg                                                 >  5:221205/1‑149 (MQ=255)                                                                                             tGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACACCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGt                                               <  1:165603/149‑1 (MQ=255)                                                                                               gAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTg                                              <  5:151751/148‑1 (MQ=255)                                                                                               gAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGa                                             >  6:169695/1‑149 (MQ=255)                                                                                                  cTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGt                                          <  6:98878/149‑1 (MQ=255)                                                                                                       ggATAAAGGTCTTATCTGATGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGAc                                     <  3:280045/149‑1 (MQ=255)                                                                                                          tAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGcc                                  >  6:371279/1‑149 (MQ=255)                                                                                                            aaGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCt                                 <  6:221205/148‑1 (MQ=255)                                                                                                             aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                               >  7:59829/1‑149 (MQ=255)                                                                                                             aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                               >  2:333649/1‑149 (MQ=255)                                                                                                             aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                               <  3:403113/149‑1 (MQ=255)                                                                                                                    aTCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATc                         >  3:288725/1‑148 (MQ=255)                                                                                                                     tCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCaa                       >  5:432328/1‑149 (MQ=255)                                                                                                                          ggCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCg                  <  6:432328/149‑1 (MQ=255)                                                                                                                           gCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                          >  6:385632/1‑140 (MQ=255)                                                                                                                           gCCTTGCTTGCACACGGGCGCGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                          <  5:385632/140‑1 (MQ=255)                                                                                                                            ccTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAg                <  5:62427/149‑1 (MQ=255)                                                                                                                            ccTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAg                <  4:288725/149‑1 (MQ=255)                                                                                                                             cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg               >  5:348168/1‑149 (MQ=255)                                                                                                                             cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg               <  8:207194/149‑1 (MQ=255)                                                                                                                                  ttgcACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACCTCAATACCGAGGTAAc           >  6:214844/1‑148 (MQ=255)                                                                                                                                     cacaCGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTTGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGt       <  3:232282/149‑1 (MQ=255)                                                                                                                                      acacGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGCGGTAACTAGtt      >  4:159407/1‑149 (MQ=255)                                                                                                                                       cacGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGttt     >  7:259183/1‑149 (MQ=255)                                                                                                                                        acGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGttt     >  4:401350/1‑148 (MQ=255)                                                                                                                                          gggcgggAGGTTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcac  >  5:175632/1‑149 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGcc                              <  7:360260/121‑1 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGcc                              >  8:360260/1‑121 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                          <  5:380224/125‑1 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                          >  6:380224/1‑125 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcac  <  3:28952/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            GATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCAC  >  NC_002947/4488846‑4489131

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