Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,915,227 G→T 10.7% intergenic (‑27/+326) glmU ← / ← atpC fused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase/ATP synthase F1 complex subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,915,2270GT10.7% 115.1 / 4.3 46intergenic (‑27/+326)glmU/atpCfused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase/ATP synthase F1 complex subunit epsilon
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base T (2/3);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAA  >  NC_000913/3915089‑3915366
                                                                                                                                          |                                                                                                                                           
gCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCCGACTGCAATTTGCGCACGAAtt                                                                                                                                         >  2:801409/1‑143 (MQ=255)
     aCCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAc                                                                                                                                    <  1:94318/143‑1 (MQ=255)
       cATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCg                                                                                                                                  >  2:657092/1‑143 (MQ=255)
        aTCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGc                                                                                                                                 >  2:463482/1‑143 (MQ=255)
             gtccccGGCAAGTGTATGCAGCGCTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCAc                                                                                                                            <  2:372135/140‑1 (MQ=255)
              ttCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGAGTGCCTTTGCCTGCGGCAAGGATAAGTAAGCTAATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCAcc                                                                                                                           >  2:133127/1‑143 (MQ=255)
                 ccGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTGATTTGAGAACGAATTTAAACGGCTTCACCtgg                                                                                                                        >  1:53027/1‑141 (MQ=255)
                 ccGGCAAGGGGATGCAGCACGTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGCATGAATTTAAACCGCTTCACCTTg                                                                                                                        <  2:716265/143‑1 (MQ=255)
                                 gcaCTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGCAAAAACTACAttttt                                                                                                        >  1:714457/1‑143 (MQ=255)
                                     tttCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACTAATTTAAACCGCTTCACCTTGAAGAAACTACATTTTTTTCa                                                                                                    >  1:671557/1‑143 (MQ=255)
                                     tttCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACTAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCa                                                                                                    >  1:142828/1‑143 (MQ=255)
                                          gcaGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTg                                                                                               <  2:258422/141‑1 (MQ=255)
                                          gAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTg                                                                                               >  1:415492/1‑143 (MQ=255)
                                           aaGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATAACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTGTTTTCATCGTGa                                                                                              >  2:406216/1‑143 (MQ=255)
                                            aGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGaa                                                                                             <  2:293792/143‑1 (MQ=255)
                                            aGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGaa                                                                                             <  1:523293/143‑1 (MQ=255)
                                              aTCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAAt                                                                                           <  1:21551/143‑1 (MQ=255)
                                                cGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATgg                                                                                         >  2:503452/1‑143 (MQ=255)
                                                        aTGCGCGTGCCGTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGa                                                                                 <  1:283338/143‑1 (MQ=255)
                                                          gcgcGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATa                                                                               <  1:657092/143‑1 (MQ=255)
                                                               tGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAACTGGACAGAGGATAAATTg                                                                          >  1:609397/1‑143 (MQ=255)
                                                                  cTTGGCCTGCGGCAAGGATCACTTCGCTCGTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc                                                                       <  2:199893/143‑1 (MQ=255)
                                                                     tGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt                                                                    >  2:317261/1‑143 (MQ=255)
                                                                           cGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                                                              <  2:197609/143‑1 (MQ=255)
                                                                           cGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa                                                              <  2:610457/143‑1 (MQ=255)
                                                                            ggCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTTAAATGGACAGAGGCTAAATTGTGCAATAACGGAt                                                             >  2:240648/1‑143 (MQ=255)
                                                                             gCAAGGGTCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACACTTTTTTCATCGTGAAATGCCCACAGGATAAATTGCTCAATAACGGAtc                                                            >  2:793758/1‑142 (MQ=255)
                                                                                        tACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa                                                 <  1:89803/143‑1 (MQ=255)
                                                                                            cTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa                                             <  2:139377/143‑1 (MQ=255)
                                                                                            cTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa                                             >  1:542265/1‑143 (MQ=255)
                                                                                             ccATAGCATTATTCGACATACGCGTCCTCACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa                                            <  2:269320/142‑1 (MQ=255)
                                                                                               aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc                                          >  2:372812/1‑143 (MQ=255)
                                                                                                        ttCAACATACGCGTCCTGACTGTAATTTGGGAACGAATTTAAACCACTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGACAGCGCCACTTTc                                 >  1:652115/1‑143 (MQ=255)
                                                                                                          cAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGa                               <  1:471790/143‑1 (MQ=255)
                                                                                                           aaCATACGCGTCCTGACTGTAATTTGAGAACGACTTTAAACCGCTTCACCTTGAAAAAACTACATTTATTTCATCGTGAAATCGACAGAGGATAAATTGTTCAATAACGGATTATCACCGCTCACGAAAACGCCAtctgaccc                              >  2:642288/1‑136 (MQ=255)
                                                                                                                 cgcgTCCTGACTGTAATTTGAGAACTAATTTAAACCGCTTCACCTTGAAGAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGt                        <  2:671557/143‑1 (MQ=255)
                                                                                                                 cgcgTCCTGACTGTAATTTGAGAACTAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGt                        <  2:142828/143‑1 (MQ=255)
                                                                                                                  gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAATGGATTATCACTGCTGACGAACGCGCCATTTTCGACCATCGTc                       >  2:7025/1‑143 (MQ=255)
                                                                                                                  gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTc                       >  1:127906/1‑143 (MQ=255)
                                                                                                                  gcgTCCTGACTGTAATTTGAAAACGAATCTACACAGCTTCACCTTGCAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGCTCACTCACGCACTATCCCCGCCAACGACAGCACCATCATccaccaacacc                       >  2:759359/1‑138 (MQ=255)
                                                                                                                   cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCGCCGCTGACGAAAGCGCCATTTTCGACCATAGTCa                      >  2:56016/1‑143 (MQ=255)
                                                                                                                      ccTGACTGTAATTGGAGAACTAATTGAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGc                   <  2:368905/143‑1 (MQ=255)
                                                                                                                         gACTGTAATTTGAGAACGGATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAGTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCtt                <  2:486269/143‑1 (MQ=255)
                                                                                                                                  ttGAGACCGAATTCAAACCGCTTCACCTTGAAAAAACTACACTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCa       <  2:165727/143‑1 (MQ=255)
                                                                                                                                  ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCa       >  1:640994/1‑143 (MQ=255)
                                                                                                                                       aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGaaaa  >  2:564290/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCATTTTCGACCATCGTCACGCCTTTTGTCTGCAGAAAA  >  NC_000913/3915089‑3915366

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