Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,587,691 C→T 8.3% T141T (ACG→ACA)  ugpQ ← glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,6910CT 8.3% 147.4 / 4.2 48T141T (ACG→ACA) ugpQglycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base T (2/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

AGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGG  >  NC_000913/3587545‑3587837
                                                                                                                                                  |                                                                                                                                                  
agaTGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcggcggc                                                                                                                                                    <  2:334518/147‑1 (MQ=255)
agaTGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcggcggc                                                                                                                                                    >  1:334518/1‑147 (MQ=255)
agaTGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGt                                                                                                                                                  >  2:658723/1‑149 (MQ=255)
         gagaCGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGc                                                                                                                                         <  4:257865/149‑1 (MQ=255)
          agaCGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGcc                                                                                                                                        <  1:658723/149‑1 (MQ=255)
           gaCGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGccc                                                                                                                                       <  4:392989/149‑1 (MQ=255)
                  cccagccGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGTGTCATACCGGCCCACAGCTc                                                                                                                                >  4:418639/1‑149 (MQ=255)
                   ccagccGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTcg                                                                                                                               >  2:233936/1‑149 (MQ=255)
                    cagccGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTcgc                                                                                                                              >  3:140859/1‑149 (MQ=255)
                       ccGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTg                                                                                                                           <  1:449537/149‑1 (MQ=255)
                       ccGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTg                                                                                                                           <  2:157833/149‑1 (MQ=255)
                         gcgcgGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgcc                                                                                                                         <  1:209603/149‑1 (MQ=255)
                               tCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAgc                                                                                                                   >  4:155010/1‑149 (MQ=255)
                                   ttCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGcca                                                                                                               <  2:350728/149‑1 (MQ=255)
                                     cgcgCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGTGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGccacc                                                                                                             <  3:418639/149‑1 (MQ=255)
                                       cgcCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAt                                                                                                           <  3:155010/149‑1 (MQ=255)
                                          cAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCTATCTCAAACGATGACAGCAGCGGCGGTGTCATAc                                                                                                                                             <  3:459649/112‑1 (MQ=255)
                                          cAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGTGTCATAc                                                                                                                                             >  4:459649/1‑112 (MQ=255)
                                             tcgtcgCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGcc                                                                                                     >  4:219563/1‑149 (MQ=255)
                                                    gcCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGcca                                                                                                               >  1:174008/1‑132 (MQ=255)
                                                    gcCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGcca                                                                                                               <  2:174008/132‑1 (MQ=255)
                                                          cATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccg                                                                                        <  4:260916/149‑1 (MQ=255)
                                                           aTCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGcc                                                                                                     >  1:160294/1‑135 (MQ=255)
                                                           aTCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGcc                                                                                                     <  2:160294/135‑1 (MQ=255)
                                                           aTCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccgg                                                                                       >  2:425156/1‑149 (MQ=255)
                                                             ccAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtg                                                                                     >  3:438341/1‑149 (MQ=255)
                                                              cAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgc                                                                                    >  1:452035/1‑149 (MQ=255)
                                                               aGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgcc                                                                                   >  4:193868/1‑149 (MQ=255)
                                                                 caacaaACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgccgg                                                                                 >  3:348568/1‑149 (MQ=255)
                                                                     aaaCCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGCTAATGGCCCGGCGCCggtggt                                                                             >  1:62638/1‑149 (MQ=255)
                                                                       aCCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGtggtgg                                                                           <  4:473303/149‑1 (MQ=255)
                                                                                cAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGAttt                                                                  <  1:233936/149‑1 (MQ=255)
                                                                                   ttCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcaca                                                                                                                                    <  4:518756/80‑1 (MQ=255)
                                                                                   ttCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcaca                                                                                                                                    >  3:518756/1‑80 (MQ=255)
                                                                                        gTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTc                                                          >  4:167220/1‑149 (MQ=255)
                                                                                          gccgccTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGc                                                        <  3:219563/149‑1 (MQ=255)
                                                                                             gccTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCcat                                                     >  1:638947/1‑149 (MQ=255)
                                                                                                     tgCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCgt                                             <  2:233900/149‑1 (MQ=255)
                                                                                                            tCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTcgcg                                      <  1:668848/149‑1 (MQ=255)
                                                                                                                aaGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCAACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCAGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCAGCCCGTGTTCGCGGCaa                                  >  4:563263/1‑149 (MQ=255)
                                                                                                                   cATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGt                               >  1:476666/1‑149 (MQ=255)
                                                                                                                           tcAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCAc                       >  2:372077/1‑149 (MQ=255)
                                                                                                                           tcAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCAc                       >  4:596531/1‑149 (MQ=255)
                                                                                                                           tcAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCAc                       >  2:407111/1‑149 (MQ=255)
                                                                                                                            cAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCAcc                      <  4:140859/149‑1 (MQ=255)
                                                                                                                             aaaCGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCt                     >  3:688646/1‑149 (MQ=255)
                                                                                                                             aaaCGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGGGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCt                     >  2:86522/1‑149 (MQ=255)
                                                                                                                               aCGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGc                   >  3:333697/1‑149 (MQ=255)
                                                                                                                                      cagcagCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCa            >  2:76555/1‑149 (MQ=255)
                                                                                                                                        gcagcGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAAc          <  1:425156/149‑1 (MQ=255)
                                                                                                                                            cggcggcgTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTa      >  2:27129/1‑149 (MQ=255)
                                                                                                                                                ggcgTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCgg  >  4:628307/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCGATATTCGCCATCATCCCGTGTTCGCGGCAACGTTCCGCCACCTGCGAAAGCAACGGTAGCGG  >  NC_000913/3587545‑3587837

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: ‑