Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	476,808	T→C	V46A (GTA→GCA)	ybaA →	DUF1428 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	476,808	0	T	C	100.0%	204.1 / NA	60	V46A (GTA→GCA)	ybaA	DUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/32);  total (28/32)

GATAAAAAGGATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGTACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAAGAGG  >  NC_000913/476711‑476903                                                                                                  |                                                                                                gATAAAAAGGATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCAc                                                                                               >  2:1008675/1‑100 (MQ=255) gATAAAAAGGATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCAc                                                                                               >  2:1484937/1‑100 (MQ=255)        aGGATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCAGAt                                                                                          >  2:800444/1‑98 (MQ=255)          gATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGc                                                                                       <  1:1521833/99‑1 (MQ=255)             gCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAg                                                                                   <  2:78840/100‑1 (MQ=255)               cTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTg                                                                                 <  1:757254/100‑1 (MQ=255)               cTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTg                                                                                 >  2:1368222/1‑100 (MQ=255)                tATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGa                                                                                >  1:1596503/1‑100 (MQ=255)                tATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGa                                                                                >  2:114717/1‑100 (MQ=255)                tATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGa                                                                                <  2:1436800/100‑1 (MQ=255)                tATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGa                                                                                <  1:1008665/100‑1 (MQ=255)                tATCGGGAAATGGACGCTAAGGCTGCGCCATAGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGa                                                                                >  2:1223067/1‑100 (MQ=255)                  tCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGAcc                                                                              <  1:1376045/100‑1 (MQ=255)                  tCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGAcc                                                                              >  2:228116/1‑100 (MQ=255)                    gggAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGa                                                                            <  2:1552472/100‑1 (MQ=255)                     ggAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGAt                                                                           <  2:780709/100‑1 (MQ=255)                      gAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGTGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGAt                                                                           <  1:1034942/99‑1 (MQ=255)                         aTGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGAtttt                                                                        >  1:411697/1‑99 (MQ=255)                           ggCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGt                                                                     <  2:1427454/100‑1 (MQ=255)                           ggCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGt                                                                     >  1:1497921/1‑100 (MQ=255)                            gCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTa                                                                    <  2:1134686/100‑1 (MQ=255)                            gCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTa                                                                    <  1:469669/100‑1 (MQ=255)                             ccGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGACAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTAt                                                                   >  1:2133/1‑100 (MQ=255)                               gCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATgg                                                                 >  1:1306302/1‑100 (MQ=255)                               gCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATgg                                                                 <  2:1684737/100‑1 (MQ=255)                                cTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGc                                                                <  1:204933/100‑1 (MQ=255)                                    ggCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTg                                                            >  2:541990/1‑100 (MQ=255)                                    ggCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTg                                                            <  2:1755485/100‑1 (MQ=255)                                     gCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGa                                                           <  2:716664/100‑1 (MQ=255)                                         cgcCATTGTTTAAAGAGTTTTGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGc                                                       <  2:304622/100‑1 (MQ=255)                                             aTTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGaa                                                   >  2:934901/1‑100 (MQ=255)                                              ttGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAag                                                  <  1:782833/100‑1 (MQ=255)                                              ttGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAag                                                  <  2:1666449/100‑1 (MQ=255)                                                gTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAagag                                                <  1:74575/100‑1 (MQ=255)                                                   tAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAAt                                             >  2:642279/1‑100 (MQ=255)                                                    aaaGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATg                                            >  2:993503/1‑100 (MQ=255)                                                      agagTTTGGCGCGCTACGTATTGTCAAATGATGGGCAAGCGATGCACCGGATGGCAAAGTGAACGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGaa                                          >  1:9550/1‑100 (MQ=255)                                                        agTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAaga                                        >  2:1111646/1‑100 (MQ=255)                                                         gTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAagag                                       <  2:87372/100‑1 (MQ=255)                                                         gTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAagag                                       <  2:868543/100‑1 (MQ=255)                                                          tttGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAgg                                      <  1:1787961/100‑1 (MQ=255)                                                               cgcgcTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTc                                 >  1:1113746/1‑100 (MQ=255)                                                                 cgcTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCtt                               <  2:173928/100‑1 (MQ=255)                                                                   cTTCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTa                             >  1:275208/1‑100 (MQ=255)                                                                    ttCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCt                                <  1:1461522/96‑1 (MQ=255)                                                                     tCGTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGc                           <  2:575918/100‑1 (MQ=255)                                                                       gTATTGTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTg                         <  2:561520/100‑1 (MQ=255)                                                                            gTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTg                    >  2:360943/1‑100 (MQ=255)                                                                            gTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTg                    <  2:442968/100‑1 (MQ=255)                                                                            gTCGAATGCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTg                    >  2:271843/1‑100 (MQ=255)                                                                                   gCTGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATAccc             >  2:414191/1‑100 (MQ=255)                                                                                     tGGGCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCtt           >  1:514975/1‑100 (MQ=255)                                                                                      gggCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTc          <  2:305343/100‑1 (MQ=255)                                                                                       ggCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCa         >  1:830783/1‑100 (MQ=255)                                                                                       ggCCAGCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCa         >  1:1748371/1‑100 (MQ=255)                                                                                            gCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAaga    >  2:1799663/1‑100 (MQ=255)                                                                                            gCGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAaga    <  2:1497932/100‑1 (MQ=255)                                                                                             cGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAagag   <  2:1140205/100‑1 (MQ=255)                                                                                             cGATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAagag   >  2:1483347/1‑100 (MQ=255)                                                                                              gATGCACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAAGAgg  <  1:360941/100‑1 (MQ=255)                                                                                                  |                                                                                                GATAAAAAGGATGCCTATCGGGAAATGGCCGCTAAGGCTGCGCCATTGTTTAAAGAGTTTGGCGCGCTTCGTATTGTCGAATGCTGGGCCAGCGATGTACCGGATGGCAAAGTGACCGATTTTCGTATGGCGGTGAAAGCGGAAGAGAATGAAGAGGTTGTCTTTAGCTGGATTGAATACCCTTCAAAAGAGG  >  NC_000913/476711‑476903

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