Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,709,221	G→C	11.3%	A98A (GCC→GCG)	PP_4166 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,709,221	0	G	C	11.3%	172.6 / 12.9	62	A98A (GCC→GCG)	PP_4166	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base C (4/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

GCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAGGG  >  NC_002947/4709074‑4709368                                                                                                                                                    |                                                                                                                                                    gCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCggg                                                                                                                                                    >  2:297645/1‑149 (MQ=255) gCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                                     <  8:276802/148‑1 (MQ=255)   ccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                                     <  8:168510/146‑1 (MQ=255)   ccATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCgg                                                                                                                                                     >  7:168510/1‑146 (MQ=255)           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                                                                                                         >  6:346091/1‑146 (MQ=255)                                               cATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGc                                                                                                         >  3:351965/1‑146 (MQ=255)                                                aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                     >  5:46722/1‑149 (MQ=255)                                                aTCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGa                                                                                                     >  1:357829/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                                                                                             >  6:376213/1‑140 (MQ=255)                                                                 ggTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGt                                                                                             <  5:376213/140‑1 (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)                                                                                                                                                    |                                                                                                                                                    GCCCATTGTCGGCATCGCGGGGCAGTTCGGCAATGTCCATAAAGGTCATCGGCGAGGGGGATCAGGTCATGCCGAGCGGTCGAAACTGCTCAGCGACCGACGCAGGCAGGCCTGCATGTAGGCCACTTGGCTCTCCAGCGCCTCGCGGGCCAGACGCTGGTCCTGGTAGTCCAGTTGCAGGCGGCGCAAGGCCAGACAACTGGTTTGCAGCAGGTGCGAGACCCGCAGGCAGGCATCCTGCAGCGTACGCAGGTCCTGCTCGCATTGTTGCCTCGCAGCAGGGTCCAGCCCAGGG  >  NC_002947/4709074‑4709368

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