Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,352,826	G→A	P67L (CCG→CTG)	arcB ←	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,352,826	0	G	A	98.7%	248.4 / ‑4.3	74	P67L (CCG→CTG)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (41/32);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

CTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTACCACAAT  >  NC_000913/3352729‑3352920                                                                                                  |                                                                                               ctcgcgcATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGc                                                                                              >  2:1052595/1‑100 (MQ=255)   cgcgcATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGt                                                                                            <  1:680405/100‑1 (MQ=255)    gcgcATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTa                                                                                           >  2:14197/1‑100 (MQ=255)     cgcATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTaa                                                                                          >  2:662127/1‑100 (MQ=255)          ctcctcCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGc                                                                                     <  2:1091952/100‑1 (MQ=255)            cctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCaa                                                                                   >  2:921409/1‑100 (MQ=255)            cctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCaa                                                                                   >  1:302452/1‑100 (MQ=255)            cctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGCCGACCGCCCACAGCGTAATCAGCaa                                                                                   >  1:1464587/1‑100 (MQ=255)            cctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGACAGAAAGTAGACCGCCCACAGCGTAATCAGCaa                                                                                   >  2:419540/1‑100 (MQ=255)             ctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCaaa                                                                                  >  1:1629005/1‑100 (MQ=255)              tccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCa                                                                                    <  1:1507610/97‑1 (MQ=255)              tccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAAc                                                                                 >  2:1574536/1‑100 (MQ=255)                cAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACGGCGTAATCAGCAAAcc                                                                                <  1:121041/99‑1 (MQ=255)                cAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCa                                                                               >  1:1583890/1‑100 (MQ=255)                  gTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCaaa                                                                             <  2:722191/100‑1 (MQ=255)                     tttGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAaagaa                                                                          <  1:1491663/100‑1 (MQ=255)                      ttGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAaagaag                                                                         <  1:308202/100‑1 (MQ=255)                      ttGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAaagaag                                                                         >  2:715099/1‑100 (MQ=255)                       tGCACCAGCCGTAACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAagaaga                                                                        >  2:1115425/1‑100 (MQ=255)                        gCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGAt                                                                       <  2:1239067/100‑1 (MQ=255)                         cACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATa                                                                      >  2:1494027/1‑100 (MQ=255)                           ccAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGCCCACCGATAGAAAGTAGACCCCCCACAGCGTAATCAGCAAACCAAAGAAGATa                                                                      >  2:42364/1‑98 (MQ=255)                           ccAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATa                                                                      <  1:281227/98‑1 (MQ=255)                           ccAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGa                                                                    >  1:1360250/1‑100 (MQ=255)                           ccAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGa                                                                    >  1:1062093/1‑100 (MQ=255)                           ccAGCCGTGACAGACGTTGTCGTGACGCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACAAAAAAAGATAGa                                                                    >  2:111998/1‑100 (MQ=255)                             aGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAAc                                                                  >  2:1208146/1‑100 (MQ=255)                                cGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGaa                                                               <  1:188300/100‑1 (MQ=255)                                cGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGaa                                                               >  2:975446/1‑100 (MQ=255)                                 gTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAAt                                                              <  1:921407/100‑1 (MQ=255)                                 gTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAAt                                                              <  2:972184/100‑1 (MQ=255)                                 gTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAAt                                                              <  2:527137/100‑1 (MQ=255)                                  tGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATa                                                             <  1:892778/100‑1 (MQ=255)                                       gACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATc                                                        <  1:14197/100‑1 (MQ=255)                                         cGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCaa                                                      >  2:891210/1‑100 (MQ=255)                                         cGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCaa                                                      <  1:1389696/100‑1 (MQ=255)                                           ttGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATg                                                    <  1:547312/100‑1 (MQ=255)                                                gTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTc                                               <  2:1109918/100‑1 (MQ=255)                                                 tGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTc                                               >  1:1177267/1‑99 (MQ=255)                                                 tGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCg                                              <  2:157904/100‑1 (MQ=255)                                                  gACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCg                                              <  2:283172/99‑1 (MQ=255)                                                   aCTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGAc                                            <  1:1566610/100‑1 (MQ=255)                                                    ctcctcCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGAc                                            >  1:39750/1‑99 (MQ=255)                                                     tcttttaGTTGCTCGATGACTAGCGATAGAAAGTAGTCTGTGCACAGAGTAATTAGCATATTAATGAAGATAGAACGTATAACATCAATGCTTTCg                                              <  1:1491662/90‑1 (MQ=255)                                                      cctccAGTTGCTCGACGACCACCAATAGAAAATAGACCGCCCACAGCGTAATAAGCAAACCAAAGAAGATAGAACGAAAAACATCAATACTTTCGa                                             >  1:151428/1‑96 (MQ=255)                                                       ctccAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGa                                        >  2:1035390/1‑100 (MQ=255)                                                       ctccAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGa                                        >  2:911818/1‑100 (MQ=255)                                                         ccAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGAcc                                      >  1:1257095/1‑100 (MQ=255)                                                         ccAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGAcc                                      >  1:1058623/1‑100 (MQ=255)                                                         ccAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGAcc                                      >  2:1105534/1‑100 (MQ=255)                                                           aGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATCACATCAATGCTTTCCACCTGACCAt                                    >  1:362025/1‑100 (MQ=255)                                                           aGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTCTCGACCTGACCa                                     >  1:1396784/1‑99 (MQ=255)                                                            gTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATg                                   <  1:111998/100‑1 (MQ=255)                                                               gCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATgcag                                >  2:75045/1‑100 (MQ=255)                                                               gCTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATgcag                                <  1:1199921/100‑1 (MQ=255)                                                                cTCGACGACCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACAAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATgcagc                               >  1:1087610/1‑100 (MQ=255)                                                                  cgacgaCCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCAc                             <  2:1456880/100‑1 (MQ=255)                                                                  cgacgaCCACCGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCAc                             >  2:10567/1‑100 (MQ=255)                                                                       accaccGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATgg                        <  1:1511861/100‑1 (MQ=255)                                                                       accaccGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATgg                        <  1:10567/100‑1 (MQ=255)                                                                       accaccGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATgg                        <  2:496549/100‑1 (MQ=255)                                                                       accaccGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATgg                        <  2:302453/100‑1 (MQ=255)                                                                           ccGATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTAc                    >  1:371010/1‑100 (MQ=255)                                                                             gATAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCg                  <  1:1671629/100‑1 (MQ=255)                                                                              aTAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGc                 >  1:1217942/1‑100 (MQ=255)                                                                              aTAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGc                 >  2:422301/1‑100 (MQ=255)                                                                               tAGAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGcc                <  2:1062098/100‑1 (MQ=255)                                                                                 gAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCAt              >  2:843034/1‑100 (MQ=255)                                                                                 gAAAGTAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCAt              >  1:781657/1‑100 (MQ=255)                                                                                      tAGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTa         >  2:550092/1‑100 (MQ=255)                                                                                       aGACCGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTAc        <  1:959149/100‑1 (MQ=255)                                                                                          ccGCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACAGCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTACcac     <  2:102220/100‑1 (MQ=255)                                                                                            gCCCACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTACCACaa   >  1:584960/1‑100 (MQ=255)                                                                                             cccACAGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTACCACAAt  >  1:1261979/1‑100 (MQ=255)                                                                                                  |                                                                                               CTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGACCGCCCACGGCGTAATCAGCAAACCAAAGAAGATAGAACGAATAACATCAATGCTTTCGACCTGACCATGCAGCACCATGGTTACCGCCATTTGTACCACAAT  >  NC_000913/3352729‑3352920

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