Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,648,020	C→T	intergenic (‑232/‑43)	dinQ ← / → agrA	UV inducible membrane toxin DinQ/small RNA AgrA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,648,020	0	C	T	100.0%	146.9 / NA	40	intergenic (‑232/‑43)	dinQ/agrA	UV inducible membrane toxin DinQ/small RNA AgrA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/23);  total (17/23)

TTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCAC  >  NC_000913/3647944‑3648097                                                                             |                                                                              ttATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGt                                                                           >  1:641772/1‑81 (MQ=255)   aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTcc                                                                         >  1:641780/1‑81 (MQ=255)   aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTc                                                                          <  2:641728/80‑1 (MQ=255)     gcgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCtt                                                                       >  1:641778/1‑81 (MQ=255)          ttGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAAt                                                                  <  1:641730/81‑1 (MQ=255)            gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATa                                                                 <  1:641731/80‑1 (MQ=255)             ccAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGt                                                               >  1:641769/1‑81 (MQ=255)              cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGtt                                                              <  2:641731/81‑1 (MQ=255)               agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGttt                                                             <  2:641732/81‑1 (MQ=255)                 agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCa                                                           <  1:641734/81‑1 (MQ=255)                      aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATg                                                      <  1:641735/81‑1 (MQ=255)                          gCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTAt                                                  >  1:641783/1‑81 (MQ=255)                           ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTAtt                                                 <  1:641736/81‑1 (MQ=255)                            cTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTa                                                <  2:641736/81‑1 (MQ=255)                             tCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAg                                               >  1:641775/1‑81 (MQ=255)                              cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAg                                               >  2:641775/1‑80 (MQ=255)                               agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCa                                             >  2:641790/1‑81 (MQ=255)                                 agCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCAtt                                           >  1:641797/1‑81 (MQ=255)                                 agCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCAtt                                           >  1:641796/1‑81 (MQ=255)                                  gCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTg                                          <  1:641738/81‑1 (MQ=255)                                   cTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGa                                         <  2:641738/81‑1 (MQ=255)                                     aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATg                                       >  1:641773/1‑81 (MQ=255)                                      aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGAt                                        <  2:641739/79‑1 (MQ=255)                                      aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGAt                                        <  1:641742/79‑1 (MQ=255)                                      aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATg                                       <  1:641741/80‑1 (MQ=255)                                      aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGc                                      >  1:641795/1‑81 (MQ=255)                                      aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGc                                      <  2:641742/81‑1 (MQ=255)                                               tCATACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAgatga                             <  1:641744/81‑1 (MQ=255)                                                 aTACAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAgatgat                            <  2:641744/80‑1 (MQ=255)                                                   aCAGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGc                         >  2:641809/1‑81 (MQ=255)                                                     aGTACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTa                       <  1:641746/81‑1 (MQ=255)                                                       tACATGCAGGTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTAtc                     <  1:641747/81‑1 (MQ=255)                                                               ggTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTc               >  1:641786/1‑79 (MQ=255)                                                               ggTTATAAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCa             <  2:641747/81‑1 (MQ=255)                                                                   ataAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGt           >  2:641770/1‑79 (MQ=255)                                                                   ataAAACCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc          <  1:641749/80‑1 (MQ=255)                                                                     aaaaCCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa         >  1:641781/1‑79 (MQ=255)                                                                     aaaaCCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa         >  2:641802/1‑79 (MQ=255)                                                                      aaaCCAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCagag      <  2:641749/81‑1 (MQ=255)                                                                          cAGTACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccac  <  1:641751/81‑1 (MQ=255)                                                                             |                                                                              TTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCAC  >  NC_000913/3647944‑3648097

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