Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,806,042 (C)6→7 coding (42/759 nt) ydiY ← acid‑inducible putative outer membrane protein YdiY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,0361.C96.4% 92.4 / ‑3.3 28coding (48/759 nt)ydiYacid‑inducible putative outer membrane protein YdiY
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (16/11);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑‑CCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGA  >  NC_000913/1805904‑1806165
                                                                                                                                     |                                                                                                                                    
gCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑‑ccccccc                                                                                                                            >  2:99842/1‑139 (MQ=255)
                                 tGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCaa                                                                                             >  1:567869/1‑140 (MQ=255)
                                   cGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAg                                                                                          >  2:126887/1‑140 (MQ=255)
                                    gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                                                         >  2:146025/1‑140 (MQ=255)
                                      tCAGATTACCTTCGAACTGTTTATTAGCGCCGGAAAGATAATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACCCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCtt                                                                                       <  2:320506/140‑1 (MQ=255)
                                        aGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCAAACATAATGATAGCAGGAACTGTCTTCAAAAGCTTCa                                                                                     >  1:220529/1‑140 (MQ=255)
                                               cTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCa                                                                              <  1:99842/140‑1 (MQ=255)
                                                   gAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAagag                                                                          >  2:421729/1‑140 (MQ=255)
                                                    aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGTAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGt                                                                         >  1:557315/1‑140 (MQ=255)
                                                         ttttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTa                                                                    >  1:592793/1‑140 (MQ=255)
                                                          tttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                                                   <  2:511515/140‑1 (MQ=255)
                                                                   cTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATATCGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACaaaaaa                                                          >  2:410334/1‑140 (MQ=255)
                                                                      gCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAGCTGTCTTCTAAAACTTCATTTATCAGGAGTTGGTACAACAAAAAAaga                                                       >  2:325516/1‑140 (MQ=255)
                                                                                ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcgcg                                             >  1:547811/1‑140 (MQ=255)
                                                                                       aGTAAAAACGGAATCATCGGCAGCTGCATTCAGGGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCgg                                      >  2:46121/1‑140 (MQ=255)
                                                                                        gTAAAAACGGAATCACCGGCAGATGCATTCAGTGACGAAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGa                                     <  1:325516/140‑1 (MQ=255)
                                                                                         tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaa                                    >  2:542394/1‑140 (MQ=255)
                                                                                         tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaa                                    >  2:666490/1‑140 (MQ=255)
                                                                                         tAAAAACGGAATCATCGGAAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaa                                    <  1:46121/140‑1 (MQ=255)
                                                                                          aaaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaaa                                   >  1:419098/1‑140 (MQ=255)
                                                                                          aaaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATGTATCAAGAATCCGTACGACAAAAAAGGAGACCATCGCGGTCCCGGaaa                                   >  2:193712/1‑140 (MQ=255)
                                                                                           aaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAAc                                  >  1:345306/1‑140 (MQ=255)
                                                                                                    tcatcGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAg                         <  2:539212/140‑1 (MQ=255)
                                                                                                        cGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATc                     <  1:421729/140‑1 (MQ=255)
                                                                                                                gCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTa             <  1:598063/140‑1 (MQ=255)
                                                                                                                 cATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAGAGACCATCGCGGTCCCGGAAACCTTCTTAAGGATCAAAGATTAGc           <  1:126887/140‑1 (MQ=255)
                                                                                                                    tCAGTGACGCAAACACG‑‑CCCCCTCGCCAACATAACGAGGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGt         <  1:193712/140‑1 (MQ=255)
                                                                                                                     cAGTGACGCAAACATG‑‑CACCCCGGCCAACATAACGGTGGCAGGAACTGACTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTc        <  2:287634/140‑1 (MQ=255)
                                                                                                                          aCGCAAACATGCCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTg    <  2:567869/140‑1 (MQ=255)
                                                                                                                           cGCAAACATGCCCCCCCGCCAATATCACGTTGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGa  <  1:271010/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATG‑‑CCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGA  >  NC_000913/1805904‑1806165

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