Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,523,564	G→A	8.2%	R437C (CGC→TGC)	yhjG ←	protein YhjG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,523,564	0	G	A	8.2%	176.6 / 4.5	61	R437C (CGC→TGC)	yhjG	protein YhjG
Reads supporting (aligned to +/- strand):  ref base G (32/24);  new base A (3/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATTGGCGCGAAGGTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGGCAGCACCTTG  >  NC_002947/5523416‑5523712                                                                                                                                                     |                                                                                                                                                     aTTGGCGCGAAGGTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCg                                                                                                                                                      >  5:7172/1‑149 (MQ=255)     gcgcGAAGGTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTaa                                                                                                                                                   >  6:735/1‑148 (MQ=255)             gTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAgg                                                                                                                                               <  8:161950/144‑1 (MQ=255)             gTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAgg                                                                                                                                               >  7:161950/1‑144 (MQ=255)              tAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCt                                                                                                                                         <  5:90314/149‑1 (MQ=255)                  gAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTc                                                                                                                                       >  3:183325/1‑147 (MQ=255)                  gAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGa                                                                                                                                     <  6:288633/149‑1 (MQ=255)                        ttGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGcacca                                                                                                                               <  5:735/149‑1 (MQ=255)                           cttgagcttgaAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACAt                                                                                                                            >  7:316510/1‑149 (MQ=255)                            ttgagcttgaAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATg                                                                                                                           >  3:75928/1‑149 (MQ=255)                                  ttgaAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGa                                                                                                                     >  1:119412/1‑149 (MQ=255)                                       gCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGt                                                                                                                <  6:344995/149‑1 (MQ=255)                                         cACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCa                                                                                                              <  6:7172/149‑1 (MQ=255)                                         cACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCa                                                                                                              <  4:75928/149‑1 (MQ=255)                                                    aaCTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAg                                                                                                   >  2:89536/1‑149 (MQ=255)                                                    aaCTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAg                                                                                                   >  3:37332/1‑149 (MQ=255)                                                    aaCTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAg                                                                                                   >  3:53613/1‑149 (MQ=255)                                                         ctggagcccctgcgctttggcgtggccggcggcagcctggcctccaacatccgcctggacgggcgcagcgtgcccttgcaggggcTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTgcg                                                                                              >  3:243663‑M2/86‑149 (MQ=255)                                                         ctggagcccctgcgctttggcgtggccggcggcagcctggcctccaacatccgcctggacgggcgcagcgtgcccttgcaggggcTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTgcg                                                                                              >  4:303436‑M2/86‑149 (MQ=255)                                                         cgcgcgCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGACAGTTgcg                                                                                              >  7:102276/1‑149 (MQ=255)                                                                ccctgcgctttggcgtggccggcggcagcctggcctccaacatccgcctggacgggcgcagcgtgcccttgcaggggcTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCtgtg                                                                                       <  4:243663‑M2/71‑1 (MQ=255)                                                                CCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCtgtg                                                                                       <  5:104177‑M1/149‑1 (MQ=255)                                                                  CTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCa                                                                                     >  5:309155‑M1/1‑149 (MQ=255)                                                                      aaGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGAt                                                                                 <  4:183325/149‑1 (MQ=255)                                                                            aCGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCtt                                                                           >  6:133397/1‑149 (MQ=255)                                                                              gtggccggcggcagcctggcctccaacatccgcctggacgggcgcagcgtgcccttgcaggggcTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGc                                                                         <  3:303436‑M2/85‑1 (MQ=255)                                                                              gCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGc                                                                         <  4:53613/149‑1 (MQ=255)                                                                              gCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGc                                                                         >  1:41237/1‑149 (MQ=255)                                                                                  ctcccGTCCAGTCTGATGTTCGAGGCCAGGCTGCCGCCGGCCACGCAAAAGCGCAGGGGCTCCAGGCGTAACAGGCGGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAAGGGCAGTTGCGCGCTCTGCACGATGCGCTTGCCa                                                                       <  6:284477/145‑1 (MQ=255)                                                                                  cgcCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCa                                                                       <  2:41237/147‑1 (MQ=255)                                                                                   gcCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATg                                                                                <  1:131074/137‑1 (MQ=255)                                                                                   gcCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATg                                                                                >  2:131074/1‑137 (MQ=255)                                                                                   gcCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGc                                                                     >  5:311362/1‑148 (MQ=255)                                                                                    cccGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCa                                                                    >  5:64881/1‑148 (MQ=255)                                                                                     ccGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCa                                                                    >  8:207523/1‑147 (MQ=255)                                                                                           aGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCATAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCaaa                                                                  <  8:155225/143‑1 (MQ=255)                                                                                           aGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCATAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCaaa                                                                  <  8:150904/143‑1 (MQ=255)                                                                                           aGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCATAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCaaa                                                                  >  7:155225/1‑143 (MQ=255)                                                                                           aGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCATAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCaaa                                                                  >  7:150904/1‑143 (MQ=255)                                                                                             gCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACAt                                                          <  8:9480/149‑1 (MQ=255)                                                                                               ggATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACAt                                                          >  4:47318/1‑147 (MQ=255)                                                                                                   gTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGt                                                    >  7:46262/1‑149 (MQ=255)                                                                                                       gAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAgg                                                                >  4:80645/1‑133 (MQ=255)                                                                                                       gAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAgg                                                                <  3:80645/133‑1 (MQ=255)                                                                                                       gAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCAt                                                >  4:151576/1‑149 (MQ=255)                                                                                                          gCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGc                                             >  6:344701/1‑149 (MQ=255)                                                                                                          gCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGc                                             <  2:119412/149‑1 (MQ=255)                                                                                                           ccAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGcc                                            <  4:37332/149‑1 (MQ=255)                                                                                                            cAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGccc                                           <  8:46262/149‑1 (MQ=255)                                                                                                                   ccgccgGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCaa                                      <  6:104177/147‑1 (MQ=255)                                                                                                                      ccgGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTc                                 <  3:151576/149‑1 (MQ=255)                                                                                                                         gCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGt                              <  3:262507/149‑1 (MQ=255)                                                                                                                            aCGCCAAAGCGCAGGGGCTCCAGGCGTAACAGCCCGTCTTCGAGCACCACATGGCCCGACACGTCAGAGAACCGCAGATGCGCGCAGTGCACGATGCGCTTGCCAGCAAAGGTGACAAAAGCGTCCATTGCCCGCCAACGCTCGGTGc                            >  2:251415/1‑148 (MQ=255)                                                                                                                             cGCCAAAGCGCAGGGGCTCCAGGCATAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCg                           >  6:6734/1‑148 (MQ=255)                                                                                                                                 aaaGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACg                      >  2:338080/1‑149 (MQ=255)                                                                                                                                 aaaGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACg                      >  1:238949/1‑149 (MQ=255)                                                                                                                                 aaaGCACAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACg                      >  7:25055/1‑149 (MQ=255)                                                                                                                                  aaGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCTACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGc                     <  8:316510/149‑1 (MQ=255)                                                                                                                                                cAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGgcagca       >  4:203370/1‑149 (MQ=255)                                                                                                                                                cAGGCGTAACAGGCCGTCTTCGAGCACAACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATAGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGgcagca       >  7:142212/1‑149 (MQ=255)                                                                                                                                                 aGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGGCAGCAc      >  8:66495/1‑149 (MQ=255)                                                                                                                                                   gCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGGCAGCACCt    >  2:208084/1‑149 (MQ=255)                                                                                                                                                     gTAACAGGCGGTCTTGGAGCACCACATGGGCCGTCAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGGCAGCACCTtg  <  8:142212/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATTGGCGCGAAGGTAGGGAACAGTTGCTTGAGCTTGAAGCCACGCGCAGTCAACTTCGCGCGCCCCTGCAAGGGCACGCTGCGCCCGTCCAGGCGGATGTTGGAGGCCAGGCTGCCGCCGGCCACGCCAAAGCGCAGGGGCTCCAGGCGTAACAGGCCGTCTTCGAGCACCACATGGGCCGACAGGTCAGTGAACGGCAGTTGCGCGCTGTGCACGATGCGCTTGCCAGCAAAGGTGACATCAGCGTCCATTGCCCGCCAACGCTCGGTGCGAAACGCCTCCACTGGCAGCACCTTG  >  NC_002947/5523416‑5523712

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