Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,709,227	A→G	10.9%	R96R (CGT→CGC)	PP_4166 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,709,227	0	A	G	10.9%	160.7 / 9.7	65	R96R (CGT→CGC)	PP_4166	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (4/3);  total (35/30)
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.29e-01

GGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAGGGTTGTC  >  NC_002947/4709084‑4709373                                                                                                                                                |                                                                                                                                                   ggCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTg                                                                                                                                               <  1:297645/149‑1 (MQ=255)  gCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTgg                                                                                                                                              <  5:43800/149‑1 (MQ=255)     tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTggaga                                                                                                                                           <  7:132650/149‑4 (MQ=255)     tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGTCTGGCCCGCGAGGCGCTggaga                                                                                                                                           <  2:124355/149‑4 (MQ=255)     tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTcc                                                                                                                                           >  6:237317/1‑149 (MQ=255)     tCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTcc                                                                                                                                           <  5:124241/149‑1 (MQ=255)        cgGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                                        <  1:44684/149‑1 (MQ=255)          gggCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTa                                                                                                                                      <  3:170742/149‑1 (MQ=255)              aGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTc                                                                                                                                   >  4:364375/1‑148 (MQ=255)               gTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGCGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGCCGCGGGTCCTGGTAGTc                                                                                                                                   >  4:308286/1‑147 (MQ=255)                ttCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAg                                                                                                                                >  2:110469/1‑149 (MQ=255)                      aaTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGt                                                                                                                               >  5:390084/1‑144 (MQ=255)                      aaTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGt                                                                                                                               <  6:390084/144‑1 (MQ=255)                       aTGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAgg                                                                                                                         <  8:257851/149‑1 (MQ=255)                             aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGGTGGTCCTGGTAGTCCAGTTGCAGgcggcg                                                                                                                    >  6:380188/1‑148 (MQ=255)                             aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                                        >  3:282451/1‑128 (MQ=255)                             aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTgg                                                                                                                                        <  4:282451/128‑1 (MQ=255)                             aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTc                                                                                                                                   <  1:304576/133‑1 (MQ=255)                             aTAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACAGACGCAGGCAGGCCTGCATGTAGGCAACTTGGCTCTCCAGCGCCTCGAGGGACAGACGCTGGTCCTGGTAGTc                                                                                                                                   >  2:304576/1‑133 (MQ=255)                                     cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                                              >  3:351965/1‑146 (MQ=255)                                     cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                                              >  6:346091/1‑146 (MQ=255)                                      aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          >  1:357829/1‑149 (MQ=255)                                      aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          >  5:46722/1‑149 (MQ=255)                                        cGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGAc                                                                                                         >  3:20394/1‑148 (MQ=255)                                                gggATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGttt                                                                                                <  1:110469/149‑1 (MQ=255)                                                 ggATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTg                                                                                               <  6:46722/149‑1 (MQ=255)                                                       ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGt                                                                                                  <  5:376213/140‑1 (MQ=255)                                                       ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGt                                                                                                  >  6:376213/1‑140 (MQ=255)                                                        gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          >  3:92347/1‑131 (MQ=255)                                                        gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          <  4:92347/131‑1 (MQ=255)                                                        gTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAgg                                                                                        >  8:282697/1‑149 (MQ=255)                                                           aTGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGc                                                                                     <  8:245555/149‑1 (MQ=255)                                                              ccGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCga                                                                                   <  2:243904/148‑1 (MQ=255)                                                                    ggTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          >  7:94700/1‑119 (MQ=255)                                                                    ggTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                          <  8:94700/119‑1 (MQ=255)                                                                      tCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGGATGCAGCAGGTGCGAGACCCgcag                                                                          >  6:77525/1‑149 (MQ=255)                                                                           aCTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGcaggcagg                                                                     >  6:241747/1‑149 (MQ=255)                                                                                tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATcc                                                                >  2:337097/1‑149 (MQ=255)                                                                                tCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATcc                                                                >  1:366302/1‑149 (MQ=255)                                                                                  aGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTg                                                              <  7:282697/149‑1 (MQ=255)                                                                                      ctacatgcaggcctGCCTGCGTGTAGGCCACTTGGCTCTCCAGCGCATCGCCGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGc                                                          >  4:192233/7‑149 (MQ=255)                                                                                      ctacatgcaggcctGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGc                                                          >  1:163585/7‑149 (MQ=255)                                                                                      accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCt                                                               <  3:332111/144‑1 (MQ=255)                                                                                      accgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCt                                                               >  4:332111/1‑144 (MQ=255)                                                                                        cgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCCTCCTGCAGCg                                                         >  6:323042/1‑148 (MQ=255)                                                                                        cgacGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCg                                                         >  3:339911/1‑148 (MQ=255)                                                                                               ggcaggcCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAgg                                                 <  5:346091/149‑1 (MQ=255)                                                                                                     ccTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGc                                           >  3:124619/1‑149 (MQ=255)                                                                                                              aGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgt                                  >  2:146774/1‑149 (MQ=255)                                                                                                               ggCCACTTTGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGTTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                                 <  1:134905/149‑1 (MQ=255)                                                                                                               ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                                 <  5:237317/149‑1 (MQ=255)                                                                                                               ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                                 <  4:84007/149‑1 (MQ=255)                                                                                                               ggCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgtt                                 <  4:351965/149‑1 (MQ=255)                                                                                                                gCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCAttgttg                                >  8:87009/1‑149 (MQ=255)                                                                                                                  cACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGcc                              >  7:9968/1‑149 (MQ=255)                                                                                                                   aCTTGGGTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCt                             <  1:268894/149‑1 (MQ=255)                                                                                                                       ggCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCg                           <  4:20394/147‑1 (MQ=255)                                                                                                                       ccaggaccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCg                           >  8:206140/7‑147 (MQ=255)                                                                                                                         aggaccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGATGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGCCCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGGTCGCATTGTTGCCTCgctc                        >  8:169009/5‑146 (MQ=255)                                                                                                                             ccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAggg                   >  4:201110/1‑149 (MQ=255)                                                                                                                             ccAGCGTCTGGCCCGCGAGGCGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAggg                   >  1:310520/1‑149 (MQ=255)                                                                                                                              cAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAgg                    <  2:163585/147‑1 (MQ=255)                                                                                                                                         gggCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAggg       <  3:192233/149‑1 (MQ=255)                                                                                                                                              aGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTTCTCGCATAGTTGCCTCGCAGCAGGGTCCAGCCCACGGATGTc  >  4:83876/1‑149 (MQ=255)                                                                                                                                              aGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAGGGTTGTc  <  5:380188/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAGGGTTGTC  >  NC_002947/4709084‑4709373

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