Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,181,024	A→G	intergenic (+99/‑221)	rplL → / → rpoB	50S ribosomal subunit protein L12/RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,181,024	0	A	G	100.0%	195.9 / NA	57	intergenic (+99/‑221)	rplL/rpoB	50S ribosomal subunit protein L12/RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/32);  total (25/32)

TCTGGAAGAAGCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTC  >  NC_000913/4180889‑4181159                                                                                                                                        |                                                                                                                                        tctGGAAGAAGCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGt                                                                                                                                     >  2:1466990/1‑140 (MQ=255) tctGGAAGAAGCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGt                                                                                                                                     <  2:241145/140‑1 (MQ=255)      aagaagCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGac                                                                                                                                <  1:1292894/140‑1 (MQ=255)           gCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTActg                                                                                                                           <  2:22077/140‑1 (MQ=255)             tGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTActgct                                                                                                                         >  1:341360/1‑140 (MQ=255)               gcgcTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCtgt                                                                                                                       <  2:1326650/140‑1 (MQ=255)               gcgcTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCtgt                                                                                                                       >  2:1290990/1‑140 (MQ=255)                    gaagttgaagttAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCtt                                                                                                                  >  2:1412371/1‑140 (MQ=255)                          gaagttAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATg                                                                                                            <  1:587270/140‑1 (MQ=255)                                aaaTAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGtt                                                                                                      >  1:1160396/1‑140 (MQ=255)                                 aaTAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGttt                                                                                                     >  2:643884/1‑140 (MQ=255)                                  aTAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTc                                                                                                    >  2:64764/1‑140 (MQ=255)                                      gCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATc                                                                                                >  2:159203/1‑140 (MQ=255)                                        cAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              <  1:1268456/140‑1 (MQ=255)                                        cAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              <  1:956181/140‑1 (MQ=255)                                        cAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              <  1:702998/140‑1 (MQ=255)                                        cAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              <  1:712669/140‑1 (MQ=255)                                         aaCCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcgac                                                                                             <  1:1101301/140‑1 (MQ=255)                                          aCCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcgacg                                                                                            <  2:1160396/140‑1 (MQ=255)                                               tCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTa                                                                                       <  2:1111772/140‑1 (MQ=255)                                               tCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTa                                                                                       >  2:831068/1‑140 (MQ=255)                                               tCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTa                                                                                       >  2:694529/1‑140 (MQ=255)                                                ccGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTaa                                                                                      >  1:744970/1‑140 (MQ=255)                                                 cGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              >  2:608288/1‑131 (MQ=255)                                                 cGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATcga                                                                                              <  1:608288/131‑1 (MQ=255)                                                    ttGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAAtata                                                                                  <  2:341360/140‑1 (MQ=255)                                                    ttGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAAtata                                                                                  <  2:593875/140‑1 (MQ=255)                                                      gCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATtgacg                                                                                            <  2:214794/128‑1 (MQ=255)                                                      gCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATtgacg                                                                                            >  1:214794/1‑128 (MQ=255)                                                      gCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACt                                                                                >  2:316200/1‑140 (MQ=255)                                                        aGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGc                                                                              >  2:457395/1‑140 (MQ=255)                                                             gagaAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAg                                                                         <  1:1466990/140‑1 (MQ=255)                                                              agaAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACt                                                                                <  1:600492/132‑1 (MQ=255)                                                              agaAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACt                                                                                >  2:600492/1‑132 (MQ=255)                                                              agaAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAgg                                                                        >  2:956742/1‑140 (MQ=255)                                                                 aaTCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACg                                                                     <  1:831068/140‑1 (MQ=255)                                                                  aTCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGt                                                                    <  1:948777/140‑1 (MQ=255)                                                                    cAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTcc                                                                  >  2:931954/1‑140 (MQ=255)                                                                       gCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGtt                                                               <  1:227224/140‑1 (MQ=255)                                                                           aTGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCtgt                                                           >  1:310528/1‑140 (MQ=255)                                                                                          ttAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaa                                            <  2:299150/140‑1 (MQ=255)                                                                                           tAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaaa                                           <  1:931954/140‑1 (MQ=255)                                                                                           tAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaaa                                           <  1:711440/140‑1 (MQ=255)                                                                                           tAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaaa                                           <  1:316200/140‑1 (MQ=255)                                                                                           tAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATACTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGaaa                                           <  1:1290990/140‑1 (MQ=255)                                                                                                 ccAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGttt                                     >  2:1075473/1‑140 (MQ=255)                                                                                                  cAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTa                                    <  1:159203/140‑1 (MQ=255)                                                                                                    gCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAg                                  <  1:1412371/140‑1 (MQ=255)                                                                                                       tttttGCGCTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCACAATGAAATGGTTTAAGCGt                               <  2:744970/140‑1 (MQ=255)                                                                                                               cTGTAAGGCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                       >  1:1377569/1‑140 (MQ=255)                                                                                                                     ggCGCCAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCa                 <  2:1377569/140‑1 (MQ=255)                                                                                                                          cAGTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCg            >  1:658004/1‑140 (MQ=255)                                                                                                                            gTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGa          >  1:463127/1‑140 (MQ=255)                                                                                                                            gTAGCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGa          <  2:310528/140‑1 (MQ=255)                                                                                                                                 gTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGTGGAAAgtg     <  1:64764/140‑1 (MQ=255)                                                                                                                                   ttCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGtt   >  2:562778/1‑140 (MQ=255)                                                                                                                                    tCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc  >  2:1038360/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        TCTGGAAGAAGCTGGCGCTGAAGTTGAAGTTAAATAAGCCAACCCTTCCGGTTGCAGCCTGAGAAATCAGGCTGATGGCTGGTGACTTTTTAGTCACCAGCCTTTTTGCGCTGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTC  >  NC_000913/4180889‑4181159

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: ‑