Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,244,133	T→A	6.3%	M117L (ATG→TTG)	hmgA ←	homogentisate 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,244,133	0	T	A	6.3%	212.4 / 3.0	74	M117L (ATG→TTG)	hmgA	homogentisate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base A (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

GCTCGGTGGCGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGC  >  NC_002947/5243986‑5244279                                                                                                                                                    |                                                                                                                                                    gCTCGGTGGCGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATg                                                                                                                                                    >  6:47751/1‑149 (MQ=255)     ggTGGCGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGc                                                                                                                                                 <  5:47751/148‑1 (MQ=255)          cGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCTGCCCATGGCTGCCaa                                                                                                                                           <  3:24538/149‑1 (MQ=255)          cGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCTGCCCATGGCTGCCaa                                                                                                                                           <  4:144368/149‑1 (MQ=255)           gATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAAc                                                                                                                                          >  8:26149/1‑149 (MQ=255)           gATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAAc                                                                                                                                          >  2:326719/1‑149 (MQ=255)           gATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAAc                                                                                                                                          >  6:372506/1‑149 (MQ=255)             tACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCt                                                                                                                                        >  8:52290/1‑149 (MQ=255)              aCGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCtt                                                                                                                                       <  7:282436/149‑1 (MQ=255)                 cAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGa                                                                                                                                    <  1:207958/149‑1 (MQ=255)                  aGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAATAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGAt                                                                                                                                   <  7:150592/149‑1 (MQ=255)                  aGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGAt                                                                                                                                   <  6:269930/149‑1 (MQ=255)                         ccTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCATGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAAtcg                                                                                                                            >  4:280685/1‑149 (MQ=255)                            gTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAAtcggtc                                                                                                                         >  1:325542/1‑149 (MQ=255)                               ccGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAAtcggtcgg                                                                                                                       <  6:255495/148‑1 (MQ=255)                                  gCACCAGTATCAGCTCACCGGATGCGTTGAAGAACACCCGTTCCATGTACCTGTTGGCGCTGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGGGTTGGCAGCCATGGGCAGCCAACCTTCGATTAAATCGGTCTGTTca                                                                                                                   <  4:182892/149‑1 (MQ=255)                                        gTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAgg                                                                                                              >  2:256485/1‑148 (MQ=255)                                           gcagcTCACAGTCTGAGTTGAAGAACACCCGTTCCATGGACAGGTTGGCGCGGTAGATGTAGATGCTCAAGCCAGCCGGTTTTTCCGCTCAGGCGTTGGCATCAATGGGCAGCAAACCTTCGATTAAATCGGTCTGTTAAGCAGGaa                                                                                                            >  1:177988/1‑147 (MQ=255)                                                  aCCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCg                                                                                                  >  7:198859/1‑149 (MQ=255)                                                       atGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTc                                                                                             >  1:207037/2‑149 (MQ=255)                                                           gTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAgc                                                                                         >  2:102229/1‑149 (MQ=255)                                                                agaaCACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAg                                                                                     <  6:312555/148‑1 (MQ=255)                                                                 gaaCACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCATGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGc                                                                                   <  3:280685/149‑1 (MQ=255)                                                                 gaaCACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGc                                                                                   >  3:91363/1‑149 (MQ=255)                                                                 gaaCACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGc                                                                                   <  2:207037/149‑1 (MQ=255)                                                                  aaCACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCg                                                                                  >  5:106981/1‑149 (MQ=255)                                                                    cacCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGt                                                                                >  1:138133/1‑149 (MQ=255)                                                                      cccGTTCCATGGACCGGTTGGCGCTGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTg                                                                              <  8:295533/149‑1 (MQ=255)                                                                        cGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTggg                                                                            >  2:354295/1‑149 (MQ=255)                                                                        cGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTggg                                                                            >  6:344210/1‑149 (MQ=255)                                                                         gTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTgggg                                                                           <  8:44124/149‑1 (MQ=255)                                                                             cATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGa                                                                       >  6:203875/1‑149 (MQ=255)                                                                              aTGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGAt                                                                      >  7:272494/1‑149 (MQ=255)                                                                              aTGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGAt                                                                      >  6:85417/1‑149 (MQ=255)                                                                              aTGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGAt                                                                      <  4:23969/149‑1 (MQ=255)                                                                              aTGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGAt                                                                      <  2:372545/149‑1 (MQ=255)                                                                                ggACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGc                                                                    >  8:152635/1‑149 (MQ=255)                                                                                  aCCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCa                                                                  <  1:300151/149‑1 (MQ=255)                                                                                    cGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCAcc                                                                >  3:86443/1‑149 (MQ=255)                                                                                      gTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCa                                                              >  3:9567/1‑149 (MQ=255)                                                                                       ttGGCGCGGTAGATGTAGATGCTCATGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTACAGCGCAGGCGGTTGTGGTTGATGCCACCCAg                                                             >  6:260355/1‑149 (MQ=255)                                                                                            gcgGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGcc                                                        >  6:247832/1‑149 (MQ=255)                                                                                               gtagatgtagatGCTCACGCCTGCCGGTTTTTCCGCGCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCccgc                                                     <  7:253803/149‑1 (MQ=255)                                                                                                  gatgtagatgCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCccgcc                                                    >  4:73762/1‑147 (MQ=255)                                                                                                        gatgCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCt                                            >  7:272406/1‑149 (MQ=255)                                                                                                         atgCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTGCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTg                                           <  5:344210/149‑1 (MQ=255)                                                                                                           gCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTgg                                          <  7:52290/148‑1 (MQ=255)                                                                                                             tCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCgg                                                                                 >  6:47514/1‑107 (MQ=255)                                                                                                             tCACGACAGCCGGTTTTACCGCTCAGGCGTTGTAAGCCAAGGCCACCCAACCTACCATGAAAAAGGACGGTTCAGCACCAAT‑CGGCTGCGGGCTCCAGCGCAGGCCGTTAGGGTTGATCACACCCATTGGCACGCCGCGCGGCTGgcgc                                       >  8:122533/1‑149 (MQ=255)                                                                                                               gcGCCAACCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGgcgcgc                                     <  5:47009/148‑1 (MQ=255)                                                                                                               aCGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGgcgcgc                                     >  3:204767/1‑149 (MQ=255)                                                                                                               aCGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGgcgcgc                                     <  8:272494/149‑1 (MQ=255)                                                                                                                cGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGgcgcgc                                     >  5:250281/1‑148 (MQ=255)                                                                                                                 gCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGCTCAGCAGTAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCa                                   <  5:52740/149‑1 (MQ=255)                                                                                                                   catgctgTTTTTCCGATCATGCGTTTGAATCCATGGGCAGCCAACCTTCGATGAAATCGGTAGGTTCAGAATGAAT‑CGTCTGCGGTCTCCATCTCAGGAGGTTGTGGTTGATGACACCCAGTGGCCCGCCGAGCTGCTTGCGCGCCAgg                                 <  2:177988/143‑1 (MQ=255)                                                                                                                     gCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGgcg                               <  2:325542/149‑1 (MQ=255)                                                                                                                       cGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCt                             <  5:85417/149‑1 (MQ=255)                                                                                                                       cGGTTTTTCCGCGCCGGCGTTGGCAGCCCTGGGCAGCCAACCTTCGATGAAATCGGTCGGTGCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCt                             <  6:288935/149‑1 (MQ=255)                                                                                                                          tttttCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCgg                          <  7:152635/149‑2 (MQ=255)                                                                                                                            tttCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                         <  5:247832/148‑1 (MQ=255)                                                                                                                            tttCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGaa                         <  6:106981/148‑1 (MQ=255)                                                                                                                                     cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGgc                                >  4:358619/1‑132 (MQ=255)                                                                                                                                     cGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGt                 <  8:272406/147‑1 (MQ=255)                                                                                                                                       gCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGc               >  5:310000/1‑147 (MQ=255)                                                                                                                                        ggTTGGCTGCCCATGGCTGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCg                           >  6:154408/2‑134 (MQ=255)                                                                                                                                        ggTTGGCTGCCCATGGCTGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCg                           <  5:154408/133‑1 (MQ=255)                                                                                                                                         gTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAgg           <  5:203875/149‑1 (MQ=255)                                                                                                                                             gCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCGGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAg        >  6:262544/1‑148 (MQ=255)                                                                                                                                             gCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAg        >  6:87970/1‑148 (MQ=255)                                                                                                                                              cAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGaa      <  7:26149/149‑1 (MQ=255)                                                                                                                                              cAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGaa      <  4:86443/149‑1 (MQ=255)                                                                                                                                                gCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAgg    <  1:326719/149‑1 (MQ=255)                                                                                                                                                   aTGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAATCCGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAACGGc  >  1:265007/1‑149 (MQ=255)                                                                                                                                                   aTGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGc  >  1:187852/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCTCGGTGGCGATACGCAGGCGGCCCTGTTCCGGCACCAGTAGCAGCTCACCGTCTGCGTTGAAGAACACCCGTTCCATGGACCGGTTGGCGCGGTAGATGTAGATGCTCACGCCAGCCGGTTTTTCCGCTCCGGCGTTGGCAGCCATGGGCAGCCAACCTTCGATGAAATCGGTCGGTTCAGCAGGAAT‑CGGCTGCGGGCTCCAGCGCAGGCGGTTGGGGTTGATGCCACCCAGTGGCCCGCCGAGCGGCTGGCGCGCCAGGCGCTCGAAGCGTGGGTGCAAGGCAGAAGGGC  >  NC_002947/5243986‑5244279

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