Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,472,611 G→A 15.5% G144G (GGC→GGT btuB ← vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,472,6110GA15.5% 113.1 / 12.5 45G144G (GGC→GGTbtuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (3/4);  total (20/25)
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.77e-01

TTGCTTCCCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAA  >  AP009048/3472543‑3472681
                                                                    |                                                                      
ttGCTTCCCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGcc                                                                      <  1:3802195/71‑1 (MQ=255)
              tgctgAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAAc                                                        >  1:2754331/1‑71 (MQ=255)
                  gAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCAt                                                     >  1:3576507/1‑70 (MQ=255)
                  gAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCAt                                                     >  1:175722/1‑70 (MQ=255)
                   aaaTTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATa                                                    <  1:3161077/70‑1 (MQ=255)
                      tttCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACa                                                <  1:1877741/71‑1 (MQ=255)
                      tttCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACa                                                <  1:740064/71‑1 (MQ=255)
                         ccgttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAg                                               <  1:4267718/69‑1 (MQ=255)
                          cgttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCAt                                                              >  1:1868634/1‑53 (MQ=255)
                          cgttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCg                                             <  1:3682337/70‑1 (MQ=255)
                            ttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCg                                             <  1:333277/68‑1 (MQ=255)
                            ttccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAg                                          <  1:1928278/71‑1 (MQ=255)
                             tccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTa                                                                    <  1:2036961/44‑1 (MQ=255)
                             tccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTAt                                                                   >  1:1165898/1‑45 (MQ=255)
                              ccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATa                                                    <  1:2922426/59‑1 (MQ=255)
                              ccgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAgcg                                        >  1:4472316/1‑71 (MQ=255)
                               cgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgcg                                        <  1:4384579/70‑1 (MQ=255)
                               cgGGTTCATCGCGCGTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAgcgc                                       <  1:2813526/71‑1 (MQ=255)
                                gggTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATaa                                                   >  1:2049283/1‑58 (MQ=255)
                                gggCTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGcgcg                                      <  1:3695200/71‑1 (MQ=255)
                                   ttCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGcgcg                                      <  1:2332830/68‑1 (MQ=255)
                                    tCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGcc                                   >  1:155958/1‑70 (MQ=255)
                                          cgcgTCGTGATGATATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGat                            <  1:2438413/71‑1 (MQ=255)
                                           gcgTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGata                           >  1:3317586/1‑71 (MQ=255)
                                            cgtcgtGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgc                                         >  1:387711/1‑56 (MQ=255)
                                             gtcgtGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGatat                          <  1:2261983/70‑1 (MQ=255)
                                               cgtGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAg                                          <  1:972212/52‑1 (MQ=255)
                                               cgtGATGATATTCAACACCCCGCCTATTGCATCGGAACCATAAACAGCg                                             >  1:911255/1‑49 (MQ=255)
                                                  gatgatATTCACCACCCCGCCTATTGCATCGGAACCATaaa                                                  >  1:347252/1‑41 (MQ=255)
                                                   atgatATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAg                                          >  1:1750751/1‑48 (MQ=255)
                                                   atgatATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCaa                    >  1:4346298/1‑70 (MQ=255)
                                                   atgatATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCaa                    >  1:3244239/1‑70 (MQ=255)
                                                    tgatATTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAaca                  >  1:3475371/1‑71 (MQ=255)
                                                      atatTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAgcg                                        <  1:2603604/47‑1 (MQ=255)
                                                      atatTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAacac                 <  1:3764373/70‑1 (MQ=255)
                                                       tatTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGc               <  1:195642/71‑1 (MQ=255)
                                                       tatTCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGc               >  1:3091842/1‑71 (MQ=255)
                                                        atTCACCACCCCGCCTATTGCATCGGAACCATAAACAGc                                              >  1:4550399/1‑39 (MQ=255)
                                                         ttCACCACCCCACCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTg             <  1:2869519/71‑1 (MQ=255)
                                                            accaccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGAc          <  1:2095846/71‑1 (MQ=255)
                                                            accaccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGAc          <  1:609372/71‑1 (MQ=255)
                                                              caccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACaa        <  1:2915821/71‑1 (MQ=255)
                                                              caccCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACaa        <  1:4239146/71‑1 (MQ=255)
                                                                  ccGCATATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAgcg     >  1:2931004/1‑70 (MQ=255)
                                                                   cGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGa           >  1:4139522/1‑63 (MQ=255)
                                                                    gCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCaa  >  1:4097727/1‑71 (MQ=255)
                                                                    |                                                                      
TTGCTTCCCCACCCTGCTGAAATTTCCGTTCCGGGTTCATCGCGCGTCGTGATGATATTCACCACCCCGCCTATTGCATCGGAACCATAAACAGCGGAGCGCGGCCCACGGATATATTCAACACGCTGGACAAGCGCAA  >  AP009048/3472543‑3472681

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