Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,488,993	C→T	7.5%	intergenic (+346/+225)	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,993	0	C	T	7.5%	193.1 / 3.3	66	intergenic (+346/+225)	PP_3981/PP_3982	transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (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.78e-01

GATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGCTAA  >  NC_002947/4488846‑4489140                                                                                                                                                    |                                                                                                                                                    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                                                                                         <  7:301569/127‑1 (MQ=255)                                                                                  cgcgCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCAc                                                                                         >  8:301569/1‑127 (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                                        <  3:403113/149‑1 (MQ=255)                                                                                                             aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                                        >  7:59829/1‑149 (MQ=255)                                                                                                             aGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGc                                        >  2:333649/1‑149 (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                         <  4:288725/149‑1 (MQ=255)                                                                                                                            ccTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAg                         <  5:62427/149‑1 (MQ=255)                                                                                                                             cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg                        <  8:207194/149‑1 (MQ=255)                                                                                                                             cttgcttgCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAgg                        >  5:348168/1‑149 (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                                   >  6:380224/1‑125 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACAt                                   <  5:380224/125‑1 (MQ=255)                                                                                                                                          gggcgggAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTcac           <  3:28952/149‑1 (MQ=255)                                                                                                                                                    cTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAATACCTCACATCACACTTGGCTGAACTGGGCCTAAAACTCTCCATAGCTGATGTCAGACGCCTGCAAAATCAATACAGAGGTAACGAGTTTCACCACCGCTaa  >  2:181092/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GATGCCCAAGGCTCGACATTGCGCCATGACGACTGGAAAATGCAGGTTGTCGAAAATCCACCGATCGCCAAGATCAAGCATCGCGCCCAGCGTGGAGCTGAAGGATAAAGGGCTTATCTGAGGCCTTGCTTGCACACGGGCGGGAACCTCCCGCCCATATGCGCGTACTTCATAGCGCGACAATCCACTCCAAGACCTCACATCTCACTTGGCTGAACTGGGCCTAAAACTCTCCATAGGTGATGTCTGACGCCTGCCACATCAATACCGAGGTAACGAGTTTCACCACCGCTAA  >  NC_002947/4488846‑4489140

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