Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,832,060	A→T	6.9%	intergenic (+35/+117)	pvdS → / ← pvdY	alternative sigma factor/hydroxyproline acetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,832,060	0	A	T	6.9%	179.2 / 2.8	58	intergenic (+35/+117)	pvdS/pvdY	alternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAA  >  NC_002947/4831913‑4832203                                                                                                                                                      |                                                                                                                                                gCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGa                                                                                                                                                 <  5:196518/148‑1 (MQ=255) gCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGAt                                                                                                                                                <  5:93646/149‑1 (MQ=255)      ccGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGctc                                                                                                                                           >  6:328225/1‑149 (MQ=255)        gCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCtctc                                                                                                                                         >  6:112473/1‑149 (MQ=255)         cAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCg                                                                                                                                        >  5:254252/1‑149 (MQ=255)          aGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCGAGAGCGATCCAGGTTCGCTCTGGc                                                                                                                                       <  6:7615/149‑1 (MQ=39)             aGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTTCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCtt                                                                                                                                     <  1:148688/148‑1 (MQ=255)               gACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGATCTACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCtttt                                                                                                                                   <  1:230943/149‑1 (MQ=255)                  aTTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTttgctt                                                                                                                               >  1:289395/1‑149 (MQ=255)                   ttGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCAAGAGCGAACCTGGATCGCTCTCGCTTttgcttg                                                                                                                              <  6:128509/149‑1 (MQ=255)                         aGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcgc                                                                                                                        <  7:156054/149‑1 (MQ=255)                               tGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTAt                                                                                                                  >  4:39668/1‑149 (MQ=255)                                 ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                                         >  4:45680/1‑140 (MQ=255)                                 ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                                         <  3:45680/140‑1 (MQ=255)                                       tCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGtttt                                                                                                          >  5:74027/1‑149 (MQ=255)                                            gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCa                                                                                                                            <  7:19990/126‑1 (MQ=255)                                            gACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GAAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCa                                                                                                                            >  8:19990/1‑126 (MQ=255)                                              ccTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                                        <  7:89785/144‑1 (MQ=255)                                              ccTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                                        >  8:89785/1‑144 (MQ=255)                                               cTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATAc                                                                                                  >  6:39725/1‑149 (MQ=255)                                                   gTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGc                                                                                               <  6:137011/148‑1 (MQ=255)                                                      aaCTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCaa                                                                                            <  3:169541/148‑1 (MQ=255)                                                         ttCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAAc                                                                                        >  3:74260/1‑149 (MQ=255)                                                          tCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACt                                                                                       >  2:59398/1‑149 (MQ=255)                                                                  cgcgATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAATGCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAgcg                                                                                >  1:256486/1‑149 (MQ=255)                                                                              gTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccg                                                                    >  1:283542/1‑148 (MQ=255)                                                                              gTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccg                                                                    >  4:73615/1‑148 (MQ=255)                                                                               tGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccgc                                                                   >  6:222676/1‑148 (MQ=255)                                                                                     gCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTAt                                                            >  8:221352/1‑149 (MQ=255)                                                                                      ccGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                                           <  5:328225/149‑1 (MQ=255)                                                                                      ccGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                                           <  5:222676/149‑1 (MQ=255)                                                                                         cAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTc                                                        >  3:212789/1‑149 (MQ=255)                                                                                                   aGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgtg                                              <  2:289395/149‑1 (MQ=255)                                                                                                        ccAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAaca                                         <  5:39725/149‑1 (MQ=255)                                                                                                        ccAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAaca                                         >  3:243830/1‑149 (MQ=255)                                                                                                         cAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAacac                                        <  2:256486/149‑1 (MQ=255)                                                                                                           ggCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcacca                                      <  5:112473/149‑1 (MQ=255)                                                                                                            gCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcaccac                                     <  7:167887/149‑1 (MQ=255)                                                                                                            gCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAAcacca                                      >  2:313951/1‑148 (MQ=255)                                                                                                                 aaGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGGCGTTATTTGTGAAACACCACACTa                                 <  3:313531/148‑1 (MQ=255)                                                                                                                   gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              >  5:131002/1‑149 (MQ=255)                                                                                                                   gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              >  4:116368/1‑149 (MQ=255)                                                                                                                   gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              <  4:243830/149‑1 (MQ=255)                                                                                                                    tagcgcTGCAA‑GCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGTGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCTTAAGTTATGGTCGTTATTTGTGa                                             <  7:201918/129‑1 (MQ=21)                                                                                                                    tagcgcTGCAA‑GCAAAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGTGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTTTTTGTGa                                             >  8:201918/6‑134 (MQ=21)                                                                                                                       gACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc                            >  3:258920/1‑148 (MQ=255)                                                                                                                         cTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTAc                         <  3:39668/149‑1 (MQ=255)                                                                                                                            caa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                               <  3:201696/124‑1 (MQ=255)                                                                                                                            caa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTtgt                                               >  4:201696/1‑124 (MQ=255)                                                                                                                             aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc                            <  2:202745/142‑1 (MQ=255)                                                                                                                             aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc                            >  1:202745/1‑142 (MQ=255)                                                                                                                             aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCa                     <  4:258920/149‑1 (MQ=255)                                                                                                                             aa‑gcaaAAGCGAGAGCGATCCAGGTTCGCTCTGGCGTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGt                       >  4:255042/1‑147 (MQ=38)                                                                                                                                 caaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCTCAGCGCAAAAAACGCAGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCACCTACGTCAGCt                  >  4:180240/1‑149 (MQ=255)                                                                                                                                        cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              <  5:286956/130‑1 (MQ=255)                                                                                                                                        cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              >  6:286956/1‑130 (MQ=255)                                                                                                                                          gagCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGAt         >  5:177999/1‑149 (MQ=255)                                                                                                                                               aaCCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGaaa    <  1:150815/149‑1 (MQ=255)                                                                                                                                                 ccTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGaaaaa  <  1:59398/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                GCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAA  >  NC_002947/4831913‑4832203

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