Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,616 G→C 7.8% A856P (GCA→CCA)  ydbH → PF11739 family protein YdbH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,445,6160GC 7.8% 140.4 / 4.0 53A856P (GCA→CCA) ydbHPF11739 family protein YdbH
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base C (2/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCT  >  NC_000913/1445474‑1445758
                                                                                                                                              |                                                                                                                                               
tctcGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAgg                                                                                                                                                 <  2:375608/143‑1 (MQ=255)
 ctcGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGc                                                                                                                                                >  2:363428/1‑143 (MQ=255)
     gTTTATTACCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATgc                                                                                                                                            <  1:524964/143‑2 (MQ=255)
        tattaACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTg                                                                                                                                         <  1:167264/143‑1 (MQ=255)
          ttaACCATGCAGACCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                       <  1:68467/143‑1 (MQ=255)
            aaCAATGCAGGCCAAAACTACAGGTACCAGTCGCGTGGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAga                                                                                                                                     <  1:415985/143‑2 (MQ=255)
             aCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCa                                                                                                                                    >  2:86348/1‑143 (MQ=255)
             aCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCa                                                                                                                                    >  1:108513/1‑143 (MQ=255)
             aCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCa                                                                                                                                    >  1:321660/1‑143 (MQ=255)
                       ccAACATTACAGGTACCTGTCGCGTTGATGGTAAAAGTGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATATTCTTCAGTCATGGCTGGAGCAGAACGCACGt                                                                                                                          <  2:447584/143‑1 (MQ=255)
                         aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCt                                                                                                                        >  1:237610/1‑143 (MQ=255)
                          aCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGCAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCACGGCTGGAGCAGAACGCACGTCTg                                                                                                                       >  2:547056/1‑143 (MQ=255)
                              tACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACACGGTAGCTGGGTTTTCGCATGAACGATAATCTCCAGGCATGGCTGGAGGAGAACGCACGTCTGCCgg                                                                                                                   >  2:241717/1‑143 (MQ=255)
                                aaGGTGCGAGTCGCGTTGATGGTAAAAGCGGTATGGTAAACCTTAACTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGa                                                                                                                 <  2:512360/142‑1 (MQ=255)
                                  tggaCCAGGCGCGTTCATGCTAAAAGAGCGACGGTACACCTTAACTCCCATCATGAAGTGGATATTTTTACTCTGTCGCGCCGTTTACGCTTTGGCTATAATCGCGAAGCATGTCTGGAGCAGAACGCACGTCTGCCGGGaaa                                                                                                               <  2:204076/140‑1 (MQ=255)
                                    tACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATg                                                                                                             >  2:449922/1‑143 (MQ=255)
                                     aCCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                            <  1:240541/143‑1 (MQ=255)
                                     aCCAGGCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                            <  1:501694/143‑1 (MQ=255)
                                         gTCGCGTTGATAGTAAGAGCGCTACGGTCAGTCTCGTGGTTCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCTATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGTGAAATGACTGt                                                                                                        <  1:98989/143‑1 (MQ=255)
                                             cgTTGATGGTAAAGGCGGTACGGTAAACCTTGATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                    <  2:514717/143‑1 (MQ=255)
                                               ctGAGGGTGAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCGCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                  <  2:507122/142‑1 (MQ=255)
                                                 gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCCGCCGGGAAATGACTGTCCGCAAgg                                                                                                >  1:224918/1‑143 (MQ=255)
                                                  agcgTAAAAGCGGTACGGTAAACATTAATTACCATCATGAATAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTGGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACTTCTTTCGTGAAATGACTGTCCGCAAGGa                                                                                               <  2:86296/140‑1 (MQ=255)
                                                    ggTATACGCGTTACGGTAAACTTTAATTACCATCATGAAGATAATATTTTTATGTTGTTTTTTTGTTTTCTCTTTGGCGATAAGCTACAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaa                                                                                             <  1:48487/143‑1 (MQ=255)
                                                     gTAAATGCGGTACGGTAAACCTTAATTACCAGCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                            <  1:291736/143‑1 (MQ=255)
                                                      tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGCAATGACTGTCCGCAAGGAAAag                                                                                           >  2:395475/1‑143 (MQ=255)
                                                         aaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAagaga                                                                                        >  1:457968/1‑142 (MQ=255)
                                                          atgggtacggtaACACTTAATTACCATCATGCAGAGAATATTTTGACGCGGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCTGGAAATGACTGTCCGCAAGGAAAAGAgtg                                                                                       <  2:201636/140‑1 (MQ=255)
                                                          aGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGCAATGACTGTCCGCAAGGAAAagcggg                                                                                       >  1:419942/1‑139 (MQ=255)
                                                          aGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAagaggg                                                                                       >  1:356514/1‑141 (MQ=255)
                                                             ggtacggtaAGCCTGTATTGCCGTGATGGAGGGAAGATTTTTACGCTGTGGCGCAGGGTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGTCTGTCCGCAAGGAAAAGAGTGTGa                                                                                    <  2:126671/143‑1 (MQ=255)
                                                                  ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGCATAGCCCCAAGGAAAAGACTGTCAGGa‑aa                                                                              >  2:259287/1‑143 (MQ=255)
                                                                        ccTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAAt                                                                        >  2:189960/1‑143 (MQ=255)
                                                                            aaTTACCGTCATAAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGACATGGCTGGAGCAGAACGCACGTCTGCCGAGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGaaa                                                                    <  2:542381/143‑1 (MQ=255)
                                                                              ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAAt                                                                  <  2:420949/143‑1 (MQ=255)
                                                                              ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAAt                                                                  >  1:536332/1‑143 (MQ=255)
                                                                                             gaATATTTTTACGCTGTGGCGCAGTTTACGCGTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGtt                                                   <  1:453617/143‑1 (MQ=255)
                                                                                                    tttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCa                                            <  1:408348/143‑1 (MQ=255)
                                                                                                     ttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGAGTGAGGA‑AAAACAATGAAAATTTTACTGGATGCGTTGACGTCAt                                           >  2:272412/1‑143 (MQ=255)
                                                                                                       aCGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCt                                         >  2:460437/1‑143 (MQ=255)
                                                                                                           tgtgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                     <  1:453119/143‑1 (MQ=255)
                                                                                                             tgGCGCAGGTTACGAGTGGGGGATAATATCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGACGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATg                                   <  2:518153/143‑1 (MQ=255)
                                                                                                                     tttACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGACGAG‑GGGTACTTATAGAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTgg                           >  1:251510/1‑143 (MQ=255)
                                                                                                                          ttttgggCGTTAATGTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAGTGACTGTCCGCAAGGAAGAGTGTGTGAGGA‑AAAACGATGAAAATTTTACTGGCTGCGTTGACTTCATCTTTTATGCTGGTTGGCTGTa                      <  2:194048/138‑1 (MQ=255)
                                                                                                                             ttggataTAATCTCCAGGCGTGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGc                   <  2:321660/137‑1 (MQ=255)
                                                                                                                             ttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGc                   <  1:86348/143‑1 (MQ=255)
                                                                                                                             ttGGCGATAATCTCCAGGCAGGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGc                   <  2:108513/143‑1 (MQ=255)
                                                                                                                                gctaTAATCTCCAGGCATGGCGGGAGGGGGGTGCACGTTTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTc                <  1:424086/140‑1 (MQ=255)
                                                                                                                                gCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAAAGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGTCTGTACGCCTc                >  1:279563/1‑143 (MQ=255)
                                                                                                                                 cGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTCCGCCTca               >  2:427228/1‑142 (MQ=255)
                                                                                                                                      tgctcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCAttc          >  2:365984/3‑142 (MQ=255)
                                                                                                                                        ctcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa        <  1:461586/143‑1 (MQ=255)
                                                                                                                                         tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCCGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTCAGACTCGCATTGAAg       >  2:103002/1‑143 (MQ=255)
                                                                                                                                            aGGCATGGCTGGCGCAGAACGCACGACGGCCGCGAAATGACTGTCCGCCAGGAAAAGTGTGTGAGGA‑AAAACACTGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCg    <  2:237610/143‑1 (MQ=255)
                                                                                                                                              gCATGGCTGAAGCAGAACGTACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAACATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCt  <  1:197783/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGA‑AAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCT  >  NC_000913/1445474‑1445758

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