Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	427,753	G→T	8.6%	L16M (CTG→ATG) ‡	PP_0350 ←	ferrichrome‑iron receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	427,753	0	G	T	8.6%	179.4 / 4.0	57	L16M (CTG→ATG) ‡	PP_0350	ferrichrome‑iron receptor
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base T (3/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGA  >  NC_002947/427607‑427900                                                                                                                                                    |                                                                                                                                                    gATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAgg                                                                                                                                                    <  4:218040/148‑1 (MQ=255) gATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGCGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAggg                                                                                                                                                   >  6:15005/1‑149 (MQ=255)      ccAAGGCAATGCCGGTCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCAt                                                                                                                                               >  7:141512/1‑148 (MQ=255)      ccAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATc                                                                                                                                              <  8:25204/149‑1 (MQ=255)        aaGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGCGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAg                                                                                                                                            >  3:180339/1‑149 (MQ=255)         aGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAgg                                                                                                                                           >  3:140825/1‑149 (MQ=255)            cAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTg                                                                                                                                        <  3:137168/149‑1 (MQ=255)              aTGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTTGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCTAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTggg                                                                                                                                      <  2:303882/149‑1 (MQ=255)                   ggCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCa                                                                                                                                 <  4:140825/149‑1 (MQ=255)                             gcGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGAt                                                                                                                        >  3:155436/1‑148 (MQ=255)                                     gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAATGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGc                                                                                                                 <  5:20497/147‑1 (MQ=255)                                     gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCaa                                                                                                                     <  6:55786/143‑1 (MQ=255)                                     gTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACgg                                                                                                                  <  1:300547/146‑1 (MQ=255)                                     gTTGAGGGTGGTGGCCAGGGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCaa                                                                                                                     <  8:60220/143‑1 (MQ=255)                                          gggTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGc                                                                                                                 <  1:177451/142‑1 (MQ=255)                                          gggTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGc                                                                                                                 >  2:177451/1‑142 (MQ=255)                                                 ggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGata                                                                                                                            <  7:307956/124‑1 (MQ=255)                                                 ggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGata                                                                                                                            >  8:307956/1‑124 (MQ=255)                                                 ggCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATgcggc                                                                                                   >  2:47302/1‑149 (MQ=255)                                                  gCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATgcggcg                                                                                                  >  3:273341/1‑149 (MQ=255)                                                  gCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATgcggcg                                                                                                  <  4:273341/149‑1 (MQ=255)                                                    ttttGGCGCGGTGGGCCTGTTGTTTGCGCGGGTAGTGGTTTTGCCGGCTGCGAACAAGGCACGGCTGGGCCTCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGcggcgg                                                                                                 <  4:262789/145‑1 (MQ=255)                                                       tGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTcc                                                                                             >  3:229405/1‑149 (MQ=255)                                                                 ggCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACggg                                                                                     <  5:40810/147‑1 (MQ=255)                                                                 ggCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGCCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGt                                                                                    >  5:142636/1‑148 (MQ=255)                                                                          tAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttgc                                                                          <  1:16574/149‑1 (MQ=255)                                                                          tAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGtgcttg                                                                           <  4:229405/148‑1 (MQ=255)                                                                                   gACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                                 <  1:47302/149‑1 (MQ=255)                                                                                   gACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                                 >  1:82036/1‑149 (MQ=255)                                                                                    aCGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTg                                                                <  8:141512/149‑1 (MQ=255)                                                                                        ggATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGat                                                                                                                             >  5:25770/1‑84 (MQ=255)                                                                                        ggATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGat                                                                                                                             <  6:25770/84‑1 (MQ=255)                                                                                          aTTGCTCGGCGGCGAACGAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGACACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagc                                                          >  5:122585/1‑149 (MQ=255)                                                                                           ttGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCATGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagcg                                                         <  5:177833/149‑1 (MQ=255)                                                                                            tGCTCGGCGGCGAACATGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCcgagcg                                                         <  7:86824/148‑1 (MQ=255)                                                                                                      cgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                              <  4:180339/149‑1 (MQ=255)                                                                                                      cgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                              <  5:15005/149‑1 (MQ=255)                                                                                                      cgAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                              <  2:82036/149‑1 (MQ=255)                                                                                                       gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGaaaa                                              <  1:120720/148‑1 (MQ=255)                                                                                                       gAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAg                                             >  6:300490/1‑149 (MQ=255)                                                                                                              gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCt                                                   >  1:237230/1‑136 (MQ=255)                                                                                                              gCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCt                                                   <  2:237230/136‑1 (MQ=255)                                                                                                                aGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                                 >  1:240059/1‑120 (MQ=255)                                                                                                                aGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCtt                                                                 <  2:240059/120‑1 (MQ=255)                                                                                                                 gggCTGGTCAGCGGGGCAGCCAGGCCGAGTGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGt                                    >  1:193135/1‑149 (MQ=255)                                                                                                                          aGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGaa                           >  1:10691/1‑148 (MQ=255)                                                                                                                                 cAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                    >  6:318080/1‑148 (MQ=255)                                                                                                                                    ccAGCCTGA‑TGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTgag                 >  3:66741/1‑148 (MQ=255)                                                                                                                                     cAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGaa               >  2:69717/1‑149 (MQ=255)                                                                                                                                        gCCGA‑TGGCCAGGGCCATCAGGCTGGGGGGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGt                                            <  3:250593/117‑1 (MQ=255)                                                                                                                                        gCCGA‑TGGCCAGGGCCATCAGGCTGGGGAGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGt                                            >  4:250593/1‑117 (MQ=255)                                                                                                                                        cctgA‑TGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                    >  6:327845/4‑141 (MQ=255)                                                                                                                                        cctgA‑TGGCCCTGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATAtt                    <  5:327845/138‑1 (MQ=255)                                                                                                                                         ctgA‑TGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattat           >  2:108819/3‑149 (MQ=255)                                                                                                                                           gA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGACATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTAttt         <  1:245810/149‑1 (MQ=255)                                                                                                                                            a‑TGGCCCTGGCCATCGGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAattatt          <  4:66741/147‑1 (MQ=255)                                                                                                                                            a‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTa        <  4:155436/149‑1 (MQ=255)                                                                                                                                                gCCAGGGCCATCAGGCGGGGGCGCAGGATATGATCACCGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAAt      <  8:282674/148‑1 (MQ=255)                                                                                                                                                   aGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGa  >  4:316561/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GATGGCCAAGGCAATGCCGGCCTGGCTGGCGATCTGGTTGAGGGTGGTGGCCAGTGGCGCGCTGGGCAGGTTGTAGGCGCGCGACGTGGATTGCTCGGCGGCGAACAAGGCAGGGCTGGTCAGCGGGGCAGCCAGGCCGA‑TGGCCAGGGCCATCAGGCTGGGGCGCAGGATATGATCAACGGCACGGGACATGCGGCGGTTCCTCGACGGGTAAGTGCTTGCTGCTTCCTTGCCGAGCGAGATCTGAAAAGTGACAGGGTGAATCTGAAAAATATTGAGAATTATTTAATGGGA  >  NC_002947/427607‑427900

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