Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,229,624	C→A	10.1%	T246K (ACG→AAG)	ygjK →	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,229,624	0	C	A	10.1%	167.4 / 4.7	62	T246K (ACG→AAG)	ygjK	alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (38/17);  new base A (2/5);  total (40/22)
Fisher's exact test for biased strand distribution p-value = 8.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGAAAGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAA  >  NC_000913/3229475‑3229767                                                                                                                                                      |                                                                                                                                                ccTGAAAGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCAc                                                                                                                                                 <  1:1222360/150‑1 (MQ=255)      aaGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGctct                                                                                                                                            <  2:1125439/150‑1 (MQ=255)      aaGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGctct                                                                                                                                            <  2:422737/150‑1 (MQ=255)      aaGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTATACCAGTAAGGCGAATATCAACGGTTCGACCACGctct                                                                                                                                            >  1:59842/1‑150 (MQ=255)          cACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTCCCGGTGCAGCCTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTacac                                                                                                                                        >  2:349581/1‑150 (MQ=255)                 ggCAAAGTGCGCGCCACGTGGGATCTGCTGACCTCCGGGGAATCCGAATATCAGGTGCATAAATCCCTGCCGGTGCAGCCTGAAAGCAACGGTAGTCGCTTTTGGAGTAAGGCGGATATCACCGGTTCGCCCCAGCGCTccaccagcgat                                                                                                                                 >  1:295875/1‑145 (MQ=255)                 ggCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGCCTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCAAGCTCTACACCACCTAt                                                                                                                                 >  1:535267/1‑150 (MQ=255)                    aaaGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTcc                                                                                                                              <  1:1701561/150‑1 (MQ=255)                        tGCGCGCCACCTGGGATCTGCTGCCCTCCGGCGAATAGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACTGTTCGACCACGCTCTACACCACCTATTCACAta                                                                                                                          >  2:536998/1‑149 (MQ=255)                           gcgcCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATctgc                                                                                                                       >  1:1507972/1‑150 (MQ=255)                            cgcCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATctgct                                                                                                                      >  2:1021892/1‑150 (MQ=255)                              ccACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGa                                                                                                                    >  2:1778292/1‑150 (MQ=255)                               cACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATAACAACGGTTCGACCCCGCTCTTCACCACCTATTCCCATCTGCTGAc                                                                                                                   >  1:1271272/1‑150 (MQ=255)                                  cTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCACCGGTAATCGCTTGACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGAcctc                                                                                                                >  2:1013541/1‑148 (MQ=255)                                          tgctgACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAg                                                                                                        <  1:1407731/150‑1 (MQ=255)                                           gctgACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGt                                                                                                       <  1:1117357/150‑1 (MQ=255)                                            cttaCCTCCCGCTAATCGGCATATGAGGTGCATAAATACCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTTATTTTTTTTTTAATTTTAAGGGGTCGACCTAGATCTACACAACCTATTCCCATCTGCTGACCGCTCAGGAAGtt                                                                                                      <  2:808613/147‑1 (MQ=255)                                            ctgACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGtt                                                                                                      >  2:922187/1‑150 (MQ=255)                                                    cGGCGAATCGGAATCTCAGGTGCATAAATCCCTGCCGGTTCAGACTGAAATCAACGGTAATCGCTTTAATTTTTTTTTTAATTTTAAGGGTTCGACCTAGATCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAaga                                                                                              <  2:795046/150‑1 (MQ=255)                                                     ggCGACTCGGAATATCAGGTGCATAAATCCCTGCCGTTGCAGACTGAAATCAACGGTAATCGCTTTAAATGTTATTATAATTTTAAGGGTTGGACCAAGATCTACACAACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAagag                                                                                             <  2:1912507/150‑1 (MQ=255)                                                         aaTCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCaaa                                                                                         >  1:1772698/1‑150 (MQ=255)                                                          aTCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAAt                                                                                        >  1:1832145/1‑150 (MQ=255)                                                              gAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTAGCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAg                                                                                    <  1:1346348/150‑1 (MQ=255)                                                                      ggTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCATTAATTATAATTTTAACGGGTCGACCTAGATCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGa                                                                            <  2:493205/150‑1 (MQ=255)                                                                         gCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGatat                                                                         >  2:671760/1‑150 (MQ=255)                                                                         gCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTAAACCACCTATTCCCATCTGATGACCGCTCAGGAAGTTATCAAAGAGCAAATGCAGATCCGCGatat                                                                         >  2:1277957/1‑150 (MQ=255)                                                                          cATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATAtt                                                                        >  1:1535519/1‑150 (MQ=255)                                                                             aaaTCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTg                                                                     >  2:68002/1‑150 (MQ=255)                                                                              aaTCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTgg                                                                    >  1:387561/1‑150 (MQ=255)                                                                               aTCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCATTAAGGCGCATATCAACGGTTCGACCACTCTCTACACCACCTATTCCCATCTGCTGACCGCTCATGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGgc                                                                   <  1:536998/150‑1 (MQ=255)                                                                                 cccTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCATTAAGGCGCATATCAACGGTTAGACAACTCTCTACACCACCTATTCCCATCTGATGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGgcgc                                                                 >  1:1159573/1‑150 (MQ=255)                                                                                   cTGCCGGTGCAGACTGAAATCAACGGGAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGt                                                               >  1:1195204/1‑150 (MQ=255)                                                                                        ggTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGc                                                          >  2:390296/1‑150 (MQ=255)                                                                                           gCAGACTGAAATCACCGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGCCCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGtt                                                       >  1:1138670/1‑150 (MQ=255)                                                                                           gCAGACTGAAATCAACGGTAATCGCTTTACGATTATTTATCATTTAAAGGGGTCGACCTAGATCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGtt                                                       <  2:1832658/150‑1 (MQ=255)                                                                                                cTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCACAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTAta                                                  >  1:1548967/1‑149 (MQ=255)                                                                                                 tGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTAtct                                                 >  1:1588549/1‑150 (MQ=255)                                                                                                 tGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGGTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTAtct                                                 <  1:816143/150‑1 (MQ=255)                                                                                                      tCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCg                                            >  1:1320077/1‑150 (MQ=255)                                                                                                       cAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGc                                           >  2:363665/1‑150 (MQ=255)                                                                                                        aaCGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGcc                                          >  2:473334/1‑150 (MQ=255)                                                                                                          cGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTc                                        >  1:385214/1‑150 (MQ=255)                                                                                                            gTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCgc                                      <  1:1875857/150‑1 (MQ=255)                                                                                                             tAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCgca                                     >  2:1584617/1‑150 (MQ=255)                                                                                                              aaTCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCgcag                                    >  1:1207227/1‑150 (MQ=255)                                                                                                              aaTCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCgcag                                    <  1:427554/150‑1 (MQ=255)                                                                                                              aaTCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCgcag                                    <  1:420446/150‑1 (MQ=255)                                                                                                                   cTTTCCCAGTCCGGCGCATATCAACGGGGCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAACTGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACg                               <  1:418408/150‑1 (MQ=255)                                                                                                                   cTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACg                               <  1:499637/150‑1 (MQ=255)                                                                                                                     ttACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCt                             <  1:398980/150‑1 (MQ=255)                                                                                                                         cAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGga                         <  1:1888525/150‑1 (MQ=255)                                                                                                                         cAGTAAGGAGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGga                         >  1:89402/1‑150 (MQ=255)                                                                                                                              aGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAAt                    <  1:390296/150‑1 (MQ=255)                                                                                                                                 cgcATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATAcc                 >  1:540957/1‑150 (MQ=255)                                                                                                                                   cATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTg               <  2:387561/150‑1 (MQ=255)                                                                                                                                       tCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTgaaga           >  2:612984/1‑150 (MQ=255)                                                                                                                                        cAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTgaagaa          >  2:1376481/1‑150 (MQ=255)                                                                                                                                        cAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTgaagaa          >  2:616685/1‑150 (MQ=255)                                                                                                                                         aaCGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTTGGAAGAATACCTGAAGaaa         >  2:1822304/1‑150 (MQ=255)                                                                                                                                           cGGTTCGACCCCGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTGTAACGCACCGCCTCGCAGCAACGCTGGGAAGACTACCTGAAGAAAgg       >  2:276024/1‑150 (MQ=255)                                                                                                                                                cGACCACGCTCTACACCCCCTATTCCCATCTGCTGACCGCTCAGGTAGGTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGATTTATCTCAAAGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTaa  >  2:1014759/1‑150 (MQ=255)                                                                                                                                                cGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTaa  >  2:609909/1‑150 (MQ=255)                                                                                                                                                cGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACAGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTaa  >  1:399312/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                CCTGAAAGTCACCTTTGGCAAAGTGCGCGCCACCTGGGATCTGCTGACCTCCGGCGAATCGGAATATCAGGTGCATAAATCCCTGCCGGTGCAGACTGAAATCAACGGTAATCGCTTTACCAGTAAGGCGCATATCAACGGTTCGACCACGCTCTACACCACCTATTCCCATCTGCTGACCGCTCAGGAAGTTAGCAAAGAGCAAATGCAGATCCGCGATATTCTGGCGCGTCCGGCGTTTTATCTCACCGCCTCGCAGCAACGCTGGGAAGAATACCTGAAGAAAGGGTTAA  >  NC_000913/3229475‑3229767

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