Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,837,985	G→T	5.9%	R527M (AGG→ATG)	PP_3388 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,837,985	0	G	T	5.9%	221.2 / 3.0	68	R527M (AGG→ATG)	PP_3388	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base T (2/2);  total (32/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 = 8.01e-01

ACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACA  >  NC_002947/3837841‑3838131                                                                                                                                                 |                                                                                                                                                   aCGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGc                                                                                                                                                  >  7:94292/1‑147 (MQ=39) aCGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCaa                                                                                                                                                <  7:105701/149‑1 (MQ=39)    tGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGcc                                                                                                                                             <  7:267955/149‑1 (MQ=39)      gCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCa                                                                                                                                           >  4:304848/1‑149 (MQ=39)        aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAg                                                                                                                                                    <  1:339446/138‑1 (MQ=40)        aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAg                                                                                                                                                    >  2:339446/1‑138 (MQ=40)        aGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCACTTGCCAGGCAAGCCCATg                                                                                                                                         >  7:282825/1‑149 (MQ=255)         gACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGc                                                                                                                                        <  8:94292/149‑1 (MQ=40)          aCCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGcc                                                                                                                                       >  4:71127/1‑149 (MQ=40)           ccGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGcc                                                                                                                                       <  8:202519/148‑1 (MQ=40)                     cTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGc                                                                                                                            >  4:313185/1‑149 (MQ=255)                     cTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGc                                                                                                                            <  5:64423/149‑1 (MQ=255)                        gctaTACGAGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCtt                                                                                                                         <  5:212159/146‑1 (MQ=255)                                 ggCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCAt                                                                                                                <  2:41468/149‑1 (MQ=255)                                  gCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATc                                                                                                               <  5:173549/149‑1 (MQ=255)                                  gCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATc                                                                                                               <  8:175197/149‑1 (MQ=255)                                   cATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCg                                                                                                                                      >  7:82287/1‑125 (MQ=255)                                   cATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCg                                                                                                                                      <  8:82287/125‑1 (MQ=255)                                     tCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCACTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAac                                                                                                            <  8:282825/149‑1 (MQ=255)                                       aCCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAaca                                                                                                           <  2:317089/148‑1 (MQ=255)                                           aCAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTg                                                                                                      >  2:109078/1‑149 (MQ=255)                                            cAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCt                                                                                                       <  3:54760/147‑1 (MQ=255)                                             aaCCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTa                                                                                                    >  6:23838/1‑149 (MQ=255)                                               ccAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTAc                                                                                                   >  6:216525/1‑148 (MQ=255)                                                 aggaggACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTTCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGaa                                                                                                <  3:71127/149‑1 (MQ=255)                                                   gaggaCCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCATGCAAGCCCATGCCGCTGCCAATGCCTTTGGgc                                                                                                                    <  2:175579/127‑1 (MQ=255)                                                   gaggaCCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCATGCAAGCCCATGCCGCTGCCAATGCCTTTGGgc                                                                                                                    >  1:175579/1‑127 (MQ=255)                                                   gaggaCCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACa                                                                                              <  2:175290/149‑1 (MQ=255)                                                       aCCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATg                                                                                          <  3:304848/149‑1 (MQ=255)                                                         cTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATg                                                                                          >  7:105992/1‑147 (MQ=255)                                                            cGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAg                                                                                       <  1:60264/147‑1 (MQ=255)                                                            cGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAg                                                                                       <  1:55380/147‑1 (MQ=255)                                                                 agaaCTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGcaca                                                                                  >  3:73148/1‑147 (MQ=255)                                                                        gcggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCATGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATg                                                                                          >  7:176462/1‑132 (MQ=255)                                                                        gcggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCATGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATg                                                                                          <  8:176462/132‑1 (MQ=255)                                                                        gcggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          <  3:313185/148‑1 (MQ=255)                                                                        gcggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAgcagca                                                                          <  5:23838/148‑1 (MQ=255)                                                                          ggcgCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACgg                                                                       >  5:149114/1‑149 (MQ=255)                                                                             gctgctGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGc                                                                      >  1:16331/1‑147 (MQ=255)                                                                                 ctgGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCAt                                                                >  2:49273/1‑149 (MQ=255)                                                                                   ggAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCg                                                              >  7:119166/1‑149 (MQ=255)                                                                                   ggAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCg                                                              >  6:72951/1‑149 (MQ=255)                                                                                    gAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCg                                                              <  3:243656/148‑1 (MQ=255)                                                                                     aaGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGc                                                            <  5:216525/149‑1 (MQ=255)                                                                                       gccgccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGa                                                                                            >  6:211845/1‑115 (MQ=255)                                                                                       gccgccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGa                                                                                            <  5:211845/115‑1 (MQ=255)                                                                                          gccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCt                                                          >  3:248044/1‑146 (MQ=255)                                                                                          gccAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCt                                                          <  4:248044/146‑1 (MQ=255)                                                                                               ggCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCt                                                  <  8:105992/149‑1 (MQ=39)                                                                                               ggCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCt                                                  <  7:217681/149‑1 (MQ=39)                                                                                                 ccGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTgg                                                <  4:73148/149‑1 (MQ=38)                                                                                                  cGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTgg                                                >  5:152793/1‑148 (MQ=37)                                                                                                  cGAACCGCCGTCAGCGTCATCCTTCGCCATGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTgg                                                >  1:19230/1‑148 (MQ=37)                                                                                                    aaCCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTAc                                             >  1:174530/1‑149 (MQ=37)                                                                                                       cgccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTc                                          <  8:119166/149‑1 (MQ=37)                                                                                                                cGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTgggg                                 <  2:16331/149‑1 (MQ=39)                                                                                                                 gTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGc                                                                      <  6:23037/111‑1 (MQ=39)                                                                                                                 gTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGc                                                                      >  5:23037/1‑111 (MQ=39)                                                                                                                 gTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGggc                                <  5:72951/149‑1 (MQ=39)                                                                                                                     tCCTTAGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATTCCTTTGGGCGCATCAACACACTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGggcggc                             >  4:33623/1‑148 (MQ=255)                                                                                                                        ttCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGt                         >  4:279576/1‑149 (MQ=39)                                                                                                                        ttCGCCCTGCAGGAAGCGCTTGCCAGGCAAGCCCCTGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGt                         >  3:183057/1‑149 (MQ=255)                                                                                                                           gCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGtgct                      >  4:308477/1‑149 (MQ=39)                                                                                                                            cccTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGtgctg                     >  8:184699/1‑149 (MQ=39)                                                                                                                             ccTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGtgctgc                    <  6:152793/149‑1 (MQ=39)                                                                                                                                    gAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGc             <  2:174530/149‑1 (MQ=39)                                                                                                                                        cgcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCAcc         >  3:296776/1‑149 (MQ=255)                                                                                                                                               cAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGGTGCGCACCGCCACGCCCAACa  <  4:183057/149‑1 (MQ=255)                                                                                                                                               cAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACa  >  8:215107/1‑149 (MQ=255)                                                                                                                                               cAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACACAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACa  >  8:33317/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACGTGGCAGACCGGGGCAACCTGGCAATACGCGGCATCACCTACAACCAGGAGGACCTGCGGGAAGAACTGGCGGCGCTGCTGGAAGCCGCCAAGGCCGAACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCAACA  >  NC_002947/3837841‑3838131

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