Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,012,700 C→T *105* (TAG→TAA fliE ← flagellar basal‑body component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,012,7000CT97.8% 146.2 / ‑3.7 46coding (315/315 nt)fliEflagellar basal‑body component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/22);  minor base . (1/0);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCA  >  NC_000913/2012572‑2012803
                                                                                                                                |                                                                                                       
tataCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCa                                                                                               <  1:293812/139‑1 (MQ=255)
tataCGGATCGTGGATGGCAGAAAACGTTCAGGATCAAGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCa                                                                                               >  2:1185171/1‑139 (MQ=255)
   aCGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGc                                                                                            <  2:509500/139‑1 (MQ=255)
                                gATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGc                                                               >  2:551701/1‑139 (MQ=255)
                                         tCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGc                                                      <  1:551701/139‑1 (MQ=255)
                                           cATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCAc                                                    >  1:1447749/1‑139 (MQ=255)
                                              tttAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTg                                                 >  1:740184/1‑139 (MQ=255)
                                               ttAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGa                                                <  2:477929/139‑1 (MQ=255)
                                                       gAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                        <  2:423215/91‑1 (MQ=255)
                                                       gAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCtcat                                                                                        >  1:423215/1‑91 (MQ=255)
                                                       gAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCAt                                        <  2:740184/139‑1 (MQ=255)
                                                        aaaTTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAAC‑‑‑‑T‑‑‑‑‑‑‑‑TACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACt                           >  2:1649013/1‑139 (MQ=255)
                                                               gtgagtTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAg                                <  2:1447749/139‑1 (MQ=255)
                                                                      tGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGa                         >  2:661744/1‑139 (MQ=255)
                                                                         cccATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGc                      <  1:1185171/139‑1 (MQ=255)
                                                                                 tGCCCCACGCGGTGAGATCGGTTGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                    <  2:1101409/85‑1 (MQ=255)
                                                                                 tGCCCCACACGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGcca                                                                    >  1:1101409/1‑85 (MQ=255)
                                                                                  gCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCat             >  2:251775/1‑139 (MQ=255)
                                                                                       aCGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatt                                                          <  1:60936/89‑1 (MQ=255)
                                                                                       aCGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatt                                                          >  2:60936/1‑89 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                             <  1:1021374/53‑1 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATg                                                                                             >  2:1021374/1‑53 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCAcc                                                   <  2:1311489/95‑1 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCAcc                                                   >  1:1311489/1‑95 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                             >  2:1371282/1‑101 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                             <  1:1371282/101‑1 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCttt                   >  1:1444917/1‑127 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCttt                   <  2:1444917/127‑1 (MQ=255)
                                                                                        cgcgGTGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGt       >  1:1450754/1‑139 (MQ=255)
                                                                                             tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACt                                                                                    >  1:213208/1‑57 (MQ=255)
                                                                                             tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACt                                                                                    <  2:213208/57‑1 (MQ=255)
                                                                                             tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacc                                                                  <  2:634949/75‑1 (MQ=255)
                                                                                             tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGccacc                                                                  >  1:634949/1‑75 (MQ=255)
                                                                                             tGAGATCGGATGTAATTAACAACAGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtcatca  >  2:1213758/1‑139 (MQ=255)
                                                                                                               acaacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                             >  1:987038/1‑78 (MQ=255)
                                                                                                               acaacaGGCTTCAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATc                                             <  2:987038/78‑1 (MQ=255)
                                                                                                                caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcgcc                                                                        <  2:300249/50‑1 (MQ=255)
                                                                                                                caacaGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATAcgcc                                                                        >  1:300249/1‑50 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                         <  2:1190672/62‑1 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCttatta                                                         >  1:1190672/1‑62 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTg                          >  2:1344664/1‑93 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTg                          <  1:1344664/93‑1 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtca     <  2:1634642/114‑1 (MQ=255)
                                                                                                                   caGGTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGtca     >  1:1634642/1‑114 (MQ=255)
                                                                                                                      gTTACAAAACTTACACCTGCATGCTCATCACTTCTTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCata            >  1:1169204/1‑104 (MQ=255)
                                                                                                                      gTTACAAAACTTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCata            <  2:1169204/104‑1 (MQ=255)
                                                                                                                                |                                                                                                       
TATACGGATCGTGGATGGCAGAAAACGTTCAGGATCAGGTATCCATTTTAAACCAGAAATTGAGTGAGTTTGCCCCATCCATGCCCCACGCGGTGAGATCGGATGTAATTAACAACAGGTTACAAAACCTACACCTGCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCA  >  NC_000913/2012572‑2012803

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