Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,136,312 C→A 7.2% L20L (CTG→CTT pitB ← metal phosphate:H(+) symporter PitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,136,3120CA 7.2% 142.5 / 4.3 56L20L (CTG→CTTpitBmetal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base A (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

AATAACACGCCAAAAAAGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAAAGTGAAATATCATTTT  >  NC_000913/3136172‑3136454
                                                                                                                                            |                                                                                                                                              
aaTAACACGCCAAAAAAGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCag                                                                                                                                              >  1:648622/1‑143 (MQ=255)
     cacGCCAAAAAAGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGACACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCATagcca                                                                                                                                         >  2:850809/1‑141 (MQ=255)
      acGCCAAAAAAGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGagcaag                                                                                                                                        >  2:261049/1‑143 (MQ=255)
         acAAAAAAGTTAAGAACGGCGGCCATCACCTCAGGAAGTGGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAgcaagcaa                                                                                                                                     <  2:567794/142‑1 (MQ=255)
          gaaaaaaGTTAAAAAAGACCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAAt                                                                                                                                    <  2:992386/142‑1 (MQ=255)
           aaaaaaGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAACAAGCAATa                                                                                                                                   >  2:836078/1‑143 (MQ=255)
              aaaGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATAGATTGCTCCGTAGAACAACACAAATGCCAGCGCAAGCAATAACa                                                                                                                                >  2:699848/1‑143 (MQ=255)
                 gtcaaaaaaAGCTTCCATCTCCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAACAAGCAATAACaaaa                                                                                                                             <  1:134256/140‑1 (MQ=255)
                  ttAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAg                                                                                                                            >  1:751268/1‑143 (MQ=255)
                  ttAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAg                                                                                                                            >  2:180226/1‑143 (MQ=255)
                  ttAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCAAGAGCAAGCAATAACAAAAg                                                                                                                            >  2:592312/1‑143 (MQ=255)
                  ttAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCAAGAGCAAGCAATAACAAAAg                                                                                                                            >  2:99820/1‑143 (MQ=255)
                   tAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTAGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGc                                                                                                                           >  1:70208/1‑143 (MQ=255)
                    aaaaaaTGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGcc                                                                                                                          >  1:50818/1‑143 (MQ=255)
                    aaaaaaTGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAAACACTGATAGCTTCGCAGAACAACACAAATGCCAGCGCAAGCAACAACAAAAGcc                                                                                                                          >  2:401636/1‑143 (MQ=255)
                     aaaaaTGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGccc                                                                                                                         >  1:641848/1‑143 (MQ=255)
                      aaaaTGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCt                                                                                                                        <  1:206215/143‑1 (MQ=255)
                        aaTGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCtgt                                                                                                                      <  1:758545/143‑1 (MQ=255)
                           gctgccATCACGAGAGCAGGTTGTGCGGGCATGGCATGATTATAAAGAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCGTTGAATGATTTGAAGAATAAGACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGta                                                                                                                   <  2:454101/140‑1 (MQ=255)
                           gccgccATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAGCGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGta                                                                                                                   <  2:47505/143‑1 (MQ=255)
                             cgccATCACCACAGCAAGTGGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGtata                                                                                                                 <  2:30735/143‑1 (MQ=255)
                                      cacaGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAgga                                                                                                        >  1:318994/1‑142 (MQ=255)
                                              gTTGTGGTTGCATGGCACGAGTATAAATAATGGCTGCTACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAAGAACAAAAGCCCTGTGTATATATCAAGGCCAACaaat                                                                                                <  1:718799/143‑1 (MQ=255)
                                                   ggTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCAAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAAtt                                                                                           <  1:592312/143‑1 (MQ=255)
                                                   ggTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCAAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAAtt                                                                                           <  1:99820/143‑1 (MQ=255)
                                                    gTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAAttt                                                                                          >  1:391671/1‑143 (MQ=255)
                                                    gTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAAAACAAATGCCAGAGCAAGCAATAAAAGAAGCACTGTGGATATATCACGGCCAACAAAAAAAatt                                                                                          >  2:871693/1‑143 (MQ=255)
                                                     ttGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGAACTGAGAGAAAACCTGTGTATATATCAAGGCCAACAAATAAAATTa                                                                                         >  2:982761/1‑143 (MQ=255)
                                                      tGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAATAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAg                                                                                        <  2:973603/143‑1 (MQ=255)
                                                         aTGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCat                                                                                     <  1:213381/143‑1 (MQ=255)
                                                           ggCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCatat                                                                                   <  2:152199/143‑1 (MQ=255)
                                                             cACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGCACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAGATAAATATAGCatatat                                                                                 >  2:543316/1‑143 (MQ=255)
                                                               ataGTAGAAATAACGACTGACACCGAATTCGTCGAATCATGAAAACCATTGATATCTACATAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCTCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTa                                                                               <  2:545181/141‑1 (MQ=255)
                                                                            cGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCAttt                                                                  >  1:370811/1‑143 (MQ=255)
                                                                             ggCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGCAGAACAACACAACTGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATCTATTAGATTACGCGAtgtt                                                                 >  2:53927/1‑143 (MQ=255)
                                                                               cTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGa                                                               <  1:261049/143‑1 (MQ=255)
                                                                               cTGCCAACGCATTCGCAGTGTCATGAAAACCATTGATTTCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATGTAGCATATATTAGATTACGCCATTTTGa                                                               <  2:635236/143‑1 (MQ=255)
                                                                                tGCCACCGCATTCGCCGGGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGaa                                                              <  2:127811/143‑1 (MQ=255)
                                                                                    aCCGCATGCGCCGTGTCATGAAAACCATTGCTGGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAGAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTa                                                          <  1:80525/143‑1 (MQ=255)
                                                                                         aTTCGTGGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAAc                                                     <  2:862461/143‑1 (MQ=255)
                                                                                           tCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACgg                                                   >  1:681621/1‑143 (MQ=255)
                                                                                             gCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATCACAAAAGCCCTGTGTATATATCAAGGCCAACAAATACGTTTAGGCTAGATCAGATTACGCCATTTCGGAGTTACGAAgggct                                                 >  1:158765/1‑141 (MQ=255)
                                                                                                   tCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTa                                           >  2:236975/1‑143 (MQ=255)
                                                                                                    cATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTAt                                          <  1:505308/143‑1 (MQ=255)
                                                                                                      tGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCa                                        <  1:986725/143‑1 (MQ=255)
                                                                                                        aaaaCCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGt                                      >  1:881325/1‑143 (MQ=255)
                                                                                                          acccATTGATTGCGTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGa                                    <  1:260867/141‑1 (MQ=255)
                                                                                                               ttgattgCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTa                               >  1:856138/1‑143 (MQ=255)
                                                                                                                    tgCTTCGTAGAACAATACAAATGCCAGAGCATGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGc                          <  2:1014065/143‑1 (MQ=255)
                                                                                                                      cTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCAt                        >  1:603629/1‑143 (MQ=255)
                                                                                                                            aGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCaa                  >  2:171240/1‑143 (MQ=255)
                                                                                                                             gAACAACACAAATGCCAGAGCAAGCAATAACAAAGGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCaaa                 <  2:580828/143‑1 (MQ=255)
                                                                                                                                  acacAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAACGAGaa            >  2:793222/1‑143 (MQ=255)
                                                                                                                                       aaTGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGTTTTTATTATCGGACGCATTATCAGTGACTTTAACGTCATGGGCAAAGTGAAATATc       >  2:635278/1‑143 (MQ=255)
                                                                                                                                       aaTGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAAAGTGAAATATc       <  2:552795/143‑1 (MQ=255)
                                                                                                                                          gCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAAAGTGAAATATCAtt    <  1:982761/143‑1 (MQ=255)
                                                                                                                                           ccAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTAAGCCATTTTGAATTTACGAACGGACGCATTATAAGTGACTTTAACGGTATGGGGAAAGTGAAATATCAttt   >  2:491175/1‑143 (MQ=255)
                                                                                                                                            cAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAAAGTGAAATATCAtttt  <  2:681621/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
AATAACACGCCAAAAAAGTTAAAAAATGCCGCCATCACCACAGCAAGTTGTGGTTGCATGGCACGAGTATAAATAACGGCTGCCACCGCATTCGCCGTGTCATGAAAACCATTGATTGCTTCGTAGAACAACACAAATGCCAGAGCAAGCAATAACAAAAGCCCTGTGTATATATCAAGGCCAACAAATAAATTTAGCATATATTAGATTACGCCATTTTGAATTTACGAACGGACGCATTATCAGTGACTTTAACGGCATGGGCAAAGTGAAATATCATTTT  >  NC_000913/3136172‑3136454

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: