Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,860,630	(C)5→4	coding (348/894 nt)	yidL →	putative DNA‑binding transcriptional regulator YidL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,860,626	0	C	.	100.0%	162.2 / NA	39	coding (344/894 nt)	yidL	putative DNA‑binding transcriptional regulator YidL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/18);  total (21/18)

CAAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGT  >  NC_000913/3860487‑3860758                                                                                                                                            |                                                                                                                                     cAAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTAc                                                                                                                                      <  2:128564/140‑1 (MQ=255)         tgGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAg                                                                                                                             >  2:370202/1‑140 (MQ=255)          gggAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGt                                                                                                                            <  2:501353/140‑1 (MQ=255)                aTTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTatgatg                                                                                                                      <  1:294601/140‑1 (MQ=255)                   cTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGat                                                                                                                   <  2:58581/140‑1 (MQ=255)                   cTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGat                                                                                                                   <  2:148959/140‑1 (MQ=255)                              tGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTgg                                                                                                        >  1:69446/1‑140 (MQ=255)                                   aaCAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGc                                                                                                   >  2:303319/1‑140 (MQ=255)                                    aCAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCa                                                                                                  <  1:303319/140‑1 (MQ=255)                                      aaaTAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCaaa                                                                                                <  1:313586/140‑1 (MQ=255)                                         tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                             <  2:69446/140‑1 (MQ=255)                                         tAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCg                                                                                             <  1:176961/140‑1 (MQ=255)                                                gCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGttattta                                                                                      >  1:249579/1‑140 (MQ=255)                                                        aCTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCg                                                                              <  1:370202/140‑1 (MQ=255)                                                                  atTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATaa                                                                    >  2:508640/1‑140 (MQ=255)                                                                      ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGTATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAg                                                                >  1:434408/1‑140 (MQ=255)                                                                      ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGTGAAATTTATAATCAg                                                                >  2:453918/1‑140 (MQ=255)                                                                      ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAg                                                                >  2:218644/1‑140 (MQ=255)                                                                      ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAg                                                                >  2:440079/1‑140 (MQ=255)                                                                      ttAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAg                                                                <  2:240263/140‑1 (MQ=255)                                                                           gCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAgg                                                               >  2:376870/1‑136 (MQ=255)                                                                           gCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAgg                                                               <  1:376870/136‑1 (MQ=255)                                                                           gCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACt                                                           >  2:364262/1‑140 (MQ=255)                                                                            ccAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTc                                                          <  2:504924/140‑1 (MQ=255)                                                                              aaTGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCAc                                                        >  2:475219/1‑140 (MQ=255)                                                                               aTGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCAcc                                                       >  2:559254/1‑140 (MQ=255)                                                                                      ttCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                                                >  1:196274/1‑140 (MQ=255)                                                                                      ttCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTg                                                >  1:48681/1‑140 (MQ=255)                                                                                        cACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGTATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCt                                              <  2:434408/140‑1 (MQ=255)                                                                                           tctTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAg                                           >  1:380813/1‑140 (MQ=255)                                                                                               aTCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaa                                       >  1:198343/1‑140 (MQ=255)                                                                                                  aCTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTtaataa                                    <  1:70809/140‑1 (MQ=255)                                                                                                   cTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAAc                                   <  2:380813/140‑1 (MQ=255)                                                                                                              tAGTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAg                        >  2:429265/1‑140 (MQ=255)                                                                                                                gTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAAc                                   <  2:503170/127‑1 (MQ=255)                                                                                                                gTCTGGGAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAAc                                   >  1:503170/1‑127 (MQ=255)                                                                                                                      gAACAGTACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATaaaa                <  1:249100/140‑1 (MQ=255)                                                                                                                            tACTGGATGGAATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAaataat          <  1:364262/140‑1 (MQ=255)                                                                                                                                   tggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAg   >  1:107718/1‑140 (MQ=255)                                                                                                                                    ggaATTTA‑CCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGt  >  1:311377/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                     CAAGGATGTGGGAAAATTCTTTTAAAGAATGGCGAACAAATAACACTGCATGGCAACTGTATTATATTCTTAAAGCCAATGGATATTCACTCTTATCACTGTGAAGGTTTAGTCTGGGAACAGTACTGGATGGAATTTACCCCCACCAGTATGATGGATATTCCCGTTGGTCAGCAAAGCGTTATTTATAATGGCGAAATTTATAATCAGGAACTCACCGAAGTTGCTGAGTTAATAACTTCACCAGAAGCAATAAAAAATAATCTGGCAGT  >  NC_000913/3860487‑3860758

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