Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,073	(A)6→5	coding (1296/1506 nt)	rbsA →	ribose ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,068	0	A	.	100.0%	215.6 / NA	51	coding (1291/1506 nt)	rbsA	ribose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/28);  total (23/28)

AACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGAAAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCATCTCAGCGGG  >  NC_000913/3934932‑3935202                                                                                                                                         |                                                                                                                                       aaCTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaa                                                                                                                                     <  1:95596/140‑1 (MQ=255)  aCTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaaa                                                                                                                                    <  2:27693/140‑1 (MQ=255)    tCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaaaga                                                                                                                                  <  1:644091/140‑3 (MQ=255)    tCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaaaga                                                                                                                                  <  2:423011/140‑3 (MQ=255)     ccGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaaagag                                                                                                                                 >  2:619110/1‑137 (MQ=255)      cGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGaaaaagaga                                                                                                                                >  1:865983/1‑136 (MQ=255)       gTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGAt                                                                                                                              >  2:4720/1‑140 (MQ=255)       gTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGAt                                                                                                                              <  1:371503/140‑1 (MQ=255)            tGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATc                                                                                                                         >  2:917821/1‑140 (MQ=255)                       aTTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAAc                                                                                                              <  2:865983/140‑1 (MQ=255)                        ttGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGACTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAAcc                                                                                                             >  2:892237/1‑140 (MQ=255)                         tGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATc                                                                                                                             <  1:516072/123‑1 (MQ=255)                         tGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATc                                                                                                                             >  2:516072/1‑123 (MQ=255)                         tGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCa                                                                                                            <  1:620644/140‑1 (MQ=255)                              tGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCa                                                                                                       >  1:779401/1‑140 (MQ=255)                               gCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCaa                                                                                                      <  2:865063/140‑1 (MQ=255)                                         ggCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGc                                                                                            <  2:1384/140‑1 (MQ=255)                                          gCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCt                                                                                           >  2:188386/1‑140 (MQ=255)                                          gCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCt                                                                                           >  1:673891/1‑140 (MQ=255)                                          gCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCt                                                                                           >  1:260810/1‑140 (MQ=255)                                           cAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCtt                                                                                          <  2:260810/140‑1 (MQ=255)                                                aGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGca                                                                                     <  1:540783/140‑1 (MQ=255)                                                aGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGca                                                                                     <  1:619110/140‑1 (MQ=255)                                                     aaaGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCAtt                                                                                <  2:779401/140‑1 (MQ=255)                                                        gTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTg                                                                             >  1:267602/1‑140 (MQ=255)                                                                   ccGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGga                                                                  <  1:709083/140‑1 (MQ=255)                                                                    cgcgGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGgag                                                                 <  2:548069/140‑1 (MQ=255)                                                                          cTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCa                                                           <  2:267602/140‑1 (MQ=255)                                                                           tgatgaCACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAg                                                          <  1:409572/140‑1 (MQ=255)                                                                           tgatgaCACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAg                                                          <  1:892237/140‑1 (MQ=255)                                                                            gatgaCACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGa                                                         <  1:304628/140‑1 (MQ=255)                                                                             atgaCACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGaa                                                        >  1:41498/1‑140 (MQ=255)                                                                                acacGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTa                                                     <  1:188386/140‑1 (MQ=255)                                                                                  acGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGga                                                                  <  2:5695/125‑1 (MQ=255)                                                                                  acGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGga                                                                  >  1:5695/1‑125 (MQ=255)                                                                                   cGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTa                                                  >  1:288390/1‑140 (MQ=255)                                                                                     ccccAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAg                                                          >  2:781670/1‑130 (MQ=255)                                                                                     ccccAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAg                                                          <  1:781670/130‑1 (MQ=255)                                                                                     ccccAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAgg                                                >  2:911256/1‑140 (MQ=255)                                                                                       ccAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCa                                              >  2:242337/1‑140 (MQ=255)                                                                                          aaGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGa                                           <  1:242337/140‑1 (MQ=255)                                                                                            gtgtTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGc                                         <  2:673891/140‑1 (MQ=255)                                                                                                tGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATc                                     >  2:209036/1‑140 (MQ=255)                                                                                                  aTCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGc                                   >  2:568188/1‑140 (MQ=255)                                                                                                            aaCCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTca                         <  2:874089/140‑1 (MQ=255)                                                                                                            aaCCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTca                         >  1:450029/1‑140 (MQ=255)                                                                                                             aCCTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTcat                        <  2:41498/140‑1 (MQ=255)                                                                                                              ccTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTcatg                       >  2:879967/1‑140 (MQ=255)                                                                                                              ccTACCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTcatg                       >  1:360831/1‑140 (MQ=255)                                                                                                                 aCCCGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTcatgcat                    >  1:635751/1‑140 (MQ=255)                                                                                                                   ccGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGa                  >  2:46779/1‑140 (MQ=255)                                                                                                                    cGTGGCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGaa                 <  2:501045/140‑1 (MQ=255)                                                                                                                        gCGTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGc             >  1:349217/1‑140 (MQ=255)                                                                                                                          gTAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCAt           <  1:917821/140‑1 (MQ=255)                                                                                                                           tAGATGTCGGCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCAtc          <  1:911256/140‑1 (MQ=255)                                                                                                                                   ggCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCATCTCAGCggg  <  1:879967/140‑1 (MQ=255)                                                                                                                                   ggCGCG‑AAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCATCTCAGCggg  <  1:4720/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       AACTCCGTCGATGGAACAGGCAATTGGTCTGCTTTCCGGTGGCAATCAGCAAAAAGTGGCGATTGCCCGCGGTCTGATGACACGCCCCAAAGTGTTGATCCTTGATGAACCTACCCGTGGCGTAGATGTCGGCGCGAAAAAAGAGATCTATCAACTGATTAACCAGTTCAAAGCCGATGGCTTGAGCATCATTCTGGTGTCATCGGAGATGCCAGAAGTATTAGGCATGAGCGATCGCATCATCGTCATGCATGAAGGGCATCTCAGCGGG  >  NC_000913/3934932‑3935202

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