Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,854	G→A	G1224G (GGG→GGA)	rhsA →	rhs element protein RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,854	0	G	A	100.0%	134.2 / NA	43	G1224G (GGG→GGA)	rhsA	rhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/14);  total (29/14)

CGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGC  >  NC_000913/3765757‑3765988                                                                                                  |                                                                                                                                       cGGGGCAGCAGTATGATCAGGAGTCCGGCCCGTCTAACATCCGCCACCCCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCa                                                                                              <  2:412596/140‑1 (MQ=11)   gggCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGt                                                                                            <  2:525902/140‑1 (MQ=16)       agcagTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACg                                                                                        <  2:807137/140‑1 (MQ=25)             atgatgAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACa                                                                                  >  2:20752/1‑140 (MQ=255)              tgatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAGTTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAg                                                                                 <  2:496849/140‑1 (MQ=255)                   aggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCTCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATACGTGGAATCCAGTTACGAATACAGATCct                                                                            >  1:379123/1‑140 (MQ=14)                   aggagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGCAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCct                                                                            >  2:589142/1‑140 (MQ=255)                     gagTCCGGCCTGTATGACAACCTCCGCCGTTATTATGGCCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCctct                                                                          <  1:404812/140‑1 (MQ=14)                        tCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTggg                                                                       >  2:517398/1‑140 (MQ=255)                               tGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTCTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGaa                                                                >  2:439028/1‑140 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                                                             <  2:199905/140‑1 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                                                             >  2:479686/1‑140 (MQ=255)                                  aTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGtt                                                             >  1:437609/1‑140 (MQ=255)                                      cAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATTCAGTTACAAATACAGACCCTCCGGAGTTAGAAGTTTTTc                                                         >  2:323019/1‑140 (MQ=17)                                           gCCACCGCTATTATGACCTGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACATATCTTCGGGGGTTGGAAGTTTTTCCTAGa                                                    >  1:274680/1‑140 (MQ=17)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                                                    >  2:685972/1‑140 (MQ=255)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                                                    >  2:114940/1‑140 (MQ=255)                                           gCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGa                                                    >  2:723208/1‑140 (MQ=255)                                             cACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTATGGGGTTAGAAGTTTTTCCTAGAcc                                                  >  1:242816/1‑140 (MQ=255)                                              aCCGCTATTATGACCTGTTGCAAGGGCGATATATCTCTCTGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCa                                                 <  2:274680/140‑1 (MQ=21)                                              aCCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCa                                                 >  1:604241/1‑140 (MQ=255)                                                   tattatGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTccc                                            <  1:517398/140‑1 (MQ=255)                                                            ccGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAAt                                   >  1:20100/1‑140 (MQ=255)                                                             cGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAAtt                                  >  1:458587/1‑140 (MQ=255)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGTGTTCGAATTTTTTCCTGGACCATTCCCCTTGCCAATTcc                                >  2:230798/1‑140 (MQ=21)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTCGAATTTTTTCCTTGACCGTTCCCCTTGCCAATTcc                                >  2:345972/1‑140 (MQ=21)                                                               cTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTcc                                >  2:701460/1‑140 (MQ=255)                                                                          atataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGGTAGAAGTTATTCTTAGGCCATTCCCTTTGCAAATTCCATGGCCCaaaa                     >  2:219731/2‑140 (MQ=21)                                                                          atataTCACTCAGGA‑CCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCa                                     >  2:69820/2‑123 (MQ=255)                                                                          atataTCACTCAGGA‑CCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCa                                     <  1:69820/122‑1 (MQ=255)                                                                           tataTCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAg                    <  2:585379/140‑1 (MQ=255)                                                                             tatCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGcc                  <  1:20752/140‑1 (MQ=255)                                                                                  cTCAGGATCCGATTGGAATGAACGGGGGATGTAATTTTTTTCAGTATCAGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGca             <  2:261077/140‑1 (MQ=255)                                                                                    cAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGcaca           >  2:552925/1‑140 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGTGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCAcag          >  1:774601/1‑140 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAg                    >  1:623913/1‑130 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAg                    <  2:623913/130‑1 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCAcag          >  1:470283/1‑140 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCAcag          >  2:453516/1‑140 (MQ=255)                                                                                     aGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCAcag          >  2:785798/1‑140 (MQ=255)                                                                                         tCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCATTGCCAAGTCCATGACCCAAACGCCCTGCACagcagc      >  1:261077/1‑140 (MQ=37)                                                                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCagcaag   <  1:701460/140‑1 (MQ=255)                                                                                            gATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCagcaag   >  1:109346/1‑140 (MQ=255)                                                                                             aTTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCagcaagc  >  1:193420/1‑140 (MQ=255)                                                                                                  |                                                                                                                                       CGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGC  >  NC_000913/3765757‑3765988

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