Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,190,398 T→G M173R (ATG→AGG)  yqiH → putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,190,3980TG98.1% 177.8 / NA 53M173R (ATG→AGG) yqiHputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (21/31);  minor base A (0/1);  total (21/32)
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

CCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAA  >  NC_000913/3190301‑3190494
                                                                                                 |                                                                                                
ccAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGc                                                                                                <  1:646745/100‑1 (MQ=255)
ccAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGc                                                                                                >  1:958695/1‑100 (MQ=255)
    aaGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTa                                                                                            <  2:51932/100‑1 (MQ=255)
      gCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATc                                                                                          <  1:628980/100‑1 (MQ=255)
      gCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATc                                                                                          <  1:511498/100‑1 (MQ=255)
         aTTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACg                                                                                       >  1:223315/1‑100 (MQ=255)
           tGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGtt                                                                                     <  2:140939/100‑1 (MQ=255)
            gAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTa                                                                                    >  1:359091/1‑100 (MQ=255)
                tAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCt                                                                                >  2:335221/1‑100 (MQ=255)
                tAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCt                                                                                >  2:112183/1‑100 (MQ=255)
                tAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCt                                                                                >  2:41267/1‑100 (MQ=255)
                tAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCt                                                                                <  1:824368/100‑1 (MQ=255)
                 agaagaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCt                                                                                <  2:134662/99‑1 (MQ=255)
                 agaagaCGATGGGGTAATGACTTACGAAAAAGTAGAAATCATTCGTAGAAATGAAGGTTAGATCCGCTTTAACAATAAGAGGCCCTATCACGTTAAGc                                                                                 >  2:443215/1‑98 (MQ=255)
                   aagaCGATTGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTggg                                                                             <  2:621087/100‑1 (MQ=255)
                   aagaCGATGTGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTggg                                                                             <  2:1273558/100‑1 (MQ=255)
                    agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGt                                                                            <  2:359088/100‑1 (MQ=255)
                    agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGt                                                                            <  2:25712/100‑1 (MQ=255)
                    agaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGt                                                                            >  1:1086104/1‑100 (MQ=255)
                     gaCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGtt                                                                           <  1:921974/100‑1 (MQ=255)
                      aCGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTta                                                                          <  2:480931/100‑1 (MQ=255)
                       cGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTta                                                                          >  1:1020719/1‑99 (MQ=255)
                            gggTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCgg                                                                    >  2:49400/1‑100 (MQ=255)
                             aatAACGACTTACAAACAAGTTGAAATCACTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAAGCCTTATCACCTTACGCTCGGTTATATCgg                                                                    <  2:286460/97‑1 (MQ=255)
                                    aCTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCAATCCGCTTTAACAATAAGAGGCCTTATCACGGTACGCTGGATTATATCGGAACaaa                                                              >  2:1024348/1‑98 (MQ=255)
                                      ttACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATggg                                                          <  2:141525/100‑1 (MQ=255)
                                      ttACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATggg                                                          <  1:335224/100‑1 (MQ=255)
                                        aCGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGt                                                        >  1:510449/1‑100 (MQ=255)
                                            aaaaGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACg                                                    <  1:385735/100‑1 (MQ=255)
                                               aGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATg                                                 <  2:928381/100‑1 (MQ=255)
                                               aGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATg                                                 <  1:519019/100‑1 (MQ=255)
                                                 ttGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGtt                                               <  1:1144550/100‑1 (MQ=255)
                                                 ttGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGtt                                               <  2:294553/100‑1 (MQ=255)
                                                 ttGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGtt                                               <  1:1245690/100‑1 (MQ=255)
                                                  tGAAATCATTCGTAGAAATTACGTTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTa                                              >  2:339810/1‑100 (MQ=255)
                                                   gAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAc                                             >  2:864346/1‑100 (MQ=255)
                                                   gAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAc                                             <  1:139000/100‑1 (MQ=255)
                                                    aaaTCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAcc                                            <  1:418009/100‑1 (MQ=255)
                                                    aaaTCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTAcc                                            <  1:737131/100‑1 (MQ=255)
                                                           tCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACg                                     <  2:1155440/100‑1 (MQ=255)
                                                            cGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGc                                    >  1:966442/1‑100 (MQ=255)
                                                            cGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGc                                    >  1:1197571/1‑100 (MQ=255)
                                                                gAAATGACGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAg                                >  1:720354/1‑100 (MQ=255)
                                                                      aCGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAAt                          >  1:152440/1‑100 (MQ=255)
                                                                       cGGTTCGATCCGCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATg                         >  1:624978/1‑100 (MQ=255)
                                                                        ggTTCGATCCGCTTTAACNATAAGAGGCCTTNTCANGTTACGCNGGGTNANATNGGAANAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATgg                        <  1:868311/100‑1 (MQ=255)
                                                                          ttggATCCTCTTTAAAAATAAGAGGCCTTGTCACGCTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTCCCACAGACGCAAAGCTTAATgg                        <  2:402867/95‑1 (MQ=255)
                                                                                  gCTTTAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCttt              >  1:57888/1‑100 (MQ=255)
                                                                                    tttAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTa            <  1:21336/100‑1 (MQ=255)
                                                                                     ttAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAg           >  2:664425/1‑100 (MQ=255)
                                                                                      tAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGc          <  2:333981/100‑1 (MQ=255)
                                                                                      tAACAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGc          <  1:864352/100‑1 (MQ=255)
                                                                                       aaCAATAAGAGGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCt         <  1:810039/100‑1 (MQ=255)
                                                                                              agagGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCaaa  >  1:307708/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CCAAAAGCCATTGAATTAGAAGACGATGGGGTAATGACTTACGAAAAAGTTGAAATCATTCGTAGAAATGACGGTTCGATCCGCTTTAACAATAAGATGCCTTATCACGTTACGCTGGGTTATATCGGAACAAATGGGGTAACGATGTTACCACAAACGCAAAGCTTAATGGTAACCCCTTTTAGCTATGCAAA  >  NC_000913/3190301‑3190494

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: