Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,622,422	C→G	A1077A (GCC→GCG)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,622,422	0	C	G	100.0%	172.4 / NA	55	A1077A (GCC→GCG)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (33/22);  total (33/22)

GCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAA  >  NC_000913/3622290‑3622506                                                                                                                                     |                                                                                     gCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGGGCTTCACGCCACTCATCAGGGCCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTcagc                                                                               >  1:1467428/1‑140 (MQ=255) gCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTcagc                                                                               >  1:945563/1‑140 (MQ=255)     aCCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTc                                                                           >  2:1118077/1‑140 (MQ=255)       cACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCcg                                                                         <  2:724716/140‑1 (MQ=255)            tACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCggcggc                                                                    >  2:891554/1‑140 (MQ=255)            tACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCggcggc                                                                    <  1:1598898/140‑1 (MQ=255)            tACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCggcggc                                                                    >  1:150287/1‑140 (MQ=255)             aCAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGgcggcg                                                                   <  1:459947/140‑1 (MQ=255)               aGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGaa                                                                 >  1:193462/1‑140 (MQ=255)                 aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAGACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGa                                                               >  1:1480795/1‑140 (MQ=255)                 aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGa                                                               >  2:1430630/1‑140 (MQ=255)                 aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGa                                                               >  1:249715/1‑140 (MQ=255)                 aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGa                                                               >  1:850626/1‑140 (MQ=255)                 aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGa                                                               >  2:1117455/1‑140 (MQ=255)                   cGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACg                                                             >  1:1533161/1‑140 (MQ=255)                     aCAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGt                                                           <  1:1388875/140‑1 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCTGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          <  2:695729/140‑1 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          >  2:1263227/1‑140 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          >  1:972511/1‑140 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          <  2:972511/140‑1 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          >  2:1695252/1‑140 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          >  2:311146/1‑140 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          <  2:413686/140‑1 (MQ=255)                      cAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTg                                                          >  1:1164531/1‑140 (MQ=255)                         gAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCa                                                       <  1:891554/140‑1 (MQ=37)                         gAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCa                                                       >  1:1120936/1‑140 (MQ=37)                         gAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCa                                                       <  2:1096196/140‑1 (MQ=37)                           gCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgt                                                     >  1:1398125/1‑140 (MQ=25)                           gCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgt                                                     >  1:1242350/1‑140 (MQ=25)                              gCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                  <  2:787565/140‑1 (MQ=21)                              gCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                  <  1:1117455/140‑1 (MQ=21)                              gCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                  <  1:175737/140‑1 (MQ=21)                              gCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                  <  1:348867/140‑1 (MQ=21)                                aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTgtggt                                                >  2:1005232/1‑140 (MQ=21)                                     gACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTcc                                           >  2:960797/1‑140 (MQ=21)                                         aTTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccgcc                                       >  2:440141/1‑140 (MQ=25)                                          tttATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccgccg                                      >  2:1466178/1‑140 (MQ=25)                                             aTCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTg                                   <  2:1164531/140‑1 (MQ=21)                                             aTCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCAAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGGGGAAGTGTGGTGTCCCCGCCGGTg                                   <  2:1303032/140‑1 (MQ=11)                                                  ccGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGt                              >  2:1054603/1‑140 (MQ=25)                                                  ccGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGt                              >  1:1381358/1‑140 (MQ=25)                                                  ccGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGt                              >  2:709580/1‑140 (MQ=37)                                                    ggggAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTATGGTGTTCCCGCCGGTGCTGGTGc                            <  1:311146/140‑1 (MQ=21)                                                           ttCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCt                     <  2:1533161/140‑1 (MQ=21)                                                            tCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTg                                   <  2:1627829/125‑1 (MQ=25)                                                            tCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTg                                   >  1:1627829/1‑125 (MQ=38)                                                            tCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTc                    >  2:1593973/1‑140 (MQ=37)                                                               cGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGAc                 <  1:1590048/140‑1 (MQ=21)                                                                  cACTCATCAGGGTCGAAACTGCCACCGGTGAGCGGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCgg              <  2:1687685/140‑1 (MQ=21)                                                                       atcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa         <  2:1539613/140‑1 (MQ=21)                                                                       atcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa         <  1:440141/140‑1 (MQ=21)                                                                       atcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGa         >  1:554081/1‑140 (MQ=37)                                                                        tcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaa        >  2:1536709/1‑140 (MQ=37)                                                                           gggTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGt     <  1:1263227/140‑1 (MQ=18)                                                                              tCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGTCTGGCGGATGCGCTTCAGGAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGaa  >  2:1512686/1‑140 (MQ=9)                                                                              tCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGc             >  2:1099810/1‑129 (MQ=17)                                                                              tCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGc             <  1:1099810/129‑1 (MQ=14)                                                                                                                                     |                                                                                     GCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAA  >  NC_000913/3622290‑3622506

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