Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,832,057	T→A	6.7%	intergenic (+32/+120)	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,057	0	T	A	6.7%	181.5 / 2.0	59	intergenic (+32/+120)	pvdS/pvdY	alternative sigma factor/hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAA  >  NC_002947/4831910‑4832203                                                                                                                                                      |                                                                                                                                                   aTGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCt                                                                                                                                                    <  7:34578/148‑1 (MQ=255)    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                                                                                                                         <  3:45680/140‑1 (MQ=255)                                    ggCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAgcg                                                                                                                         >  4:45680/1‑140 (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                                                                                                        >  8:89785/1‑144 (MQ=255)                                                 ccTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGc                                                                                                        <  7:89785/144‑1 (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                                                                    >  4:73615/1‑148 (MQ=255)                                                                                 gTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAAcgccg                                                                    >  1:283542/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:222676/149‑1 (MQ=255)                                                                                         ccGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATg                                                           <  5:328225/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                                         >  3:243830/1‑149 (MQ=255)                                                                                                           ccAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAaca                                         <  5:39725/149‑1 (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                              >  4:116368/1‑149 (MQ=255)                                                                                                                      gAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              >  5:131002/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                            >  1:202745/1‑142 (MQ=255)                                                                                                                                aa‑gcaaACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcccc                            <  2:202745/142‑1 (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                              >  6:286956/1‑130 (MQ=255)                                                                                                                                           cAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGcc                              <  5:286956/130‑1 (MQ=255)                                                                                                                                             gagCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGAt         >  5:177999/1‑149 (MQ=255)                                                                                                                                                  aaCCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGaaa    <  1:150815/149‑1 (MQ=255)                                                                                                                                                    ccTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGaaaaa  <  1:59398/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ATGGCGTGCCGCAGAAGGACATTGCCAAGGAACTGGGCGTATCGCCGACCTTGGTCAACTTCATGATCCGCGATGCGCTGGTGCATTGCCGCAAGACGGCCAGTCGCCAGGCCTGAAGAT‑GACTGCAA‑GCAAACGCCAGAGCGAACCTGGATCGCTCTCGCTTTTGCTTGCAGCGCTAGTATGGTGTTTTGCAAATACGGCCAATAACTCGCAGCGCAAAAAACGCCGCAAGTTATGGTCGTTATTTGTGAAACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAA  >  NC_002947/4831910‑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: ‑