Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,750,005	A→G	5.6%	E17E (GAA→GAG)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,750,005	0	A	G	5.6%	203.9 / 2.2	71	E17E (GAA→GAG)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (38/29);  new base G (2/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

GAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTG  >  CP009273/1749857‑1750146                                                                                                                                                     |                                                                                                                                              gAGCCACCACCACTTTCGTAATCCCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGaa                                                                                                                                               <  1:732748/149‑1 (MQ=255)    gCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCAGAagag                                                                                                                                            >  2:195755/1‑149 (MQ=255)       ccaccacTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATg                                                                                                                                         >  1:208422/1‑149 (MQ=255)        caccacTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGAACCCTCGGACCGAAAACCGAATCTGAAGAGgtgt                                                                                                                                        >  2:766711/1‑149 (MQ=255)        caccacTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGAAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGt                                                                                                                                        >  2:807990/1‑149 (MQ=255)         accacTTCCGTAATCCCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGtt                                                                                                                                       <  1:1108593/149‑1 (MQ=255)              tttCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCGTCCTTCTCAACTTACAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTGAGCTa                                                                                                                                  >  1:67977/1‑149 (MQ=255)              tttCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTa                                                                                                                                  >  2:41300/1‑149 (MQ=255)                ttgTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTaaa                                                                                                                                <  2:56185/147‑1 (MQ=255)                 cGTAATACCGGATTCGCTTTCCGGCAGTGCTCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTaaaa                                                                                                                               <  1:1019958/149‑1 (MQ=255)                   tAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCATCTTAAAGACTAAGACTGTCATGAAGAAGACCAAAATTGTTTGCACCATCGTACCGAAATCCGAATCTGAAGAGATGTTAGCTAAAATg                                                                                                                             <  1:804526/149‑1 (MQ=255)                   tAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATg                                                                                                                             >  2:414961/1‑149 (MQ=255)                      tcccGGATTCGCTTTCCTGCAGTGCGCCCATAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTg                                                                                                                          <  1:264893/147‑1 (MQ=255)                      tACCGGATTCGCTTTCCGGCAGTTCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTg                                                                                                                          <  1:1194723/149‑1 (MQ=255)                      tACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTg                                                                                                                          <  1:614054/149‑1 (MQ=255)                      aaCCGGACTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTg                                                                                                                          <  1:435686/148‑1 (MQ=255)                        ccGGATTCGTTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATAGGACAGAAATCCGAATCTGAAGAGTTGTTAGCTAAAATGCTGGa                                                                                                                        <  2:729970/149‑1 (MQ=255)                        ccGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGACTCTGAAGAGATGTTAGCTAAAATGCTGGa                                                                                                                        >  1:1058395/1‑149 (MQ=255)                        ccGGATTCGCGTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTGAGCTAAAATGCTGGa                                                                                                                        <  2:1158271/149‑1 (MQ=255)                             tccGCTTTCCGGCAGTCCGCCCAGAAAGCAAGTTTCTCCCAAGCTCCTCAACTTAAAGACTAAGACTGTCATGAATAAGACCAAAATTGTTTGCACCATCGGACCGAACACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                                                                                                                   <  1:317253/147‑1 (MQ=255)                              tCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAACAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATACTGGACGCTgg                                                                                                                  >  1:875530/1‑149 (MQ=255)                              tCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGCCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGAGGCTgg                                                                                                                  >  2:220554/1‑149 (MQ=255)                               cGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGc                                                                                                                 >  1:960471/1‑149 (MQ=255)                                 cTTTCCGGCAGTGCGCCCAGAAAGCAAGGTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAGATGCTGGACGCTGGCAt                                                                                                               >  1:685884/1‑149 (MQ=255)                                      cGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGCGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg                                                                                                          >  1:87396/1‑149 (MQ=255)                                      cGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg                                                                                                          <  1:567305/149‑1 (MQ=255)                                         ctgTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCAGGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTa                                                                                                       <  1:1063794/147‑1 (MQ=255)                                                 ccAGAAAGCAAGTTTCTCCCATCCTTCCCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTg                                                                                               <  1:455754/149‑1 (MQ=255)                                                     aaaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAGGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACt                                                                                           >  1:16649/1‑149 (MQ=255)                                                     aaaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAGGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCGGAACt                                                                                           >  1:22858/1‑149 (MQ=255)                                                            gTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCCAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtctctcc                                                                                    >  2:245249/1‑148 (MQ=255)                                                               tctcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGGAATGTTAGCTAAAATGCTGGACGCTGACAGGAACAGTATGCGTCTGAACTTTACTCATgg                                                                                 >  2:369535/1‑149 (MQ=255)                                                                 tcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTg                                                                               >  1:303339/1‑149 (MQ=255)                                                                   ccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGAc                                                                             >  1:158726/1‑149 (MQ=255)                                                                     aTCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATATGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTa                                                                           >  1:968724/1‑149 (MQ=255)                                                                         ttCTCAACTTAAAGACTAAGACTGTCATGAGAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCa                                                                       <  1:460985/149‑1 (MQ=255)                                                                           ctcAACTTAAAGATTTAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATTCTTGACGCTGACATGAACGTTATTCGCCTTAACTTCTCTCATGGTTACTATGCCGa                                                                     >  2:550662/1‑149 (MQ=255)                                                                              aaCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAaca                                                                  >  1:560088/1‑149 (MQ=255)                                                                              aaCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCAACGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAaca                                                                  >  2:1188216/1‑149 (MQ=255)                                                                                 ttAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACgg                                                               >  2:711093/1‑149 (MQ=255)                                                                                     agactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAg                                                           <  1:41300/149‑1 (MQ=255)                                                                                      gactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAgc                                                          <  2:285256/149‑1 (MQ=255)                                                                                      gactaagactGACATGAAAAAGACCAAAATTGTATGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCGgc                                                          <  2:888148/149‑1 (MQ=255)                                                                                         taagactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGACTGGACTTCTCTCATGGTGACTATGCAGAACACGGTCAgcgaa                                                       >  2:639326/1‑147 (MQ=255)                                                                                         gaagacaGGCATGAAAAAGACCAAAATCGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCa                                                       <  1:1026289/148‑1 (MQ=255)                                                                                            gactGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTc                                                    >  2:796878/1‑149 (MQ=255)                                                                                             actGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCa                                                   <  1:226439/149‑1 (MQ=255)                                                                                             actGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCa                                                   >  1:568515/1‑149 (MQ=255)                                                                                             actGTCAGGAAAAAGACCAAAACTGTTTGCACCATCGGCCCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCa                                                   <  1:712530/149‑1 (MQ=255)                                                                                                    tGAAAAAGCCCATATTTTTTTGAACCTTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTGCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTg                                            <  1:241292/149‑1 (MQ=255)                                                                                                         aaGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCCGAATCTGCGCaa                                       >  2:72744/1‑149 (MQ=255)                                                                                                         aaGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCaa                                       >  1:496667/1‑149 (MQ=255)                                                                                                          aGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAc                                      <  2:465583/149‑1 (MQ=255)                                                                                                          aGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTCCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAAc                                      <  1:929859/149‑1 (MQ=255)                                                                                                              cAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGCGATGTTAGCTAAAATGCAGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtga                                  <  1:645166/149‑1 (MQ=255)                                                                                                              cAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtga                                  <  1:770312/149‑1 (MQ=255)                                                                                                              cAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtga                                  <  2:556542/149‑1 (MQ=255)                                                                                                                 aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtgatga                               >  2:534463/1‑149 (MQ=255)                                                                                                                 aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtgatga                               >  2:534925/1‑149 (MQ=255)                                                                                                                 aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGggacga                               >  2:934741/1‑149 (MQ=255)                                                                                                                 aaTTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCAGCGCAACGggatga                               >  2:933991/1‑149 (MQ=255)                                                                                                                 aaTTGTTTGCACAATCGGAACGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGtgatgg                               >  2:762096/1‑148 (MQ=255)                                                                                                                     gtgtGCACCATCGGACCGAAACCCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCaa                           <  1:195755/146‑1 (MQ=255)                                                                                                                         gCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAACCt                       >  2:30286/1‑149 (MQ=255)                                                                                                                         gCACCATCGGACCGAAAACCGAATCCGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACt                       >  1:624304/1‑149 (MQ=255)                                                                                                                                 ggACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAc               >  2:1172824/1‑149 (MQ=255)                                                                                                                                  gACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAcc              >  1:58347/1‑149 (MQ=255)                                                                                                                                     cGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgcc           >  2:253620/1‑149 (MQ=255)                                                                                                                                      gAAAACCGAATCTGAGGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAGccgccg          <  2:22858/149‑1 (MQ=255)                                                                                                                                      gAAAACCGAATCTGAGGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAAccgccg          <  2:16649/149‑1 (MQ=255)                                                                                                                                      gAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAATccgccg          <  2:1058395/149‑1 (MQ=255)                                                                                                                                              aaTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTg  >  1:689400/1‑149 (MQ=255)                                                                                                                                              aaTCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTg  >  1:949269/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              GAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGACTATGCAGAACACGGTCAGCGCATTCAGAATCTGCGCAACGTGATGAGCAAAACTGGTAAAACCGCCGCTATCCTG  >  CP009273/1749857‑1750146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑