Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,922,931	(A)7→6	intergenic (‑491/+126)	atpI ← / ← rsmG	ATP synthase accessory factor/16S rRNA m(7)G527 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,922,925	0	A	.	100.0%	201.3 / NA	47	intergenic (‑485/+132)	atpI/rsmG	ATP synthase accessory factor/16S rRNA m(7)G527 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (24/23);  total (24/23)

TACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTT  >  NC_000913/3922786‑3923047                                                                                                                                            |                                                                                                                           tACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCa                                                                                                                            >  1:45105/1‑140 (MQ=255) tACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCa                                                                                                                            >  2:314337/1‑140 (MQ=255)      ttttAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCaaaaaa                                                                                                                       >  1:812421/1‑140 (MQ=255)      ttttAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCaaaaaa                                                                                                                       >  1:660925/1‑140 (MQ=255)       tttAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCaaaaaag                                                                                                                      <  2:45105/140‑2 (MQ=255)        ttAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCaaaaaagc                                                                                                                     >  1:395799/1‑138 (MQ=255)          aataaATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCaaaaaagcca                                                                                                                   <  1:460098/140‑5 (MQ=255)           ataaATGACATTTACGCAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCAc                                                                                                                 >  2:39095/1‑140 (MQ=255)            taaATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACg                                                                                                                <  2:563706/140‑1 (MQ=255)              aaTGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGtt                                                                                                              >  1:572340/1‑140 (MQ=255)                tGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTAt                                                                                                            >  2:771913/1‑140 (MQ=255)                     tttACACAACAAAACCCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTTCTTACC‑AAAAAAGCCACGTTATCttgt                                                                                                       <  2:396828/140‑1 (MQ=255)                     tttACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCttgt                                                                                                       <  2:660925/140‑1 (MQ=255)                       tACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCttgttg                                                                                                     >  2:2790/1‑140 (MQ=255)                        acacAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGa                                                                                                    <  1:504087/140‑1 (MQ=255)                          acaacaAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATg                                                                                                  <  1:259300/140‑1 (MQ=255)                                 aaaCCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAaga                                                                                           >  2:813853/1‑140 (MQ=255)                                 aaaCCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAaga                                                                                           >  2:425856/1‑140 (MQ=255)                                      accCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGaa                                                                                      >  1:533080/1‑140 (MQ=255)                                      accCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGaa                                                                                      >  2:571972/1‑140 (MQ=255)                                      accCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATATTGTTGATGCAAAAGAGTGaa                                                                                      >  2:988197/1‑140 (MQ=255)                                       cccATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAAc                                                                                     <  1:813853/140‑1 (MQ=255)                                        ccATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACg                                                                                    >  1:864409/1‑140 (MQ=255)                                        ccATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACg                                                                                    >  2:712100/1‑140 (MQ=255)                                         cATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGt                                                                                   <  2:477658/140‑1 (MQ=255)                                              aCATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                              <  2:533080/140‑1 (MQ=255)                                              aCATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGt                                                                              <  1:425856/140‑1 (MQ=255)                                                  ttttAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTaaa                                                                          >  2:686504/1‑140 (MQ=255)                                                     taataaTGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGt                                                                       <  2:299797/140‑1 (MQ=255)                                                       ataatGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  1:712100/140‑1 (MQ=255)                                                       ataatGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  2:572340/140‑1 (MQ=255)                                                       ataatGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  1:686504/140‑1 (MQ=255)                                                       ataatGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  1:537977/140‑1 (MQ=255)                                                       ataatGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  1:39095/140‑1 (MQ=255)                                                          atGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCa                                                                  <  1:462760/140‑1 (MQ=255)                                                            gTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGt                                                                >  1:397717/1‑140 (MQ=255)                                                              tttAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     <  1:209351/133‑1 (MQ=255)                                                              tttAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTaa                                                                     >  2:209351/1‑133 (MQ=255)                                                                tAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTtata                                                            >  1:959293/1‑140 (MQ=255)                                                                  aCAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATc                                                          <  1:571697/140‑1 (MQ=255)                                                                        aaTGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGa                                                    <  1:710794/140‑1 (MQ=255)                                                                           gatgGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGaaaa                                                 <  2:42574/140‑1 (MQ=255)                                                                                ttCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTgg                                            >  2:330531/1‑140 (MQ=255)                                                                                ttCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTgg                                            <  1:282259/140‑1 (MQ=255)                                                                                ttCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTgg                                            >  1:828106/1‑140 (MQ=255)                                                                                 tCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgt                                           >  2:720287/1‑140 (MQ=255)                                                                                  cTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgtt                                          >  1:121011/1‑140 (MQ=255)                                                                                    tAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgt                                        <  2:121011/140‑1 (MQ=255)                                                                                      gcgcCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgttg                                      >  1:211059/1‑140 (MQ=255)                                                                                       cgcCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGgttgttgt                                     >  2:507000/1‑140 (MQ=255)                                                                                         ccGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTa                                   >  2:492952/1‑140 (MQ=255)                                                                                          cGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTaa                                  <  1:691119/140‑1 (MQ=255)                                                                                                                     tAATGACTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGAt       >  2:663239/1‑140 (MQ=255)                                                                                                                          aCTTATTTTCTGCTTACC‑AAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGAtttttt  >  2:234851/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                           TACTGTTTTAATAAATGACATTTACACAACAAAAACCACCCATTGACATTTTTAATAATGTTTTAACAGCCAATGATGGTTCTTAGCGCCGATTTTTAGCAGACTGATATTTTCACTAATGACTTATTTTCTGCTTACCAAAAAAAGCCACGTTATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTT  >  NC_000913/3922786‑3923047

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