Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	331,911	A→G	intergenic (+415/‑460)	betT → / → yahA	choline transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	331,911	0	A	G	100.0%	246.9 / NA	70	intergenic (+415/‑460)	betT/yahA	choline transporter of high affinity/c‑di‑GMP‑specific phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/31);  total (39/31)

TTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCT  >  NC_000913/331778‑332043                                                                                                                                      |                                                                                                                                     ttGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGc                                                                                                                                 <  2:250524/139‑1 (MQ=255)  tGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  1:108127/139‑1 (MQ=255)  tGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  2:116114/139‑1 (MQ=255)         tAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  1:1058207/132‑1 (MQ=255)         tAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                >  2:1058207/1‑132 (MQ=255)         tAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTa                                                                                                                         >  1:187546/1‑139 (MQ=255)                tAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGACATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATAc                                                                                                                  <  2:122512/139‑1 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  1:691267/105‑1 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                >  2:691267/1‑105 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAg                                                                                                                                   >  2:899848/1‑102 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAg                                                                                                                                   <  1:899848/102‑1 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGt                                                                                                                          <  1:593830/111‑1 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGt                                                                                                                          >  2:593830/1‑111 (MQ=255)                                    tttATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCtt                                                                                              >  2:767914/1‑139 (MQ=255)                                        ttGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACt                                                                                          <  2:187546/139‑1 (MQ=255)                                           aaTGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                >  1:1313236/1‑98 (MQ=255)                                           aaTGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  2:1313236/98‑1 (MQ=255)                                            aTGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCAtt                                                                                      >  1:437504/1‑139 (MQ=255)                                                    aaTATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCAtt                                                                              <  2:437504/139‑1 (MQ=255)                                                         tGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAgg                                                                                                                                  >  1:615557/1‑82 (MQ=255)                                                         tGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAgg                                                                                                                                  <  2:615557/82‑1 (MQ=255)                                                                 ttGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCt                                                                 <  1:767914/139‑1 (MQ=255)                                                                    gCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTAtg                                                              <  2:259096/139‑1 (MQ=255)                                                                      aTTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCaa                                                                                                    >  1:1500750/1‑99 (MQ=255)                                                                      aTTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCaa                                                                                                    <  2:1500750/99‑1 (MQ=255)                                                                                   tGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGa                                                                     <  2:1323686/117‑1 (MQ=255)                                                                                   tGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGa                                                                     >  1:1323686/1‑117 (MQ=255)                                                                                   tGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGa                                                                     <  2:988623/117‑1 (MQ=255)                                                                                   tGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGa                                                                     >  1:988623/1‑117 (MQ=255)                                                                                   tGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGgaga                                               >  2:1564612/1‑139 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                <  1:844896/57‑1 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGggcaggca                                                                                                                                >  2:844896/1‑57 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGtt                                                                       >  2:792096/1‑114 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGtt                                                                       <  1:792096/114‑1 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACg                                                                   <  1:1068071/118‑1 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACg                                                                   >  2:1068071/1‑118 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGagag                                              >  1:1123384/1‑139 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGagag                                              >  2:380408/1‑139 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGagag                                              >  1:650414/1‑139 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGagag                                              >  1:1577799/1‑139 (MQ=255)                                                                                    ggTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGagag                                              >  2:1504015/1‑139 (MQ=255)                                                                                     gTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAg                                                                                                                                   <  2:269767/53‑1 (MQ=255)                                                                                     gTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAg                                                                                                                                   >  1:269767/1‑53 (MQ=255)                                                                                     gTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTg                                                                      <  1:1174067/114‑1 (MQ=255)                                                                                     gTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTg                                                                      >  2:1174067/1‑114 (MQ=255)                                                                                     gTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAgg                                             >  2:1411026/1‑139 (MQ=255)                                                                                       tATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCGTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTa                                           >  2:1064665/1‑139 (MQ=255)                                                                                        aTCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTAt                                          >  2:1183624/1‑139 (MQ=255)                                                                                          cATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCt                                                                 >  1:1388414/1‑114 (MQ=255)                                                                                          cATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCt                                                                 <  2:1388414/114‑1 (MQ=255)                                                                                          cATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGgaga                                               <  1:1432571/132‑1 (MQ=255)                                                                                          cATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGgaga                                               >  2:1432571/1‑132 (MQ=255)                                                                                               gCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATaga                                                                                                                            <  2:56753/50‑1 (MQ=255)                                                                                               gCTGAATGGATTCAGTCTTAATGAGTGGGCTTTTAAGGGGCAGGCATaga                                                                                                                            >  1:56753/1‑50 (MQ=255)                                                                                                cTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAAt                                                                                                                       <  2:916696/54‑1 (MQ=255)                                                                                                cTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAAt                                                                                                                       >  1:916696/1‑54 (MQ=255)                                                                                                    aTGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATg                                                                                                             >  1:1323941/1‑60 (MQ=255)                                                                                                    aTGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATg                                                                                                             <  2:1323941/60‑1 (MQ=255)                                                                                                        aTTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGgaga                                               <  1:582430/118‑1 (MQ=255)                                                                                                        aTTCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGgaga                                               >  2:582430/1‑118 (MQ=255)                                                                                                          tCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGCGTAATGATACGTATGCATAACCaa                                                                                                    >  1:1405413/1‑63 (MQ=255)                                                                                                          tCAGTCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCaa                                                                                                    <  2:1405413/63‑1 (MQ=255)                                                                                                              tCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCaa                                                                                                    >  1:1266740/1‑59 (MQ=255)                                                                                                              tCTTAATGAGTGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCaa                                                                                                    <  2:1266740/59‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTTATTATGTCATTTAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGt                      <  1:1343227/127‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTg                                                                             <  1:224947/72‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTg                                                                             >  2:224947/1‑72 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCt                                                                 >  2:747472/1‑84 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCt                                                                 <  1:747472/84‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTc                                      >  2:389943/1‑111 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTc                                      <  1:389943/111‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGa                          >  1:1404671/1‑123 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGa                          <  2:1404671/123‑1 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGt                      >  2:1343227/1‑127 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCAtata          >  1:1588088/1‑139 (MQ=255)                                                                                                                        tGGGTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCAtata          >  1:985652/1‑139 (MQ=255)                                                                                                                          ggTTTTTAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCAtattat        >  2:767061/1‑139 (MQ=255)                                                                                                                                tAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCt  <  2:650414/139‑1 (MQ=255)                                                                                                                                tAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCt  <  2:1123384/139‑1 (MQ=255)                                                                                                                                tAAGGGGCAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGATTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATAGAATGCGCCt  >  1:952340/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                     TTGTTATTTAAGGTTTAATTATCTGTGTGTGATATTTTATTGAATGTTTTAAATATTGTTTTTATTGGCATTGCTATAATATTGGTTATCATTTGCTGAATGGATTCAGTCTTAATGAGTGGGTTTTTAAGGGACAGGCATAGAGTAATGATACGTATGCATAACCAACATCTTTACTCATTATGTCATTGAATGTTGACGCTATGTGTTTATGAGGGAGAGGTATTTTCAGTTGATCTGGATTGTTAAATTCATATAATGCGCCT  >  NC_000913/331778‑332043

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