Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,345 T→C 17.2% pseudogene (954/2518 nt) ydbA → putative outer membrane protein, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,3450TC17.2% 62.2 / 4.9 42pseudogene (954/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (12/18);  new base C (3/4);  total (18/24)
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

CAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATA  >  NC_000913/1466204‑1466487
                                                                                                                                             |                                                                                                                                              
aaGCGCCACGGTGGACAACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATg                                                                                                                                               <  1:321850/142‑1 (MQ=255)
       aCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCACCCAGATTAAGGGtgatga                                                                                                                                        >  2:871872/1‑143 (MQ=255)
       aCGGCGGAGAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGtgatga                                                                                                                                        <  1:705880/143‑1 (MQ=255)
             gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCCTGGGTATCCAGATCGACGGTGACAAGGCCATAGTCAACAACGAGGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAGTAATGGTGATGACGCCAc                                                                                                                                  >  1:93437/1‑143 (MQ=17)
             gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCCCCCAGAGTAATGGTGATGACGGCAc                                                                                                                                  >  1:79264/1‑143 (MQ=255)
             gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATCGACGGCGACAAGGCGATCCTCAACAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGCCACCCAGATTAACGGGGACGAcccccc                                                                                                                                  >  2:91182/1‑141 (MQ=11)
                 aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCGTCGTCAATAACGAAGGCGAGCGCACCATCACCCCCGGCGGCACCGGCACCCAGGTTAATGGTGATGTCGCCACGGCa                                                                                                                              >  1:911429/1‑143 (MQ=18)
                 aCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTAAATAACGAAGGCGAGAGCACCATCACCAAAGATGGCACCGGCACCCAGATTGAGGGTGATGAGGCCACGGCa                                                                                                                              >  1:234567/1‑143 (MQ=18)
                    atggCCCCATGACCGTCACGGGTCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGCTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                           <  1:446891/141‑1 (MQ=18)
                    aGGGCACCATGACCGTCACCGATCCGGAGTCGATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                           <  2:433591/143‑1 (MQ=255)
                              gACCGTCACCGATCCGGGGACCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                                 <  2:541521/143‑1 (MQ=255)
                                   tCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTAAATATTGAAGTCTTGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCACAacc                                                                                                            >  1:278033/1‑143 (MQ=21)
                                           ccGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                    <  2:833296/143‑1 (MQ=255)
                                                  ccATGGGTATCCAGATCGACGGTGACAAGGACATCGTCAATAACGAAGGCGATAGCACCATAAACAACGGTGGAACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                             <  1:949196/143‑1 (MQ=18)
                                                  ccATGGGTACCCACATCGACGGTGACCAGGCCATCGTCACTAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATTGTGATGACGCCACGGCAAACAACAACGGCAGAACCACCGTTGACGGCAAg                                                                                             <  2:385965/143‑1 (MQ=18)
                                                  ccATGGGGATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                             <  2:98431/143‑1 (MQ=255)
                                                   cATGGGTATCCAGATCGACGGTGACAAGGCCATCGACAATAACGAAGGCGAGAGCGCCAGCACGAACGGTGGCACCGGCACCCAGAAGAATGGTGATGAGGCCACGGCCGGGAACAACGGGAAAACCACCGGTGACGGCAAgg                                                                                            >  2:415984/1‑143 (MQ=25)
                                                        gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGc                                                                                                                               >  1:7490/1‑103 (MQ=255)
                                                        gTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGc                                                                                                                               <  2:7490/103‑1 (MQ=255)
                                                         tATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                                                      <  1:888787/143‑1 (MQ=255)
                                                          aTCCAGATCGAGGGTGACAAGGCCAGCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                     <  2:186633/143‑1 (MQ=255)
                                                          aTCCAGATCGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                     <  2:278033/143‑1 (MQ=255)
                                                            ccAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                   >  1:743028/1‑143 (MQ=255)
                                                             caggacGACGGTGACAAGGCCATCGTCAAAAACGAAGGCAAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                                  <  2:229424/138‑1 (MQ=255)
                                                               gATGGACGGTGACAAGGCCATCGTCAAGAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                <  2:186597/143‑1 (MQ=255)
                                                                  cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGCCGGCAAGGATGCCCGCGGTACgg                                                                             >  1:596785/1‑143 (MQ=255)
                                                                    aCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAACGGTTGCACCGGCACCCAGATTAAGGGCGATGACGCCACGCCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAGTCCACCGGTACCGac                                                                           >  2:456104/1‑142 (MQ=17)
                                                                       gTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                              >  1:444576/1‑105 (MQ=255)
                                                                       gTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                              <  2:444576/105‑1 (MQ=255)
                                                                               gCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGACGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGCTTCCACCGGTACGGACATCAACGGtaa                                                                >  2:765420/1‑143 (MQ=255)
                                                                                  aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGtaataa                                                             <  2:968544/143‑1 (MQ=255)
                                                                                   ccGTGAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCGAGAGGCGCTAAGTAATGCTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAAc                                                            <  1:427181/142‑1 (MQ=14)
                                                                                               aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAtt                                                <  1:939833/143‑1 (MQ=255)
                                                                                                  gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                   >  1:3726/1‑105 (MQ=255)
                                                                                                  gCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                                                   <  2:3726/105‑1 (MQ=18)
                                                                                                    gagaGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGACACGGCATACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGCAATAACGGAAAGGTTATTCAGGa                                           >  1:583260/1‑143 (MQ=14)
                                                                                                            cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGCAATCGCCGGCAATAGCGGAAAAGTGATTCAGGACGGCGCTc                                   >  2:929318/1‑143 (MQ=11)
                                                                                                             aTCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                  <  2:891452/143‑1 (MQ=255)
                                                                                                                     cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCa                          <  1:887266/143‑1 (MQ=255)
                                                                                                                       gTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGc                                                                                                                               >  2:455885/1‑40 (MQ=255)
                                                                                                                                     cataTTAATGGTGAGGACGCCACGGCAAACAACAACGGCAAGACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg          <  1:871872/140‑1 (MQ=255)
                                                                                                                                        aTTAAGGGTGATGACGCCACGACAAACACCAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt       <  1:754845/143‑1 (MQ=255)
                                                                                                                                             tGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata  >  1:587410/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
CAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATA  >  NC_000913/1466204‑1466487

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