Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,944 A→G V72V (GTA→GTG)  flgG → flagellar basal‑body rod protein FlgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,9440AG96.9% 101.7 / NA 32V72V (GTA→GTG) flgGflagellar basal‑body rod protein FlgG
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (14/17);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTACTAACGGTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTC  >  NC_000913/1134813‑1135060
                                                                                                                                   |                                                                                                                    
gTCAGTACTAACGGTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGt                                                                                                              <  2:693906/140‑1 (MQ=255)
          aCGGTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGaa                                                                                                    >  2:678059/1‑140 (MQ=255)
          aCGGTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGaa                                                                                                    >  1:343033/1‑140 (MQ=255)
            ggTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACg                                                                                                  <  2:58054/140‑1 (MQ=255)
                  aaGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTaca                                                                                            <  2:505328/140‑1 (MQ=255)
                  aaGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTaca                                                                                            <  1:52556/140‑1 (MQ=255)
                  aaGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTaca                                                                                            <  1:505826/140‑1 (MQ=255)
                           cgcgcgGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAggg                                                                                   <  2:716392/140‑1 (MQ=255)
                                    tttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTaca                                                                                            <  1:457651/122‑1 (MQ=255)
                                    tttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTaca                                                                                            >  2:457651/1‑122 (MQ=255)
                                     ttGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCg                                                                         <  1:105071/140‑1 (MQ=255)
                                         aGATCTGCTTTCTCAATCCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCGCCTTACCCTCCGGATTACAAATCGGCACGGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                                                     <  2:608217/140‑1 (MQ=255)
                                            tCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGc                                                                        <  1:547591/134‑1 (MQ=255)
                                            tCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGc                                                                        >  2:547591/1‑134 (MQ=255)
                                            tCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                                                     >  2:91510/1‑137 (MQ=255)
                                            tCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                                                     <  1:91510/137‑1 (MQ=255)
                                            tCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACca                                                                  >  1:262732/1‑140 (MQ=255)
                                              tGCTTTATCAAACCATTCGCCAGCCGGGGGCATAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                                                     <  1:268703/135‑1 (MQ=255)
                                              tGCTTTATCAAACCATTCGCCAGCCGGGGGCATAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGa                                                                     >  2:268703/1‑135 (MQ=255)
                                              tGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaac                                                                >  1:386923/1‑140 (MQ=255)
                                              tGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaac                                                                >  2:668187/1‑140 (MQ=255)
                                                cTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACcaacaa                                                              >  1:400115/1‑140 (MQ=255)
                                                 tttATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCaacaac                                                             <  2:343033/140‑1 (MQ=255)
                                                          ccATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCa                                                         <  1:119174/135‑1 (MQ=255)
                                                          ccATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCa                                                         >  2:119174/1‑135 (MQ=255)
                                                               cGCCAGCCGGGGGCACAGTCTTCCGAACAAACCGCCTTTCCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTcgc                                               <  1:536681/140‑1 (MQ=255)
                                                                    gCCGGGGGCACAGTCGTCCGAACAAACCCCCTTGCCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGGCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTa                                          <  2:366943/140‑1 (MQ=255)
                                                                                gTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCt                              <  1:96571/140‑1 (MQ=255)
                                                                                   ttCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCtttt                           >  2:608340/1‑140 (MQ=255)
                                                                                    tCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCttttt                          >  1:653349/1‑140 (MQ=255)
                                                                                      cGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTcc                        <  2:400115/140‑1 (MQ=255)
                                                                                              ccaccTTACCCTCCGGATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATg                <  1:668187/140‑1 (MQ=255)
                                                                                                            ggATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTtc  >  2:259611/1‑140 (MQ=255)
                                                                                                            ggATTACAAATCGGCACGGGGGTGCGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTtc  >  2:726722/1‑140 (MQ=255)
                                                                                                                                   |                                                                                                                    
GTCAGTACTAACGGTTTTAAGCGTCAGCGCGCGGTGTTTGAAGATCTGCTTTATCAAACCATTCGCCAGCCGGGGGCACAGTCTTCCGAACAAACCACCTTACCCTCCGGATTACAAATCGGCACGGGGGTACGCCCGGTCGCCACTGAACGCTTACACAGCCAGGGAAACCTGTCGCAGACCAACAACAGCAAAGATGTCGCGATTAAAGGGCAGGGCTTTTTCCAGGTGATGTTGCCAGATGGTTC  >  NC_000913/1134813‑1135060

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