Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	427,742	G→C	6.2%	G19G (GGC→GGG)	PP_0350 ←	ferrichrome‑iron receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	427,742	0	G	C	6.2%	204.6 / 3.4	64	G19G (GGC→GGG)	PP_0350	ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base C (2/2);  total (32/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 = 8.45e-01

CAGTGCCGGGTCGATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGA  >  NC_002947/427595‑427886                                                                                                                                                     |                                                                                                                                                 cggtaccGGGTCGATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGGTGAGGGGGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGc                                                                                                                                                 <  2:330487/144‑1 (MQ=255) cAGTGCCGGGTCGATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGc                                                                                                                                                 <  3:181649/149‑1 (MQ=255)             gATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAggg                                                                                                                                     >  5:202280/1‑149 (MQ=255)                  ccAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCT‑GGCAGGTTGTACGCGGCGCG‑CGTGGATTGCGGGGCGGCGGACCAGGCAGGGCAGGTCAGGGGGGCAGCCAGGCCCATGGCCAGGGccctgc                                                                                                                               >  1:52539/1‑145 (MQ=255)                   cAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCa                                                                                                                               <  3:289369/149‑1 (MQ=255)                    aaGGCAATGCCGGCCTGGCTGGCGATCTGGTTGA‑GGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAgg                                                                                                                             >  7:43269/1‑149 (MQ=255)                        cAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGGGGGGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTg                                                                                                                          <  1:240206/149‑1 (MQ=255)                          aTGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGa                                                                                                                                              >  8:207342/1‑127 (MQ=255)                          aTGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGa                                                                                                                                              <  7:207342/127‑1 (MQ=255)                           tGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTgggg                                                                                                                       <  6:202280/149‑1 (MQ=255)                              cGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGgcgc                                                                                                                    <  2:298158/149‑1 (MQ=255)                               ggCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCa                                                                                                                   <  2:66572/149‑1 (MQ=255)                                   tggctggcGATCTGGTTGAGGGGGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGa                                                                                                                >  7:227237/1‑148 (MQ=255)                                     gctggcGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGTatat                                                                                                             <  2:95370/149‑1 (MQ=255)                                     gctggcGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGatat                                                                                                             <  1:284049/149‑1 (MQ=255)                                       tggcGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGa                                                                                                           >  5:216468/1‑149 (MQ=255)                                               tGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACgg                                                                                                    >  3:178609/1‑148 (MQ=255)                                                 gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTggg                                                                                                                        <  2:6614/126‑1 (MQ=255)                                                 gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTggg                                                                                                                        >  1:6614/1‑126 (MQ=255)                                                 gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCAc                                                                                                 <  4:178609/149‑1 (MQ=255)                                                     ctggtgggGGCCAGTGGCGCGCTGTGCATGTTGTCGGC‑GCGCGACGTGGCTTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCTGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATAAACGGCACggg                                                                                              <  8:75185/146‑1 (MQ=255)                                                                cAGTTGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGt                                                                                  >  5:251202/1‑149 (MQ=255)                                                                cAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAAc                                                                                                      >  3:135342/1‑129 (MQ=255)                                                                cAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAAc                                                                                                      <  4:135342/129‑1 (MQ=255)                                                                cAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGt                                                                                  >  4:272041/1‑149 (MQ=255)                                                                  gTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCTACGGCACGGGACATGCGGCGGTTc                                                                                <  8:227237/149‑1 (MQ=255)                                                                     gcgcgcTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCt                                                                                                                           <  5:27833/103‑1 (MQ=255)                                                                     gcgcgcTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCt                                                                                                                           >  6:27833/1‑103 (MQ=255)                                                                     gcgcgcTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTc                                                                             <  6:216468/149‑1 (MQ=255)                                                                         gcTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACg                                                                         >  4:275149/1‑149 (MQ=255)                                                                             ggCAGTTTGTAGGC‑GCGCTACGTTGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGt                                                                      <  6:251202/148‑1 (MQ=255)                                                                              gCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGt                                                                      >  3:143018/1‑147 (MQ=255)                                                                                  gTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTc                                                                             >  1:25443/1‑136 (MQ=255)                                                                                  gTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTc                                                                             <  2:25443/136‑1 (MQ=255)                                                                                     gTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttg                                                             >  3:159516/1‑149 (MQ=255)                                                                                             cgcgACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTc                                                      <  4:159516/149‑1 (MQ=255)                                                                                                       aTTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                            >  1:253751/1‑149 (MQ=255)                                                                                                       aTTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                            <  3:272041/149‑1 (MQ=255)                                                                                                           cTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                            >  4:158332/1‑145 (MQ=255)                                                                                                            tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTtgctg                                                          >  7:135952/1‑130 (MQ=255)                                                                                                            tCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTtgctg                                                          <  8:135952/130‑1 (MQ=255)                                                                                                                 ggcgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGa                                   >  5:46188/1‑148 (MQ=255)                                                                                                                    gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTaa                                                                    >  8:49658/1‑112 (MQ=255)                                                                                                                    gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGAAATGCGGCGGTTCCTCGACGGGTaa                                                                    <  7:49658/112‑1 (MQ=255)                                                                                                                     aaCAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGt                              <  2:202044/149‑1 (MQ=255)                                                                                                                       cAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                            <  3:62877/133‑1 (MQ=255)                                                                                                                       cAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                            >  4:62877/1‑133 (MQ=255)                                                                                                                       cAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGa                            >  3:259790/1‑149 (MQ=255)                                                                                                                       cAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGa                            >  7:6685/1‑149 (MQ=255)                                                                                                                        aaGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGa                            >  3:94341/1‑148 (MQ=255)                                                                                                                          ggCAGGGCTGGTCAGCGGGGCAGCAAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCATCGACGGGTAAGTGCTTGATGCTTCCTTGCCGAGCGAGATATGAAAAGTTACAg                         >  6:368067/1‑149 (MQ=255)                                                                                                                           gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAgg                        >  2:327699/1‑149 (MQ=255)                                                                                                                             aGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGt                      >  1:7773/1‑149 (MQ=255)                                                                                                                               ggCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGa                    <  4:94341/149‑1 (MQ=255)                                                                                                                                 cTGGTCAGCGGGACAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCgaga                                        <  1:43699/127‑1 (MQ=255)                                                                                                                                 cTGGTCAGCGGGACAGCCAGGCCGATGGCCAGGGCCATCAGGCTCGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCgaga                                        >  2:43699/1‑127 (MQ=255)                                                                                                                                    gTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATc                 <  3:158332/147‑1 (MQ=255)                                                                                                                                    gTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCt                >  3:107674/1‑148 (MQ=255)                                                                                                                                       aGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGaaa            <  3:275149/149‑1 (MQ=255)                                                                                                                                               aGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCTGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTga    <  5:368067/149‑1 (MQ=255)                                                                                                                                                cccAGCCTGATGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTCCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt      >  2:274320/2‑146 (MQ=255)                                                                                                                                                cccAGCCTGATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                <  4:160903/119‑1 (MQ=255)                                                                                                                                                cccAGCCTGATGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                >  3:160903/2‑120 (MQ=255)                                                                                                                                                 ccAGCCTGATGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTCCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTgaga  <  1:274320/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CAGTGCCGGGTCGATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGC‑GCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGATGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGA  >  NC_002947/427595‑427886

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