Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,705	G→T	5.3%	L162M (CTG→ATG) ‡	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,705	0	G	T	5.3%	305.7 / 4.5	93	L162M (CTG→ATG) ‡	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (38/50);  new base T (2/3);  total (40/53)
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

AAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTG  >  NC_002947/1235557‑1235853                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                     <  1:67499/148‑4 (MQ=255)                   aTCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatgg                                                                                                                                     <  5:431835/148‑4 (MQ=255)                    tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                                   >  4:101420/1‑149 (MQ=255)                    tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                                   >  7:204107/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:149376/149‑1 (MQ=255)                        aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCaa                                                                                                                               <  4:143212/149‑1 (MQ=255)                            ccAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                              >  3:321147/1‑146 (MQ=255)                              aGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGa                                                                                                                          >  6:242161/1‑148 (MQ=255)                              aGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGa                                                                                                                          >  2:124287/1‑148 (MQ=255)                                tCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                                       >  2:391966/1‑149 (MQ=255)                                tCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                                       <  3:337653/149‑1 (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                                                                                                                        >  8:308297/1‑130 (MQ=255)                                                  tGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGaaa                                                                                                                        <  7:308297/130‑1 (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                                                 <  4:348317/149‑1 (MQ=255)                                                                                                      cgGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGc                                                 <  6:238125/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)                                                                                                                                                   tCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaagaa         >  5:119535/1‑145 (MQ=255)                                                                                                                                                    cAGAATGGCAGAGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGAt    >  6:368720/1‑149 (MQ=255)                                                                                                                                                      gAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTg  <  6:2096/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTG  >  NC_002947/1235557‑1235853

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