Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,877,041 T→A 7.7% intergenic (‑89/‑311) pucM ← / → puuE 5‑hydroxyisourate hydrolase/allantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,877,0410TA 7.7% 148.1 / 3.2 52intergenic (‑89/‑311)pucM/puuE5‑hydroxyisourate hydrolase/allantoinase
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base A (2/2);  total (26/26)
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.91e-01

ATCGAGCTACCCGGGCAGCCATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCAAACAGGCAGGTTTCTTG  >  NC_002947/4876894‑4877188
                                                                                                                                                    |                                                                                                                                                   
aTCGAGCTACCCGGGCAGCCATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGttatt                                                                                                                                                     <  6:243812/148‑1 (MQ=255)
   gAGCTACCCGGGCAGCCATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAg                                                                                                                                                  >  8:161628/1‑148 (MQ=255)
         cccGGGCAGCCATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATtc                                                                                                                                            <  1:275444/148‑1 (MQ=255)
                   cATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGAt                                                                                                                                 <  8:278407/149‑1 (MQ=255)
                    aTGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGTCTTAATAACGTCATtcccg                                                                                                                                         <  3:238064/140‑3 (MQ=255)
                    aTGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGTCTTAATAACGTCATtcccg                                                                                                                                         >  4:238064/1‑138 (MQ=255)
                         cgGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTaca                                                                                                                           <  4:222990/149‑1 (MQ=255)
                         cgGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTac                                                                                                                            >  2:158708/1‑148 (MQ=255)
                             aTCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCAGATACCCAAAGCCAGTGGCGGCAGGCGGGAATGAAGTTATTATTACATTAAAACCAAGGATTGTagaaaa                                                                                                                        >  1:271350/1‑144 (MQ=255)
                                  gTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAAtttttt                                                                                                                  <  4:119402/149‑1 (MQ=255)
                                      gtgtgtGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTGATTAAGACATTCCATCAAAGGATTGTACACAATTTTTTGATc                                                                                                              >  7:135011/1‑149 (MQ=255)
                                            ggTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGTCTTAATAACGTCATtcccgc                                                                                                                                        <  3:282238/117‑4 (MQ=255)
                                            ggTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGTCTTAATAACGTCATtcccgc                                                                                                                                        >  4:282238/1‑114 (MQ=255)
                                             gTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCt                                                                                                        <  4:187626/148‑1 (MQ=255)
                                             gTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTg                                                                                                       >  4:266429/1‑149 (MQ=255)
                                             gTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTg                                                                                                       >  2:222930/1‑149 (MQ=255)
                                              tCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGt                                                                                                      >  1:14462/1‑149 (MQ=255)
                                              tCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGt                                                                                                      <  3:312676/149‑1 (MQ=255)
                                                aaaCGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTca                                                                                                    <  1:158708/149‑1 (MQ=255)
                                                  aCGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATtcc                                                                                                                                           >  4:51175/1‑108 (MQ=255)
                                                  aCGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATtcc                                                                                                                                           <  3:51175/108‑1 (MQ=255)
                                                  aCGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTcaca                                                                                                  >  2:259737/1‑149 (MQ=255)
                                                  aCGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTcaca                                                                                                  >  4:213006/1‑149 (MQ=255)
                                                      cccATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGa                                                                                                                                                 <  7:45018/98‑1 (MQ=255)
                                                      cccATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGa                                                                                                                                                 >  8:45018/1‑98 (MQ=255)
                                                        cATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGAtt                                                                                                                                >  7:184546/1‑113 (MQ=255)
                                                        cATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGAtt                                                                                                                                <  8:184546/113‑1 (MQ=255)
                                                                   gCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGcc                                                                                  >  1:308243/1‑148 (MQ=255)
                                                                    ccTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCtt                                                                                                          >  3:252736/1‑123 (MQ=255)
                                                                    ccTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCtt                                                                                                          <  4:252736/123‑1 (MQ=255)
                                                                       ctTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATa                                                                              >  5:290162/1‑148 (MQ=255)
                                                                          gTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATaaa                                                                            <  3:266429/147‑1 (MQ=255)
                                                                          gTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCCGCAGGCGGGAATGACGTTATTAAGAAATTCCATCCAAGGATTGTACACAATATTTTGATCCCTACAGTCACACCCTCCATTTGCTTCCCATaaaa                                                                           >  7:194698/1‑147 (MQ=255)
                                                                            tCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGAc                                                                        >  3:20678/1‑149 (MQ=255)
                                                                                tCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCAt                                                                    >  6:55244/1‑149 (MQ=255)
                                                                                    gcggCAAAAAATACTCACCCATCGCACCAGCTTAGCCCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGt                                                                <  8:95288/147‑1 (MQ=255)
                                                                                     ggtcAATAATTACTCATCGATCGCAGCAGGTTAGCGCCT‑GCGGCAGGCTGGAATGAGGTTTTTATGACATTCCATCCTAGGATTGTACTCAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGt                                                                <  8:194698/145‑1 (MQ=255)
                                                                                        cAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGc                                                             >  6:322924/1‑148 (MQ=255)
                                                                                          aaaaaTACTCAGCCATCGCAACATCTTATCGCCT‑GCGGCATGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCaaa                                                          <  7:161628/149‑1 (MQ=255)
                                                                                              aTACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATaca                                                      <  3:213006/149‑1 (MQ=255)
                                                                                                 cTCAGCCATCGCACGAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGc                                                   <  5:320660/149‑1 (MQ=255)
                                                                                                  tCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCa                                                  >  2:309160/1‑149 (MQ=255)
                                                                                                    cgCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCCCACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCAc                                                 >  6:165918/2‑148 (MQ=255)
                                                                                                      ccATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCAt                                              <  5:322924/149‑1 (MQ=255)
                                                                                                        aTCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATaa                                            >  2:321536/1‑149 (MQ=255)
                                                                                                            cACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGAt                                        <  1:222930/149‑1 (MQ=255)
                                                                                                             aCCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATg                                       >  7:185692/1‑149 (MQ=255)
                                                                                                                 gCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCaa                                   <  2:308243/149‑1 (MQ=255)
                                                                                                                  cTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAAc                                  <  3:214595/149‑1 (MQ=255)
                                                                                                                               ggCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGa                      <  6:290162/149‑1 (MQ=255)
                                                                                                                                    gCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCaaa                 >  2:148962/1‑149 (MQ=255)
                                                                                                                                      gggAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCAAAca               <  1:148962/149‑1 (MQ=255)
                                                                                                                                           tGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCAAAcaggcag          <  1:321536/149‑1 (MQ=255)
                                                                                                                                                   ttaAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTCTCACACCCTCCATTTGCTGCCCATAAACGAACCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCAAACAGGCAGGTTTCTTg  <  8:185692/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCGAGCTACCCGGGCAGCCATGCGCGGCATCCAGTACGTGTGTGGTCAAACGTCCCATTGCGTGTCGCCTCTTGTTCAATCTGTGGGCAAAAAATACTCAGCCATCGCACCAGCTTAGCGCCT‑GCGGCAGGCGGGAATGACGTTATTAAGACATTCCATCCAAGGATTGTACACAATTTTTTGATCCCTTCTGTCACACCCTCCATTTGCTGCCCATAAACGACCCATCGGTCGCAAATACACGCACCATAACGATGCCAACCTCATTAGCTGACCAAACAGGCAGGTTTCTTG  >  NC_002947/4876894‑4877188

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