Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,964,688	G→A	G121G (GGC→GGT)	flhA ←	flagellar biosynthesis protein FlhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,964,688	0	G	A	100.0%	140.0 / NA	42	G121G (GGC→GGT)	flhA	flagellar biosynthesis protein FlhA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/19);  total (23/19)

CGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAA  >  NC_000913/1964553‑1964821                                                                                                                                        |                                                                                                                                      cGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTAccac                                                                                                                                   <  2:80326/140‑1 (MQ=255)         gCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAgg                                                                                                                           <  1:167999/140‑1 (MQ=255)           ttACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGaa                                                                                                                         <  1:537741/140‑1 (MQ=255)                  gCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCg                                                                                                                  >  1:365576/1‑140 (MQ=255)                    aTACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGaa                                                                                                                >  2:44675/1‑140 (MQ=255)                     tACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAAc                                                                                                               <  1:706254/140‑1 (MQ=255)                             agaACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGac                                                                                                       >  1:99020/1‑140 (MQ=255)                                 cAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGaccacc                                                                                                   <  2:505873/140‑1 (MQ=255)                                    aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTc                                                                                                <  2:275355/140‑1 (MQ=255)                                      cgcgcACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTccc                                                                                              >  2:192655/1‑140 (MQ=255)                                              cgcTTCTGCGATACGCCGGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAGCGCCACGATGCCGATAGCGAAATTGCCACCAGCGAGGAAGTGGCCGAACGCTTCGACCGCCTTCCCTggcgccg                                                                                      <  2:137111/138‑1 (MQ=255)                                                  tCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCcgccgccgc                                                                                  <  1:188116/140‑1 (MQ=255)                                                      cGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcgccg                                                                              >  1:626227/1‑140 (MQ=255)                                                       gATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCgg                                                                             <  1:192655/140‑1 (MQ=255)                                                           cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATg                                                                         >  1:361226/1‑140 (MQ=255)                                                              ccGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                      >  2:856727/1‑140 (MQ=255)                                                               cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCt                                                                     <  2:365576/140‑1 (MQ=255)                                                                               gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                     >  1:314248/1‑140 (MQ=255)                                                                               gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATAc                                                     >  1:380555/1‑140 (MQ=255)                                                                                    tAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGtt                                                >  1:602921/1‑140 (MQ=255)                                                                                        gTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCCACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAg                                            <  1:735647/140‑1 (MQ=255)                                                                                          tgatgaTCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGcc                                          >  1:16004/1‑140 (MQ=255)                                                                                            atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCAc                                        >  2:835949/1‑140 (MQ=255)                                                                                               atCACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGtt                                     <  1:835949/140‑1 (MQ=255)                                                                                                 cACGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTaa                                   >  2:551200/1‑140 (MQ=255)                                                                                                   cGAGAATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGt                                 >  1:805869/1‑140 (MQ=255)                                                                                                      gaATGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                              <  1:551200/140‑1 (MQ=255)                                                                                                         tGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGa                           >  2:104829/1‑140 (MQ=255)                                                                                                         tGACAAACACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGa                           >  2:201529/1‑140 (MQ=255)                                                                                                              aaCACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCaa                      >  2:39984/1‑140 (MQ=255)                                                                                                              aaCACCACGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTCCCAGACGCaa                      >  2:629008/1‑140 (MQ=255)                                                                                                                caccacGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACa                    <  2:512636/140‑1 (MQ=255)                                                                                                                   cacGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCg                 <  2:805869/140‑1 (MQ=255)                                                                                                                    acGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGt                >  2:596679/1‑140 (MQ=255)                                                                                                                     cGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCCCCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCTCGTTAAGTGCCAGACGCAACAGCGTg               <  2:314248/140‑1 (MQ=255)                                                                                                                     cGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTg               >  2:12173/1‑140 (MQ=255)                                                                                                                     cGATACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTg               <  2:380555/140‑1 (MQ=255)                                                                                                                       aTACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGt             >  2:498811/1‑140 (MQ=255)                                                                                                                       aTACCGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGt             <  1:201529/140‑1 (MQ=255)                                                                                                                           cGATAGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAaac         <  1:44675/140‑1 (MQ=255)                                                                                                                               aGCGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAacaaca     <  2:320989/140‑1 (MQ=255)                                                                                                                                 cGAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGa   >  2:860998/1‑140 (MQ=255)                                                                                                                                  gAAATTACCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGaa  >  2:631230/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      CGCCATCTGCTTACCCGGCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAACAGAA  >  NC_000913/1964553‑1964821

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