Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,368,320 G→A 7.1% Q94Q (CAG→CAA)  arnA → fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,3200GA 7.1% 177.4 / 3.4 56Q94Q (CAG→CAA) arnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATGGTGAAACGTGC  >  NC_000913/2368176‑2368454
                                                                                                                                                |                                                                                                                                      
tCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAgctcg                                                                                                                                    >  1:505892/1‑149 (MQ=255)
    cTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTcc                                                                                                                                >  2:230398/1‑149 (MQ=255)
           cggAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGt                                                                                                                         <  1:254069/149‑1 (MQ=255)
            ggAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTg                                                                                                                        >  4:861/1‑149 (MQ=255)
             gAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGc                                                                                                                       >  1:267810/1‑149 (MQ=255)
                   ggCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTaa                                                                                                                 <  3:232734/149‑1 (MQ=255)
                    gCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAAt                                                                                                                >  3:186775/1‑149 (MQ=255)
                      aTTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCt                                                                                                              >  4:489949/1‑149 (MQ=255)
                      aTTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCt                                                                                                              <  4:149408/149‑1 (MQ=255)
                       ttCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTg                                                                                                             >  2:357411/1‑149 (MQ=255)
                       ttCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTg                                                                                                             >  2:308456/1‑149 (MQ=255)
                          cGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCAt                                                                                                          >  1:255828/1‑149 (MQ=255)
                                  gcgcCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCt                                                                                                  <  4:419043/149‑1 (MQ=255)
                                       ggATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTAc                                                                                             <  3:347994/149‑1 (MQ=255)
                                        gATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTAcc                                                                                            >  2:350120/1‑149 (MQ=255)
                                           aaCGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCt                                                                                                              >  2:565278/1‑128 (MQ=255)
                                           aaCGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCt                                                                                                              <  1:565278/128‑1 (MQ=255)
                                           aaCGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCaaa                                                                                         <  4:337361/149‑1 (MQ=255)
                                            aCGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCaaaa                                                                                        <  2:572539/149‑1 (MQ=255)
                                                tAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATc                                                                                    <  1:230398/149‑1 (MQ=255)
                                                  atcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGt                                                                                  <  1:219504/149‑1 (MQ=255)
                                                  atcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGt                                                                                  <  1:494289/149‑1 (MQ=255)
                                                  atcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAACTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGt                                                                                  >  2:193940/1‑149 (MQ=255)
                                                  atcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAACTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGt                                                                                  >  4:166144/1‑149 (MQ=255)
                                                   tcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTg                                                                                 >  3:248240/1‑149 (MQ=255)
                                                   tcatcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTg                                                                                 <  3:538804/149‑1 (MQ=255)
                                                      tcCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTc                                                                              >  3:157939/1‑149 (MQ=255)
                                                                    aaCGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGTTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAAc                                                                <  1:501272/149‑1 (MQ=255)
                                                                        cATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTggg                                                            >  1:206312/1‑149 (MQ=255)
                                                                        cATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTggg                                                            <  1:25729/149‑1 (MQ=255)
                                                                                aaCTGTCGCCAGATGTGATTTTCTCTTTTTATTCTCGCCATCTTGTTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTc                                                    >  2:569184/1‑149 (MQ=255)
                                                                                      cGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGt                                              <  2:337316/149‑1 (MQ=255)
                                                                                      cGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGt                                              <  2:464688/149‑1 (MQ=255)
                                                                                        ccAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGa                                            >  1:67321/1‑149 (MQ=255)
                                                                                                 aTTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAAc                                   <  3:489949/149‑1 (MQ=255)
                                                                                                 aTTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAAc                                   <  3:164397/149‑1 (MQ=255)
                                                                                                    ttCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTgg                                >  2:361828/1‑149 (MQ=255)
                                                                                                    ttCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTgg                                >  4:342848/1‑149 (MQ=255)
                                                                                                     tctcTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGc                               <  4:106142/149‑1 (MQ=255)
                                                                                                     tctcTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGc                               <  1:569184/149‑1 (MQ=255)
                                                                                                         tttttATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTa                           >  1:177912/1‑149 (MQ=255)
                                                                                                         tttttATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTa                           >  2:534500/1‑149 (MQ=255)
                                                                                                          ttttATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTAc                          <  1:529078/149‑1 (MQ=255)
                                                                                                                tatCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCACGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCa                    >  1:127389/1‑149 (MQ=255)
                                                                                                                     ccATCTTATTTACGACGAAATTTTGCAACTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGaa               <  3:166144/149‑1 (MQ=255)
                                                                                                                     ccATCTTATTTACGACGAAATTTTGCAACTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGaa               <  1:193940/149‑1 (MQ=255)
                                                                                                                       aTCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCACCGGTGAAACGGAAACTGGCGTTACATTGCACCGAATg             >  4:140146/1‑149 (MQ=255)
                                                                                                                       aTCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATg             >  3:138321/1‑149 (MQ=255)
                                                                                                                        tCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATgg            <  1:225683/149‑1 (MQ=255)
                                                                                                                        tCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATgg            <  2:255828/149‑1 (MQ=255)
                                                                                                                        tCTTATGTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACGTTTCAccattcgg            >  2:283765/1‑143 (MQ=255)
                                                                                                                         cTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACTGTGAAACGGAAACTGGCGTTACATTGCACCGAATGGt           <  1:350120/149‑1 (MQ=255)
                                                                                                                         cTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATGGt           <  1:465257/149‑1 (MQ=255)
                                                                                                                                  gacgaAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCGTGGTTCGCTGTTACCAAATTATCGTGGTCGCGAGTCGCTGAACTGGGTGCTGGTCAACGGTGACACGGAAACTAGCGTTACATTGCACCGTATGGTGAAACGTGc  >  3:473888/1‑149 (MQ=255)
                                                                                                                                  gacgaAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAAc                                   >  2:70145/1‑116 (MQ=255)
                                                                                                                                  gacgaAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAAc                                   <  1:70145/116‑1 (MQ=255)
                                                                                                                                  gacgaAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATGGTGAAACGTGc  <  4:304863/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
TCGTCTGGCGGCGGAAAGAGGCATTCCGGTTTATGCGCCGGATAACGTTAATCATCCGCTGTGGGTGGAACGCATTGCCCAACTGTCGCCAGATGTGATTTTCTCTTTTTATTATCGCCATCTTATTTACGACGAAATTTTGCAGCTCGCTCCCGCAGGTGCATTTAATCTGCATGGTTCGCTGTTACCAAAATATCGTGGTCGCGCGCCGCTGAACTGGGTGCTGGTCAACGGTGAAACGGAAACTGGCGTTACATTGCACCGAATGGTGAAACGTGC  >  NC_000913/2368176‑2368454

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