Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,263,500	G→A	7.4%	P54L (CCG→CTG) ‡	tdcD ←	propionate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,263,500	0	G	A	7.4%	200.1 / 3.8	67	P54L (CCG→CTG) ‡	tdcD	propionate kinase
Reads supporting (aligned to +/- strand):  ref base G (21/41);  new base A (3/2);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAAC  >  NC_000913/3263351‑3263642                                                                                                                                                      |                                                                                                                                               cTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCg                                                                                                                                                >  2:76422/1‑150 (MQ=255)  ttCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCgg                                                                                                                                               <  2:1013779/150‑1 (MQ=255)    gaTCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTg                                                                                                                                             <  1:852111/149‑1 (MQ=255)     aTCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGTGc                                                                                                                                            >  2:148979/1‑150 (MQ=255)        ggTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTgg                                                                                                                                         <  1:260434/150‑1 (MQ=255)        ggTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTgg                                                                                                                                         >  1:33564/1‑150 (MQ=255)         gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                        <  1:1305705/150‑1 (MQ=255)          taataaTGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGgggcgggcg                                                                                                                                       >  1:182912/1‑149 (MQ=255)            ataatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGctct                                                                                                                                     <  1:1154210/150‑1 (MQ=255)            ataatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctca                                                                                                                                     <  2:1301434/150‑2 (MQ=255)            ataatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCACCGCTGGctca                                                                                                                                     <  2:650037/150‑2 (MQ=255)              aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcc                                                                                                                                   <  2:913936/150‑1 (MQ=255)               atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTccc                                                                                                                                  <  1:279383/150‑1 (MQ=255)                      aCTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTa                                                                                                                           <  1:548280/150‑1 (MQ=255)                       cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAc                                                                                                                          <  2:357010/150‑1 (MQ=255)                       cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAc                                                                                                                          <  2:362101/150‑1 (MQ=255)                         cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                        <  1:48319/150‑1 (MQ=255)                         cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                        >  1:1402609/1‑150 (MQ=255)                         cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACgg                                                                                                                        <  1:1088585/150‑1 (MQ=255)                            tAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATa                                                                                                                     >  2:42646/1‑150 (MQ=255)                             aaaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATaa                                                                                                                    >  2:464/1‑150 (MQ=255)                              aaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAg                                                                                                                   <  1:766027/150‑1 (MQ=255)                                aaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGaa                                                                                                                 <  2:280895/150‑1 (MQ=255)                                    cTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCa                                                                                                             >  2:258555/1‑150 (MQ=255)                                      gccgccGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAtt                                                                                                           >  1:1368855/1‑150 (MQ=255)                                       ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATGAGAATGCAttt                                                                                                          >  2:294575/1‑150 (MQ=255)                                           cgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCg                                                                                                      <  2:1113760/150‑1 (MQ=255)                                            gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                     >  2:802000/1‑150 (MQ=255)                                               aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGtt                                                                                                  <  2:199512/150‑1 (MQ=255)                                                gCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTa                                                                                                 <  2:33564/150‑1 (MQ=255)                                                     gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                            <  2:990381/150‑1 (MQ=255)                                                          ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGTCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCgg                                                                                       <  2:109120/150‑1 (MQ=255)                                                          ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCgg                                                                                       <  2:1282348/150‑1 (MQ=255)                                                           gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGc                                                                                      <  1:703630/150‑1 (MQ=255)                                                            ccAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCa                                                                                     >  1:754376/1‑150 (MQ=255)                                                             cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCaa                                                                                    <  1:783375/150‑1 (MQ=255)                                                              aaTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAAt                                                                                   >  2:33630/1‑150 (MQ=255)                                                                 tAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAcc                                                                                <  2:369873/150‑1 (MQ=255)                                                                  aaGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCt                                                                               >  2:616265/1‑150 (MQ=255)                                                                     gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCTGCAATACCTGAc                                                                            <  2:1402609/150‑1 (MQ=255)                                                                      ccACACGGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACa                                                                           <  2:1090032/150‑1 (MQ=255)                                                                       cacaCTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACAt                                                                          >  1:616397/1‑150 (MQ=255)                                                                           cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                      >  1:296755/1‑150 (MQ=255)                                                                           cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                      <  1:500891/150‑1 (MQ=255)                                                                              tCATTTACATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGGCGGCAATACCTGACATTAATACt                                                                   <  1:1056319/150‑1 (MQ=255)                                                                              tCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACt                                                                   <  2:1270742/150‑1 (MQ=255)                                                                                aTTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                 <  2:1368855/150‑1 (MQ=255)                                                                                aTTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTc                                                                 <  2:796524/150‑1 (MQ=255)                                                                                 tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATGAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                <  1:294575/150‑1 (MQ=255)                                                                                  ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                               <  1:76422/150‑1 (MQ=255)                                                                                          cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                       >  1:653951/1‑150 (MQ=255)                                                                                           gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                      <  1:718547/150‑1 (MQ=255)                                                                                               ttCCTGTTCAAATGCATTTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATc                                                  <  1:940296/150‑1 (MQ=255)                                                                                                   aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGc                                              >  1:815015/1‑150 (MQ=255)                                                                                                          aTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGaaa                                       <  2:372028/150‑1 (MQ=255)                                                                                                            gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTAACAGTCGCTGGCATCGAGCACGGAAAAc                                     >  2:1364848/1‑150 (MQ=255)                                                                                                                   ccTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCg                              >  2:74765/1‑150 (MQ=255)                                                                                                                    cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                             <  1:288124/150‑1 (MQ=255)                                                                                                                       cAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTAAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                          <  1:298073/150‑1 (MQ=255)                                                                                                                         aTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACg                        >  2:1138544/1‑150 (MQ=255)                                                                                                                             aCCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAAcc                    <  1:616265/150‑1 (MQ=255)                                                                                                                                ttCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACAGAAAACTTAATCGAAGACGAACcaca                 <  2:653951/150‑1 (MQ=255)                                                                                                                                      gCTGTGGTGAGCCAGCGGTGCTGGCTCTGCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaat           >  1:966177/1‑150 (MQ=255)                                                                                                                                       cTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaata          <  1:148979/150‑1 (MQ=255)                                                                                                                                       cTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaata          <  2:754376/150‑1 (MQ=255)                                                                                                                                         ggggagAGCCAGCACCGCTGGCTCACCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc        >  1:979765/6‑150 (MQ=255)                                                                                                                                           ggagAGCCAGCAACGCTGGCTATCCACAATTTACGGATAAGAATGCATTTTCCGAGTTAATAACGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa      >  2:953145/4‑150 (MQ=255)                                                                                                                                            gagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCAAAGTTAATAACCaa     >  2:651911/3‑150 (MQ=255)                                                                                                                                               aGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaac  <  2:629079/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               CTTCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAAC  >  NC_000913/3263351‑3263642

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