Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	852,159	A→T	5.1%	Q313L (CAG→CTG)	hemA →	glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	852,159	0	A	T	5.1%	235.0 / 2.4	78	Q313L (CAG→CTG)	hemA	glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base T (2/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGC  >  NC_002947/852027‑852303                                                                                                                                     |                                                                                                                                                  acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCgg                                                                                                                                           >  5:143221/1‑141 (MQ=255) acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCgg                                                                                                                                           <  6:143221/141‑1 (MQ=255) acaAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAggcg                                                                                                                                   >  3:196799/1‑149 (MQ=255)      cccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCg                                                                                                                               >  6:321294/1‑148 (MQ=255)       ccATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGaa                                                                                                                             <  7:137276/149‑1 (MQ=255)       ccATCTTCATGATCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCg                                                                                                                               >  4:7950/1‑147 (MQ=255)         aTCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCTGGGCCGCGCCCTggcggctcttga                                                                                                                           <  3:22913/149‑7 (MQ=255)            ttCATGGTCGACATCGCTGTACCTCGCGATAGCGCAACCGAAGTCGGCGAGCTGGATGATGTCTCCCTGTACACCGTCGACGACCTGCACGATGGAGGGGCGGAACACCTCAAGCGCCGCCAGGGCGCGGCCCCGGCGGCCGCAGCGCt                                                                                                                        >  6:245900/1‑149 (MQ=255)              cATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGc                                                                                                                         <  7:175968/146‑1 (MQ=255)               aTGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCCCGATGTAGGGGCGGAAACCCTCACGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTgg                                                                                                                      >  6:237168/1‑148 (MQ=255)                     gACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGt                                                                                                               >  7:241161/1‑149 (MQ=255)                      aCATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg                                                                                                              >  3:214327/1‑149 (MQ=255)                       cATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg‑g                                                                                                            <  4:214327/149‑1 (MQ=255)                               tACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGaa                                                                                                    <  8:241161/149‑1 (MQ=255)                                aCCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCg                                                                                                      >  2:109470/1‑146 (MQ=255)                                       gATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg‑g                                                                                                            >  8:58501/1‑133 (MQ=255)                                       gATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTg‑g                                                                                                            <  7:58501/133‑1 (MQ=255)                                         tatCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCAtgc                                                                                          >  7:146719/1‑149 (MQ=255)                                            cGAAACCGAACTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACCACCTCCACGATGAAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCAAGGCGGCCGAACAGCTGCTATCCGTG‑GGCGCCGAACACTTCAtccacc                                                                                       >  8:200451/1‑144 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                      >  4:301412/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                      >  2:189252/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                      >  6:43315/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                      >  2:157612/1‑149 (MQ=255)                                               aaCCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcgcc                                                                                    >  1:80906/1‑149 (MQ=255)                                                    aaGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcgcctgcgc                                                                               <  5:321294/149‑1 (MQ=255)                                                     aGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTGCGGCGCCGAAGACTTCATGctgcgcctgcgc                                                                               <  1:189252/149‑1 (MQ=255)                                                         ggCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACt                                                                          >  6:46663/1‑149 (MQ=255)                                                          gCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTg                                                                         <  1:162248/149‑1 (MQ=255)                                                              gCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTg                                                                     >  7:80228/1‑149 (MQ=255)                                                                  gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                                 <  5:43315/149‑1 (MQ=255)                                                                  gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                                 <  3:7950/149‑1 (MQ=255)                                                                      atgTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                             <  5:82988/149‑1 (MQ=255)                                                                      atgTCTACCTGTAAACCGTAGACGACCTGCACGATGTAGTGTCGGAAAACCTAAAGAGCCGCCATGGCGCGGCCCAGGAGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCATGCGCGAAc                                                                           >  2:293590/1‑135 (MQ=255)                                                                      atgTCTACATGTACACCGTCGACGCCCTGCACTATGTAGTGGCGTAAAACCTCAAGAGCCGCCATGTCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAATACTTCATGCTGCGCCTGCGCTAAc                                                                           <  1:293590/135‑1 (MQ=255)                                                                         tCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGt                                                           <  7:85495/148‑1 (MQ=255)                                                                          cTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTgg                                                               >  7:60576/1‑143 (MQ=255)                                                                          cTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTgg                                                               <  8:60576/143‑1 (MQ=255)                                                                            aCCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCt                                                        >  2:159228/1‑148 (MQ=255)                                                                              cTGTACCCCGTCGACGCCCTGCCCGAGGGAGTGGCGGAAACCCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGGATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCGGCAGTGGACGTGCTCAg                                                     <  5:245900/149‑1 (MQ=255)                                                                               tGTACTCGGTCGACGACCTGCACGTTGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAgg                                                    <  8:146719/149‑1 (MQ=255)                                                                                gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTAt                                                                                                                   <  3:166810/86‑1 (MQ=255)                                                                                gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTAt                                                                                                                   >  4:166810/1‑86 (MQ=255)                                                                                gTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGAGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                   >  7:244402/1‑149 (MQ=255)                                                                                 tACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAgg                                                    <  7:200451/147‑1 (MQ=255)                                                                                 tACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGc                                                  <  8:80228/149‑1 (MQ=255)                                                                                   cacCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGcc                                                 >  5:289225/1‑148 (MQ=255)                                                                                    acCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTa                                               >  2:52340/1‑149 (MQ=255)                                                                                          gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGc                                           >  1:140205/1‑147 (MQ=255)                                                                                          gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGc                                           >  3:65116/1‑147 (MQ=255)                                                                                          gacgacCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCa                                         <  3:301412/149‑1 (MQ=255)                                                                                               ccTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCa                                      >  7:234915/1‑147 (MQ=255)                                                                                                cTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGCCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAg                                   >  3:147318/1‑149 (MQ=255)                                                                                                   cACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGt                                  >  3:286541/1‑147 (MQ=255)                                                                                                     cGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgc                              <  1:52340/149‑1 (MQ=255)                                                                                                      gATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGTGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgcg                             <  8:244402/149‑1 (MQ=255)                                                                                                      gATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAgcg                             >  1:95320/1‑149 (MQ=255)                                                                                                            gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAg                                                     <  3:240230/119‑1 (MQ=255)                                                                                                            gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAg                                                     >  4:240230/1‑119 (MQ=255)                                                                                                               gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                    >  1:218791/1‑149 (MQ=255)                                                                                                               gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCACGCTGCGCCTGCGCGAACTGGCTGCAGTGGAGGTGCTCATGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                    >  5:222252/1‑149 (MQ=255)                                                                                                               gCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCACGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCATGGCCTATCGCCAGCAAAGTGAGCGCCTTCGcga                    <  6:222252/149‑1 (MQ=255)                                                                                                                    aaaCCTCAAGAGCCGCCAGGGCGCGGCCCATGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACgaag               <  6:61975/149‑1 (MQ=255)                                                                                                                       ccTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAAc            <  4:147318/149‑1 (MQ=255)                                                                                                                        cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACt           >  2:316172/1‑149 (MQ=255)                                                                                                                        cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACt           >  1:312679/1‑149 (MQ=255)                                                                                                                        cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCGGGGATACGTG‑GGCGCCGACGACTTCATGCTGCGCCTGCGCGTACTGGATGCTGTGGCAGTGCGCAGGGCCTAGAGCGCGCGAAGTGAGCTCCTGCGCGACGACGAGCt           >  6:188327/1‑149 (MQ=255)                                                                                                                           aaGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTg          <  1:170660/147‑1 (MQ=255)                                                                                                                            agagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaa       <  8:234915/149‑1 (MQ=255)                                                                                                                             gagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                             <  8:275548/94‑1 (MQ=255)                                                                                                                             gagCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                             >  7:275548/1‑94 (MQ=255)                                                                                                                              agCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaaaa     <  2:95320/149‑1 (MQ=255)                                                                                                                              agCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCaaaa     <  2:312679/149‑1 (MQ=255)                                                                                                                                ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                   >  2:181103/1‑101 (MQ=255)                                                                                                                                ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAggg                                                   <  1:181103/101‑1 (MQ=255)                                                                                                                                ccgccTGGGCCGCGCCCTGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAgg   >  7:272782/1‑149 (MQ=255)                                                                                                                                 cgccAGGGCGCGGCGCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGc  <  4:196799/149‑1 (MQ=255)                                                                                                                                   ccAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCCCCGAACACTTCATGCTGCGCCTGCCCGAACTCGCTGCAGTGGAc                                                             >  6:222980/1‑88 (MQ=255)                                                                                                                                   acAGGGCGCGGCCCTGGCGGCCTAATAGCTGGTATCGGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCTGTGGAc                                                             <  5:222980/87‑1 (MQ=255)                                                                                                                                     |                                                                                                                                                  ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGC  >  NC_002947/852027‑852303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑