Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,697	G→T	5.1%	A164A (GCC→GCA)	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,697	0	G	T	5.1%	314.1 / 5.0	97	A164A (GCC→GCA)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (39/53);  new base T (2/3);  total (41/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

TGGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAA  >  NC_002947/1235550‑1235843                                                                                                                                                     |                                                                                                                                                   tGGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACg                                                                                                                                                    <  3:374223/148‑1 (MQ=255) tGGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACg                                                                                                                                                    <  1:258962/148‑1 (MQ=255)   gCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGcc                                                                                                                                                 <  2:60162/149‑1 (MQ=255)   gCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGcc                                                                                                                                                 <  5:129156/149‑1 (MQ=255)     ccGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCa                                                                                                                                                >  6:436050/1‑148 (MQ=255)      cGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACgg                                                                                                                                                   >  5:162682/1‑144 (MQ=255)       gAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCAt                                                                                                                                               >  3:292863/1‑147 (MQ=255)       gAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCAt                                                                                                                                               >  6:10434/1‑147 (MQ=255)        aaaTCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAg                                                                                                                                            <  2:306747/149‑1 (MQ=255)         aaTCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGa                                                                                                                                           <  2:319631/149‑1 (MQ=255)          aTCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGaa                                                                                                                                          <  3:305984/149‑1 (MQ=255)               gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATgg                                                                                                                                       >  2:25329/1‑147 (MQ=255)               gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGGCGGTCAGGGCCATGATCACGGCCATCAGAATGGCa                                                                                                                                     <  6:191137/149‑1 (MQ=255)                     ggTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGAt                                                                                                                                >  3:304083/1‑148 (MQ=255)                     ggTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATc                                                                                                                               <  7:218179/149‑1 (MQ=255)                       ttCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGC‑GGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                            <  1:252118/149‑3 (MQ=255)                          aTCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatgg                                                                                                                           <  5:431835/148‑4 (MQ=255)                          aTCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatgg                                                                                                                           <  1:67499/148‑4 (MQ=255)                           tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                         >  7:204107/1‑149 (MQ=255)                           tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                         >  4:101420/1‑149 (MQ=255)                             gCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACtt                                                                                                                        >  3:247498/1‑148 (MQ=255)                             gCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACtt                                                                                                                        >  3:106033/1‑148 (MQ=255)                             gCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTc                                                                                                                       <  5:436050/149‑1 (MQ=255)                               aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCaa                                                                                                                     <  4:143212/149‑1 (MQ=255)                               aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCaa                                                                                                                     <  4:149376/149‑1 (MQ=255)                                   ccAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                    >  3:321147/1‑146 (MQ=255)                                     aGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGa                                                                                                                >  2:124287/1‑148 (MQ=255)                                     aGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGa                                                                                                                >  6:242161/1‑148 (MQ=255)                                       tCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                             <  3:337653/149‑1 (MQ=255)                                       tCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                             >  2:391966/1‑149 (MQ=255)                                        cAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGTCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACg                                                                                                            >  8:85632/1‑149 (MQ=255)                                           tGGGCCCATGGCGCGGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCct                                                                                                         <  6:37438/149‑1 (MQ=255)                                               cccATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggct                                                                                                     <  4:313091/149‑1 (MQ=255)                                                 cATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggctgg                                                                                                   >  4:61749/1‑149 (MQ=255)                                                  aTGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGGCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCtggt                                                                                                  <  5:274081/149‑1 (MQ=255)                                                      cgcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                             >  5:210213/1‑134 (MQ=255)                                                      cgcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                             <  6:210213/134‑1 (MQ=255)                                                       gcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGa                                                                                               >  7:391862/1‑147 (MQ=255)                                                        cTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAg                                                                                            <  4:321147/149‑1 (MQ=255)                                                         tGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGaaa                                                                                                              <  7:308297/130‑1 (MQ=255)                                                         tGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGaaa                                                                                                              >  8:308297/1‑130 (MQ=255)                                                           ccAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCg                                                                                          <  3:168770/148‑1 (MQ=255)                                                            cAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGTACAGCGAc                                                                                        >  3:197607/1‑149 (MQ=255)                                                              aTCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGTCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGTCgg                                                                                      <  1:391966/149‑1 (MQ=255)                                                               tCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACggg                                                                                     <  8:204107/149‑1 (MQ=255)                                                                 aCTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTg                                                                                   >  5:238125/1‑149 (MQ=255)                                                                   tGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCg                                                                                 <  5:63662/149‑1 (MQ=255)                                                                    gCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCgg                                                                                >  1:52696/1‑149 (MQ=255)                                                                    gCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCgg                                                                                >  5:378331/1‑149 (MQ=255)                                                                      aaTCAGGAGGCCGATCAGCAGCGGGGCCATGGCGCCCCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGt                                                                              <  5:417492/149‑1 (MQ=255)                                                                      aaTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGt                                                                              <  4:304083/149‑1 (MQ=255)                                                                      aaTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGt                                                                              <  4:292863/149‑1 (MQ=255)                                                                           gcagGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTg                                                                         <  1:25329/149‑1 (MQ=255)                                                                              ggCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGgag                                                                      <  6:210573/149‑1 (MQ=255)                                                                              ggCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCACGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGgag                                                                      <  5:203654/149‑1 (MQ=255)                                                                               gCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGt                                                                                                                                   >  8:203787/1‑87 (MQ=255)                                                                               gCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGt                                                                                                                                   <  7:203787/87‑1 (MQ=255)                                                                                 cGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAac                                                                   >  3:251439/1‑149 (MQ=255)                                                                                  gATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGTGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGtcggaaga                                                                   >  1:59825/1‑142 (MQ=255)                                                                                  gATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCggg                                                                               <  2:59825/136‑1 (MQ=255)                                                                                  gATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAaca                                                                  >  2:297459/1‑149 (MQ=255)                                                                                   aTCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAacac                                                                 >  5:94913/1‑149 (MQ=255)                                                                                    tCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACAcc                                                                <  1:193682/149‑1 (MQ=255)                                                                                       gcagcGGGGCCATGGTGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGTACAGNGACGGGTGCGGGTAGGTGGAGAACACCGAg                                                             <  1:35125/149‑1 (MQ=255)                                                                                       gcagcGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAg                                                             <  1:297459/149‑1 (MQ=255)                                                                                          gcGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGc                                                           >  8:192185/1‑148 (MQ=255)                                                                                               gCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTc                                                     >  5:2096/1‑149 (MQ=255)                                                                                                     gcgcCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCt                                               <  6:378331/149‑1 (MQ=255)                                                                                                      cgccgcGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTg                                              >  7:339015/1‑149 (MQ=255)                                                                                                       gccgcGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGc                                             <  6:94913/149‑1 (MQ=255)                                                                                                         cgcgcgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCt                                           >  6:409160/1‑149 (MQ=255)                                                                                                          gcgcgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAAGACCGAGGCCAGCTCCAGGCTGCCtt                                          >  7:171991/1‑149 (MQ=255)                                                                                                          gcgcgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCtt                                          <  4:247498/149‑1 (MQ=255)                                                                                                          gcgcgGCAG‑GCCGTTGTTGGCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGGGGAGAACACCGAGGCCAGCTCCAGGCTGCCtt                                          <  6:201538/149‑1 (MQ=255)                                                                                                           cgcgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTc                                         <  6:423931/149‑1 (MQ=255)                                                                                                             cgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGc                                       <  6:238125/149‑1 (MQ=255)                                                                                                             cgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGc                                       <  4:348317/149‑1 (MQ=255)                                                                                                                 aG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAGCAGAAACGCCTGGCTGGTGGATTGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTAcc                                   <  1:124287/149‑1 (MQ=255)                                                                                                                 aG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTAcc                                   >  4:384934/1‑149 (MQ=255)                                                                                                                       gttgttGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                              <  8:339015/149‑1 (MQ=255)                                                                                                                       gttgttGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                              <  3:101420/149‑1 (MQ=255)                                                                                                                              tcgtcgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGc                       <  6:435768/149‑1 (MQ=255)                                                                                                                                gtcgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGt                     >  5:96069/1‑149 (MQ=255)                                                                                                                                 tcgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAAGACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTg                    <  8:171991/149‑1 (MQ=255)                                                                                                                                        gggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGt             <  8:391862/149‑1 (MQ=255)                                                                                                                                        gggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGt             <  7:192185/149‑1 (MQ=255)                                                                                                                                         ggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGTGTAGGTGGAGAACACCGAGGCCAGCTACAGGCTGCCTTCTCTACCGCGCAGCATTGCGTGGGCTTGtt            <  4:197607/149‑1 (MQ=255)                                                                                                                                         ggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGtt            >  2:87161/1‑149 (MQ=255)                                                                                                                                           ccATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTc           >  6:247328/1‑148 (MQ=255)                                                                                                                                            ggtGATCACTGCCATTCTGATGGCCGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCGACCGCGCAGCATGGCGTGGGCTTGTTCg          >  5:387522/3‑148 (MQ=255)                                                                                                                                             gtGATCACTGCCATTCTGATGGCCGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGa         >  5:307388/2‑148 (MQ=255)                                                                                                                                             aTGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGa         >  2:337117/1‑148 (MQ=255)                                                                                                                                                aTCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTACGAAAt      <  5:10434/148‑1 (MQ=255)                                                                                                                                                 tCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATc     >  3:198610/1‑148 (MQ=255)                                                                                                                                                  cACGGCCATCAGAATGTCAGTGATCACCACTTCAACCAGAAAGGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCg    <  3:345256/148‑1 (MQ=255)                                                                                                                                                  cACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgcagc                            >  4:46824/1‑124 (MQ=255)                                                                                                                                                  cACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgcagc                            <  3:46824/124‑1 (MQ=255)                                                                                                                                                   aCGGTCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaa  <  7:85632/149‑1 (MQ=255)                                                                                                                                                   aCGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaa  <  8:17725/149‑1 (MQ=255)                                                                                                                                                   aCGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCga   <  5:242161/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAA  >  NC_002947/1235550‑1235843

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