Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,922,931	(A)7→8	intergenic (‑491/+126)	atpI ← / ← rsmG	ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent, glucose‑inhibited cell‑division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,922,924	1	.	A	100.0%	155.6 / NA	48	intergenic (‑484/+133)	atpI/rsmG	ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent, glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (22/26);  total (22/26)

TACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTAT  >  NC_000913/3922808‑3923039                                                                                                                      |                                                                                                                    tACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttg                                                                                                >  2:797072/1‑139 (MQ=255)  acacAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgt                                                                                               <  1:217813/139‑1 (MQ=255)    acaacaAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgttg                                                                                             >  2:956308/1‑139 (MQ=255)                  ccATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGa                                                                               >  1:828762/1‑139 (MQ=255)                     ttGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACg                                                                            <  1:797072/139‑1 (MQ=255)                                aataatGTTTTGACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCTCGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAAt                                                                 <  2:790858/139‑1 (MQ=255)                                    atGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTAt                                                                                                    <  2:624026/100‑1 (MQ=255)                                    atGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTAt                                                                                                    >  1:624026/1‑100 (MQ=255)                                    atGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                             <  1:550928/139‑1 (MQ=255)                                      gTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAAcc                                                           >  2:277721/1‑139 (MQ=255)                                       ttttAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCa                                                          >  2:1104396/1‑139 (MQ=255)                                              aGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTatat                                                   <  1:277721/139‑1 (MQ=255)                                               gCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTTATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTa                                                              <  2:750764/127‑1 (MQ=255)                                               gCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                      >  1:415715/1‑119 (MQ=255)                                               gCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                      <  2:415715/119‑1 (MQ=255)                                               gCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTa                                                              >  1:750764/1‑127 (MQ=255)                                                 cAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGt                                                                           >  2:499625/1‑112 (MQ=255)                                                 cAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGt                                                                           <  1:499625/112‑1 (MQ=255)                                                 cAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTatat                                                   >  1:1039693/1‑136 (MQ=255)                                                 cAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTatat                                                   <  2:1039693/136‑1 (MQ=255)                                                  aaTGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgt                                                                                               >  2:587602/1‑91 (MQ=255)                                                  aaTGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgt                                                                                               <  1:587602/91‑1 (MQ=255)                                                   atgatgGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTa                                              <  2:187040/139‑1 (MQ=255)                                                             ttAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTgg                                    <  1:946202/139‑1 (MQ=255)                                                             ttAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTgg                                    <  1:956308/139‑1 (MQ=255)                                                               aGCGCCGATTTTTAGCTGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTCgtt                                  <  2:828762/139‑1 (MQ=255)                                                               aGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGa                                                                                            <  1:882679/81‑1 (MQ=255)                                                               aGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGa                                                                                            >  2:882679/1‑81 (MQ=255)                                                               aGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCa                                                                                        <  1:12347/85‑1 (MQ=255)                                                               aGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCa                                                                                        >  2:12347/1‑85 (MQ=255)                                                               aGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGCAGAAAACCTGgtt                                  >  2:673244/1‑139 (MQ=255)                                                                 cgcCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGt                                                               <  2:1087073/108‑1 (MQ=255)                                                                 cgcCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGt                                                               >  1:1087073/1‑108 (MQ=255)                                                                    cGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgttgt                             >  1:495819/1‑139 (MQ=255)                                                                     gATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGCAGAAAACCTGgttgttgtt                            <  1:673244/139‑1 (MQ=255)                                                                       tttttAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTaa                          <  2:495819/139‑1 (MQ=255)                                                                               aGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTaa                  <  1:173707/139‑1 (MQ=255)                                                                                gACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgttg                                                                                             >  1:903122/1‑63 (MQ=255)                                                                                gACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCttgttg                                                                                             <  2:903122/63‑1 (MQ=255)                                                                                           tCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                             <  2:1078533/84‑1 (MQ=255)                                                                                           tCACTAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                             >  1:1078533/1‑84 (MQ=255)                                                                                               tAATGACTTATTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTAt  >  1:1216126/1‑139 (MQ=255)                                                                                                        aTTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTa                                                              >  2:979040/1‑70 (MQ=255)                                                                                                        aTTTTCTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTa                                                              <  1:979040/70‑1 (MQ=255)                                                                                                             cTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgttgtt                            <  2:769212/99‑1 (MQ=255)                                                                                                             cTGCTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgttgtt                            >  1:769212/1‑99 (MQ=255)                                                                                                                cTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                      <  2:165863/54‑1 (MQ=255)                                                                                                                cTTACCAAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                      >  1:165863/1‑54 (MQ=255)                                                                                                                      |                                                                                                                    TACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTAT  >  NC_000913/3922808‑3923039

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