Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,712	A→C	5.5%	T159T (ACT→ACG)	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,712	0	A	C	5.5%	287.5 / 6.2	90	T159T (ACT→ACG)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base A (37/48);  new base C (2/3);  total (39/51)
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.42e-01

GGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACA  >  NC_002947/1235564‑1235860                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                          >  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                                                                                                                              >  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                                                                                                                              <  6:210213/134‑1 (MQ=255)                                        cgcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                                              >  5:210213/1‑134 (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                                                                                                 >  5:378331/1‑149 (MQ=255)                                                      gCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCgg                                                                                                 >  1:52696/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                                               <  3:101420/149‑1 (MQ=255)                                                                                                         gttgttGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                                               <  8:339015/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)                                                                                                                                                      aGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACa  <  3:61749/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAG‑GCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACA  >  NC_002947/1235564‑1235860

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