Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,917,690	T→A	8.1%	I30N (ATC→AAC) ‡	PP_4328 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,917,690	0	T	A	8.1%	221.0 / 6.4	74	I30N (ATC→AAC) ‡	PP_4328	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base A (3/3);  total (39/35)
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.91e-01

ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCAGCCAG  >  NC_002947/4917544‑4917838                                                                                                                                                   |                                                                                                                                                     attattACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATcc                                                                                                                                                    <  3:227012/149‑1 (MQ=255)     ttaCACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGgc                                                                                                                                                <  2:299412/149‑1 (MQ=255)         acGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                            <  6:142966/149‑1 (MQ=255)         acGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAg                                                                                                                                            >  3:22349/1‑149 (MQ=255)           gcgcCGCTGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGcc                                                                                                                                          >  6:144562/1‑149 (MQ=255)           gcgcCGCCGCGCTGCCAAGAGGCCTGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGcc                                                                                                                                          >  7:43397/1‑149 (MQ=255)                gccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGt                                                                                                                                     <  7:33867/149‑1 (MQ=255)                 ccgcGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGtt                                                                                                                                    <  4:22349/149‑1 (MQ=255)                   gcgcTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGc                                                                                                                                  <  2:151980/149‑1 (MQ=255)                       tGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAGGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAg                                                                                                                               >  6:183750/1‑148 (MQ=255)                        gCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                            >  2:233867/1‑134 (MQ=255)                        gCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGGATCAGGCGCAg                                                                                                                                            <  1:233867/134‑1 (MQ=255)                         ccAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCt                                                                                                                             >  4:191336/1‑148 (MQ=255)                              agGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGc                                                                                                                       <  1:29930/149‑1 (MQ=255)                                           cAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCTCGCAAg                                                                                                          <  5:326876/149‑1 (MQ=255)                                                    taCAATCATGACTGAAATCAATTCTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCTTGAAGGCGGGCATGTCCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTTGCTCCTGGCGAAACCGCGCAAGCCGAGGt                                                                                                   <  4:192346/147‑1 (MQ=255)                                                    taCAATAATGACTGAAATCAATAGTATCGGCGCACAAACAGCTATCAACAGCCAAGCCAAGAAGGCGGGCATGACCGGCGAACTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTACTGGCGAAACCGCGCAAGACGAGGtta                                                                                                 >  3:192346/1‑148 (MQ=255)                                                        aTCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATc                                                                                             <  1:54588/149‑1 (MQ=255)                                                          cATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGt                                                                                           >  1:10606/1‑149 (MQ=255)                                                            tGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                         <  3:191336/149‑1 (MQ=255)                                                            tGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                         <  5:144562/149‑1 (MQ=255)                                                            tGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTg                                                                                         <  6:185596/149‑1 (MQ=255)                                                              aCTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAGGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGc                                                                                        <  5:183750/148‑1 (MQ=255)                                                               cTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGt                                                                                      >  7:5159/1‑149 (MQ=255)                                                                 gAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCt                                                                                                                             <  6:196964/108‑1 (MQ=255)                                                                 gAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCt                                                                                                                             >  5:196964/1‑108 (MQ=255)                                                                     tCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACg                                                                                >  7:269856/1‑149 (MQ=255)                                                                       aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTgg                                                                                                                        <  4:71242/107‑1 (MQ=255)                                                                       aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                          <  6:192911/137‑1 (MQ=255)                                                                       aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                          >  5:192911/1‑137 (MQ=255)                                                                       aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACggg                                                                              >  3:212005/1‑149 (MQ=255)                                                                       aaTAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACggg                                                                              >  1:50464/1‑149 (MQ=255)                                                                            tctcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                          >  1:188078/1‑132 (MQ=255)                                                                            tctcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGtt                                                                                          <  2:188078/132‑1 (MQ=255)                                                                            tctcGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAg                                                                          >  2:302894/1‑148 (MQ=255)                                                                                ggCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCATGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGa                                                                        <  8:43397/146‑1 (MQ=255)                                                                                  cgcACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCTCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTc                                                                   <  8:5159/149‑1 (MQ=255)                                                                                   gcACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAg                                                                                                     >  2:186927/1‑114 (MQ=255)                                                                                   gcACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAg                                                                                                     <  1:186927/114‑1 (MQ=255)                                                                                       aaaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTAt                                                                                              >  1:35031/1‑117 (MQ=255)                                                                                       aaaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTAt                                                                                              <  2:35031/117‑1 (MQ=255)                                                                                       aaaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTg                                                              >  7:157937/1‑149 (MQ=255)                                                                                        aaCCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTgg                                                             <  2:10606/149‑1 (MQ=255)                                                                                          ccGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTg                                                           >  7:249518/1‑149 (MQ=255)                                                                                           cgcgATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGc                                                          >  8:130785/1‑149 (MQ=255)                                                                                                  aaCAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGc                                                    <  8:249518/148‑1 (MQ=255)                                                                                                    cAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGa                                                 >  3:104785/1‑149 (MQ=255)                                                                                                       ccaagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATg                                               >  8:238545/1‑148 (MQ=255)                                                                                                         aagccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCATGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAg                                            <  3:197060/149‑1 (MQ=255)                                                                                                            ccaTGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGc                                          >  6:32886/1‑148 (MQ=255)                                                                                                               tGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCGGCTGGTGCTACg                                                      >  4:198648/1‑133 (MQ=255)                                                                                                               tGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCGGCTGGTGCTACg                                                      <  3:198648/133‑1 (MQ=255)                                                                                                               tGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCACGccgtt                                      >  7:278525/1‑146 (MQ=255)                                                                                                                   ggCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCg                                  <  2:50464/149‑1 (MQ=255)                                                                                                                            gACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGAcc                          >  1:116317/1‑148 (MQ=255)                                                                                                                             aCCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGAcc                          >  3:269098/1‑147 (MQ=255)                                                                                                                               cGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCGGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGc                      >  7:160075/1‑149 (MQ=255)                                                                                                                               cGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGc                      <  4:104785/149‑1 (MQ=255)                                                                                                                                 gCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCtt                    >  4:106625/1‑149 (MQ=255)                                                                                                                                  cGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTc                   <  8:157937/149‑1 (MQ=255)                                                                                                                                     tCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAgg                >  5:290063/1‑149 (MQ=255)                                                                                                                                         ctgcGCCTGGATCAGGCGCAGCCTGGGTTGGTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGCGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCaa            >  1:29525/1‑149 (MQ=255)                                                                                                                                         ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGt                                                            <  3:45307/101‑1 (MQ=255)                                                                                                                                         ctgcGCCTGGATCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGt                                                            >  4:45307/1‑101 (MQ=255)                                                                                                                                         ctgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAg                 <  8:269856/144‑1 (MQ=255)                                                                                                                                          tgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAAc           >  4:63708/1‑149 (MQ=255)                                                                                                                                          tgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAAc           <  2:116317/149‑1 (MQ=255)                                                                                                                                          tgcGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAAc           >  1:135428/1‑149 (MQ=255)                                                                                                                                           gcgcCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTg                                                  >  7:25951/1‑109 (MQ=255)                                                                                                                                           gcgcCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTg                                                  <  8:25951/109‑1 (MQ=255)                                                                                                                                            cgcCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCACCGGCGTGCAGACCCAGCTTCAGGCCAACgc         >  6:94867/1‑149 (MQ=255)                                                                                                                                              ccTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgcca       >  2:128203/1‑149 (MQ=255)                                                                                                                                                tGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgcca       <  5:287369/147‑1 (MQ=255)                                                                                                                                                tGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACgcca       >  6:287369/1‑147 (MQ=255)                                                                                                                                                   tCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGccagccag  >  2:37878/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ATTATTACACGCGCCGCCGCGCTGCCAAGAGGCCGGGCTCGACAGGTTGTATACAATCATGACTGAAATCAATAGTCTCGGCGCACAAACCGCGATCAACAGCCAAGCCATGAAGGCGGGCATGACCGGCGATCTGCTGCGCCTGATCCAGGCGCAGCCTGGGTTGCTAGCTCCTGGCGAAACCGCGCAAGCCGAGGTTTTATCGTTGCGTCAGACGGGCCAGGACTTTCAGCTGGTGCTACGGCTGATGCAGGCCAACGGCGTGCAGACCCAGCTTCAGGCCAACGCCAGCCAG  >  NC_002947/4917544‑4917838

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