Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,460,226	Δ1 bp	6.7%	coding (136/1446 nt)	PP_3075 ←	transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,460,226	0	T	.	6.7%	190.3 / 7.5	60	coding (136/1446 nt)	PP_3075	transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base . (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGGCTGCCGGGGATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGAAAACTCATCGTCGGGGGT  >  NC_002947/3460081‑3460371                                                                                                                                                   |                                                                                                                                                   aGGCTGCCGGGGATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCgatgat                                                                                                                                                 <  5:108240/149‑1 (MQ=255)  ggCTGCCGGGGATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGAtc                                                                                                                                                >  5:867/1‑149 (MQ=255)    cTGCCGGGGATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGAtc‑tc                                                                                                                                             <  6:158307/149‑1 (MQ=255)            gATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcg                                                                                                                                     >  8:226388/1‑149 (MQ=255)             aTCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcga                                                                                                                                    <  2:372025/149‑1 (MQ=255)              tCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcgag                                                                                                                                   >  6:188965/1‑149 (MQ=255)                aCTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcgagcg                                                                                                                                 >  5:343745/1‑149 (MQ=255)                aCTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcgagc                                                                                                                                  <  7:226388/148‑1 (MQ=255)                 cTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGA‑GATCATCGAACagtccag                                                                                                                                  <  7:214952/147‑6 (MQ=255)                 cTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGA‑GATCATCGAACagtccag                                                                                                                                  >  8:214952/1‑142 (MQ=255)                    tttCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGATTGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGAt                                                                                                                              <  1:278382/149‑1 (MQ=255)                    tttCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGAt                                                                                                                              <  6:101805/148‑1 (MQ=255)                           cGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGtt                                                                                                                      <  6:199653/149‑1 (MQ=255)                                      tGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTgg                                                                                                           >  1:134726/1‑149 (MQ=255)                                        gCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTgggg                                                                                                         <  5:188965/149‑1 (MQ=255)                                              gCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTc                                                                                                   >  7:117344/1‑149 (MQ=255)                                              gCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTc                                                                                                   <  5:238597/149‑1 (MQ=255)                                              gCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGA‑GATCATCGAACagtc                                                                                                                                     >  1:16475/1‑113 (MQ=255)                                              gCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGA‑GATCATCGAACagtc                                                                                                                                     <  2:16475/115‑3 (MQ=255)                                                 gCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCtt                                                                                                <  6:867/149‑1 (MQ=255)                                                   gTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcgagc                                                                                                                                  <  1:236405/113‑1 (MQ=255)                                                   gTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACagcgagc                                                                                                                                  >  2:236405/1‑113 (MQ=255)                                                   gTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGa                                                                                              <  3:198427/149‑1 (MQ=255)                                                     cGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAgg                                                                                            <  2:134726/149‑1 (MQ=255)                                                       gccaggccaAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGgct                                                                                          >  5:308802/1‑149 (MQ=255)                                                           ggccaAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTTAGgctgtcg                                                                                      <  6:343745/149‑1 (MQ=255)                                                              caAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGtt                                                                                                                      >  6:330147/1‑114 (MQ=255)                                                              caAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCtt                                                                                                <  8:127932/136‑1 (MQ=255)                                                              caAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCtt                                                                                                >  7:127932/1‑136 (MQ=255)                                                              caAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGgctgtcgctg                                                                                   >  7:59028/1‑149 (MQ=255)                                                               aaaGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGgctgtcgctgt                                                                                  >  6:234646/1‑149 (MQ=255)                                                                aatcgtcgtGCGTAGCGTTAGTTCATCCAGACTATGCTTGCCTGGCGGTCGA‑CGATCACCTTGCCCTCTCTGGACTGTTCTATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGTTAGTCCTTGAGgctgtcgctgtc                                                                                 <  2:317247/146‑1 (MQ=255)                                                                aaGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGgctgtcgctgtc                                                                                 <  5:371245/149‑1 (MQ=255)                                                                        ggCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATgg                                                                         <  8:59028/149‑1 (MQ=255)                                                                          cGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGt                                                                       <  6:308802/149‑1 (MQ=255)                                                                                 ttcgttcATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATc                                                                <  8:117344/149‑1 (MQ=255)                                                                                  tcgttcATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGt                                                                       <  5:87156/141‑1 (MQ=255)                                                                                  tcgttcATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGt                                                                       >  6:87156/1‑141 (MQ=255)                                                                                   cgttcATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATcc                                                               <  7:252433/148‑1 (MQ=255)                                                                                           cAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCAt                                                       <  5:234646/148‑1 (MQ=255)                                                                                            aGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCt                                                              <  3:326948/140‑1 (MQ=255)                                                                                            aGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCt                                                              >  4:326948/1‑140 (MQ=255)                                                                                              aCGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTAt                                                   <  3:101155/149‑1 (MQ=255)                                                                                                     gTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTg                                                 <  5:259900/144‑1 (MQ=255)                                                                                                     gTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTg                                                 >  6:259900/1‑144 (MQ=255)                                                                                                           ggCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGaa                                                           >  6:291780/1‑128 (MQ=255)                                                                                                           ggCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGaa                                                           <  5:291780/128‑1 (MQ=255)                                                                                                              ggTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCcgcg                                    >  5:207085/1‑148 (MQ=255)                                                                                                                   a‑cgaTCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCACTTGCCGTACCTGGGGGTAGTCATTGAGGCTGTCGAAGGCGTGCATGGGTGCACATCCTGAATCATGTATTGCCACACGCAAGCGAAAagg                              >  7:165734/1‑149 (MQ=255)                                                                                                                        tCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAaggccag                          >  5:170404/1‑149 (MQ=255)                                                                                                                            cGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGa                                                                                              >  7:294497/1‑77 (MQ=255)                                                                                                                            cGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGa                                                                                              <  8:294497/77‑1 (MQ=255)                                                                                                                            cGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCt                      >  4:176124/1‑149 (MQ=255)                                                                                                                                cccTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAaggcc                            <  6:322490/139‑1 (MQ=255)                                                                                                                                cccTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAaggcc                            >  5:322490/1‑139 (MQ=255)                                                                                                                                cccTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGaa                   >  3:159983/1‑148 (MQ=255)                                                                                                                                 ccTAGCTGGACTGTTCGATTATC‑TCGAACATCGAGCGTATCGCCAGTTGCCGTACCTGGTGTTATTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTCAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGaaaa                 <  7:313420/149‑1 (MQ=255)                                                                                                                                  cTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGAAAAc                >  6:113411/1‑149 (MQ=255)                                                                                                                                            tGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGaa                   <  3:266997/136‑1 (MQ=255)                                                                                                                                            tGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGAAAACTCATCGTCgg      >  1:238865/1‑149 (MQ=255)                                                                                                                                            tGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCACGCGAAAAGGCCAGGCCTGaa                   >  4:266997/1‑136 (MQ=255)                                                                                                                                                 gatgatC‑TCGAACAGCGAGCGGATCGCAAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGAAAACTCATCGTCGGGGGt  >  2:213353/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGGCTGCCGGGGATCACTGTTTCGCACGGCTGGCCGATGGCACTGGCAGCGTCGGCCAGGCCAAAGCGTCGGGCGTAGCGTTCGTTCATCCAGACGATGCGTGCCTGGCGGTCGA‑CGATCACCGTGCCCTCGCTGGACTGTTCGATGATC‑TCGAACAGCGAGCGGATCGCCAGTTGCCGTACCTGGGGGTAGTCCTTGAGGCTGTCGCTGTCGTGCATGGGTGAACATCCTGAATCATGTATTGGCACAAGCCCGCGAAAAGGCCAGGCCTGAAAACTCATCGTCGGGGGT  >  NC_002947/3460081‑3460371

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