Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 169,820 T→C 11.6% intergenic (+784/+2903) PP_RS08750 → / ← fhuA_1 tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440169,8200TC11.6% 155.3 / 12.2 60intergenic (+784/+2903)PP_RS08750/fhuA_1tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base C (2/5);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTT  >  KT2440/169687‑169923
                                                                                                                                     |                                                                                                       
cctccAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTa                                                                                                        <  3:303286/135‑1 (MQ=32)
       tgCTAAATACTACTGACTGACCGATAGAGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGACCCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACt                                                                                                 <  1:44258/134‑1 (MQ=12)
        gCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                <  1:155863/135‑1 (MQ=32)
        gCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                >  4:192127/1‑135 (MQ=32)
        gCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                <  6:290134/135‑1 (MQ=32)
        gCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTg                                                                                                <  8:110241/135‑1 (MQ=32)
             aTACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTAc                                                                                           <  1:226213/135‑1 (MQ=32)
               actactGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCt                                                                                         >  3:282238/1‑135 (MQ=32)
               actactGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCt                                                                                         >  8:381740/1‑135 (MQ=32)
                      actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                   >  2:349743/1‑134 (MQ=32)
                      actgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGta                                                                                   >  3:225761/1‑134 (MQ=32)
                       ctgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGta                                                                                   <  2:354820/133‑1 (MQ=11)
                       ctgacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGtat                                                                                  <  8:102352/134‑1 (MQ=32)
                         gacCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAAt                                                                               >  3:308206/1‑135 (MQ=32)
                            cGATAGTGAACCAGTACAGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                            <  4:225009/135‑1 (MQ=17)
                              attgTGAACCAGTACCGTGCGGGAAAGCCGAAACCAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  3:263328/132‑1 (MQ=11)
                              aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATggg                                                                            >  7:314039/1‑133 (MQ=32)
                              aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  8:226311/135‑1 (MQ=32)
                              aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  6:237444/135‑1 (MQ=32)
                              aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  3:310179/135‑1 (MQ=32)
                              aTAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  1:294378/135‑1 (MQ=32)
                                      ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                  <  6:205512/135‑1 (MQ=11)
                                      ccAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACtt                                                                  <  5:367434/135‑1 (MQ=11)
                                       cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          >  4:218308/1‑126 (MQ=11)
                                       cAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCACGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                                                          <  3:218308/126‑1 (MQ=11)
                                            ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                              >  3:320734/1‑133 (MQ=32)
                                            ccGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTatat                                                              >  5:165627/1‑133 (MQ=32)
                                               tGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGt                                                         <  4:375783/135‑1 (MQ=32)
                                                 aGGGAAAGGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGta                                                        >  8:70474/1‑133 (MQ=11)
                                                 aGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  7:70047/1‑135 (MQ=32)
                                                 aGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  6:357019/1‑135 (MQ=32)
                                                  gggAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTgg                                                       >  5:91997/1‑134 (MQ=32)
                                                  gggAAAGGCGAAAAGAACCAAGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAATCAGTCGGAGCCTACATTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCTACTAATATTCAGTgg                                                       >  8:21603/1‑134 (MQ=11)
                                                    gAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGAGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGa                                                    >  6:263115/1‑135 (MQ=17)
                                                       aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGTCGAGc                                                  >  8:271791/1‑134 (MQ=17)
                                                       aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                                                  >  4:222878/1‑134 (MQ=32)
                                                       aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGc                                                  >  8:169891/1‑134 (MQ=32)
                                                       aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                                                 <  5:3209/135‑1 (MQ=32)
                                                       aGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGGACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCt                                                 <  4:63021/135‑1 (MQ=17)
                                                         gCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTa                                               >  3:47477/1‑135 (MQ=32)
                                                          cGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTaa                                              <  1:349743/135‑1 (MQ=32)
                                                               aGAACCCCTGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGaa                                         <  7:11072/135‑1 (MQ=17)
                                                                gAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGa                                          >  6:132988/1‑133 (MQ=32)
                                                                      cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGa                                  >  1:307201/1‑135 (MQ=32)
                                                                      cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGa                                  >  7:139107/1‑135 (MQ=32)
                                                                      cGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGa                                  >  5:391824/1‑135 (MQ=32)
                                                                              ggAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAg                           >  3:276311/1‑134 (MQ=32)
                                                                                aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCg                         >  7:242725/1‑134 (MQ=32)
                                                                                aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  8:314039/135‑1 (MQ=32)
                                                                                aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  7:381740/135‑1 (MQ=32)
                                                                                aGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  8:142949/135‑1 (MQ=32)
                                                                                aGTGAAATAGAACCTGAAACCGTGGGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGGCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGa                        <  5:211693/135‑1 (MQ=17)
                                                                                  tGAATTAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAATCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGaaa                      <  2:309222/135‑1 (MQ=17)
                                                                                     aaTAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGc                    <  3:286219/134‑1 (MQ=32)
                                                                                        aGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCAAGt                >  2:178701/1‑135 (MQ=17)
                                                                                            ccTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTa            <  4:225761/135‑1 (MQ=32)
                                                                                                 aaCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTAcc                                                                                          >  3:233675/1‑52 (MQ=35)
                                                                                                 aaCCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTAcc                                                                                          <  4:233675/52‑1 (MQ=35)
                                                                                                     gtgtGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGt   <  3:192127/135‑1 (MQ=32)
                                                                                                      tgtgCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCTAGATTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGtt  <  4:308206/135‑1 (MQ=14)
                                                                                                                                     |                                                                                                       
CCTCCAAGGCTAAATACTACTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGAGAGGGGAGTGAAATAGAACCTGAAACCGTGTGCGTACAAGCAGTGGGAGCCTACTTTGTTAGGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCAGTGGCGAGCTTAACCGAATAGGGGAGGCGTAGCGAAAGCGAGTCTTAATAGGGCGTT  >  KT2440/169687‑169923

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: