Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,245	T→C	T116A (ACG→GCG)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,245	0	T	C	100.0%	201.3 / NA	58	T116A (ACG→GCG)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/29);  total (29/29)

CCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAA  >  NC_000913/3508097‑3508391                                                                                                                                                     |                                                                                                                                                   cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCgcgc                                                                                                                                                  >  1:563287/1‑151 (MQ=255) cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCgcgc                                                                                                                                                  >  1:491434/1‑151 (MQ=255)         gCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATa                                                                                                                                          <  1:179818/151‑1 (MQ=255)           tGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATAcc                                                                                                                                        >  2:428208/1‑151 (MQ=255)                gacgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGtt                                                                                                                                   <  1:630646/151‑1 (MQ=255)                gacgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGtt                                                                                                                                   <  2:369965/151‑1 (MQ=255)                  cgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCa                                                                                                                                 <  2:376600/151‑1 (MQ=255)                    tATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAt                                                                                                                               >  1:848736/1‑151 (MQ=255)                     aTGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAtt                                                                                                                              >  2:34857/1‑151 (MQ=255)                     aTGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAtt                                                                                                                              >  2:489526/1‑151 (MQ=255)                     aTGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAt                                                                                                                               <  1:372231/150‑1 (MQ=255)                      tGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAAtt                                                                                                                              <  1:527388/150‑1 (MQ=255)                        cGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCa                                                                                                                           <  1:882336/151‑1 (MQ=255)                               aTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCa                                                                                                                    <  2:895366/151‑1 (MQ=255)                                   ggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCAtctc                                                                                                                >  2:178715/1‑151 (MQ=255)                                          ggTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGcc                                                                                                         <  1:462951/151‑1 (MQ=255)                                           gTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGcc                                                                                                         >  1:285795/1‑150 (MQ=255)                                               ggTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTc                                                                                                    >  1:409663/1‑151 (MQ=255)                                                gTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTcc                                                                                                   >  2:634605/1‑151 (MQ=255)                                                    tGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGc                                                                                                >  2:285434/1‑150 (MQ=255)                                                     gcgcCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCAc                                                                                              <  1:221951/151‑1 (MQ=255)                                                       gcCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGc                                                                                            >  1:473696/1‑151 (MQ=255)                                                            cgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCgg                                                                                       <  2:628430/151‑1 (MQ=255)                                                            cgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCgg                                                                                       >  2:714336/1‑151 (MQ=255)                                                               agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCt                                                                                                  >  2:154760/1‑137 (MQ=255)                                                               agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCt                                                                                                  <  1:154783/137‑1 (MQ=255)                                                               agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTc                                                                                    <  1:174485/151‑1 (MQ=255)                                                                gcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTc                                                                                    <  2:569972/150‑1 (MQ=255)                                                                   gCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGcc                                                                                 >  2:628166/1‑150 (MQ=255)                                                                    cAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcac                                                                               >  1:130256/1‑151 (MQ=255)                                                                      ataaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACAt                                                                             <  2:89650/151‑1 (MQ=255)                                                                        aaataCCTCCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGt                                                                           >  2:655361/1‑151 (MQ=255)                                                                                tctcTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCggg                                                                     <  1:303926/149‑1 (MQ=255)                                                                                tctcTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGa                                                                    <  2:473649/150‑1 (MQ=255)                                                                                 ctctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaa                                                                   >  1:168234/1‑150 (MQ=255)                                                                                   ctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaa                                                                 >  1:320883/1‑150 (MQ=255)                                                                                    ttCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAc                                                               >  1:443925/1‑151 (MQ=255)                                                                                      ccAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaa                                                                 <  1:789750/147‑1 (MQ=255)                                                                                      ccAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgc                                                             >  2:210707/1‑151 (MQ=255)                                                                                      ccAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgc                                                             >  1:592616/1‑151 (MQ=255)                                                                                      ccAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcg                                                              <  2:568744/150‑1 (MQ=255)                                                                                             cgTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTgg                                                      >  1:299009/1‑151 (MQ=255)                                                                                               tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                                                     >  1:442332/1‑150 (MQ=255)                                                                                                     ttCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaataa                                                <  2:781237/149‑1 (MQ=255)                                                                                                     ttCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAAc                                               <  1:477675/150‑1 (MQ=255)                                                                                                           cGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTAc                                         <  1:316410/150‑1 (MQ=255)                                                                                                            gCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAg                                       <  1:726488/151‑1 (MQ=255)                                                                                                               ggTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAgg                                      <  2:469592/149‑1 (MQ=255)                                                                                                                  gCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGccacc                                 >  1:894058/1‑151 (MQ=255)                                                                                                                      aTCTCCGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGt                             >  1:390900/1‑151 (MQ=255)                                                                                                                          ccGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTg                           <  2:123748/149‑1 (MQ=255)                                                                                                                          ccGCGATGATCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGAt                         >  1:381086/1‑151 (MQ=255)                                                                                                                               atgatCCCGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCg                    >  2:20920/1‑151 (MQ=255)                                                                                                                                    cccGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGtctc                <  2:292293/150‑1 (MQ=255)                                                                                                                                     ccGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTcgc              >  2:844846/1‑151 (MQ=255)                                                                                                                                      cGGCTTTCAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTcgcg             >  2:496207/1‑151 (MQ=255)                                                                                                                                             cAGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTAc       <  2:593174/150‑1 (MQ=255)                                                                                                                                              aGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACAt     >  1:34417/1‑151 (MQ=255)                                                                                                                                              aGCTCCGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACAt     <  2:787648/151‑1 (MQ=255)                                                                                                                                                  ccGCGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCaa  <  1:754840/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAA  >  NC_000913/3508097‑3508391

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑