Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	852,171	A→T	5.6%	Q317L (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,171	0	A	T	5.6%	246.7 / 2.0	83	Q317L (CAG→CTG)	hemA	glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base T (2/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCCTGCTG  >  NC_002947/852027‑852316                                                                                                                                                 |                                                                                                                                                   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                                                                                                   >  2:157612/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                                   >  2:189252/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                                   >  6:43315/1‑149 (MQ=255)                                             gAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGctgcg                                                                                                   >  4:301412/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                                                                              <  3:7950/149‑1 (MQ=255)                                                                  gatgatGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGt                                                                              <  5:43315/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                                                                            <  8:60576/143‑1 (MQ=255)                                                                          cTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTgg                                                                            >  7:60576/1‑143 (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                                                                  >  4:240230/1‑119 (MQ=255)                                                                                                            gTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAg                                                                  <  3:240230/119‑1 (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                        >  1:312679/1‑149 (MQ=255)                                                                                                                        cTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACt                        >  2:316172/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)                                                                                                                                            ggCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgc    <  4:65116/149‑1 (MQ=255)                                                                                                                                             gCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                                          >  1:51072/1‑78 (MQ=255)                                                                                                                                             gCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGAc                                                                          <  2:51072/78‑1 (MQ=255)                                                                                                                                             gCCCAGGCGGCCGAAGAGCTCGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgct   <  2:80906/149‑1 (MQ=255)                                                                                                                                              cccAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgctg  <  5:46663/149‑1 (MQ=255)                                                                                                                                              cccAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgctg  <  1:159228/149‑1 (MQ=255)                                                                                                                                              cccAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCAAACTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCctgctg  >  1:226051/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACAAGCCCATCTTCATGGTCGACATCGCTGTACCTCGCGATATCGAAACCGAAGTCGGCGAGCTGGATGATGTCTACCTGTACACCGTCGACGACCTGCACGATGTAGTGGCGGAAAACCTCAAGAGCCGCCAGGGCGCGGCCCAGGCGGCCGAAGAGCTGGTATCCGTG‑GGCGCCGAAGACTTCATGCTGCGCCTGCGCGAACTGGCTGCAGTGGACGTGCTCAGGGCCTATCGCCAGCAAAGTGAGCGCCTTCGCGACGAAGAACTGCAAAAGGCCTTGCGCCTGCTG  >  NC_002947/852027‑852316

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