Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 25,865 C→T 10.8% H14Y (CAC→TAC)  fkpB → FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927325,8650CT10.8% 95.2 / 3.4 37H14Y (CAC→TAC) fkpBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/2);  total (19/18)
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.00e-01

CTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGAAAGTGGGCGATAA  >  CP009273/25723‑26004
                                                                                                                                              |                                                                                                                                           
cTGACGCATGTTATCCGGCCCACAGATTGCTGCGAAATCGTAGGGCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATACAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCa                                                                                                                                       <  1:527/149‑1 (MQ=255)
   aCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGTACt                                                                                                                                          >  2:59546/1‑143 (MQ=255)
   aCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGTACt                                                                                                                                          >  2:59478/1‑143 (MQ=255)
        tCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaaa                                                                                                                               <  1:160036/149‑1 (MQ=255)
            aTCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTACATCTCAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcg                                                                                                                           <  1:334067/149‑1 (MQ=255)
                 gCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTc                                                                                                                            <  2:441730/143‑1 (MQ=255)
                    tACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGc                                                                                                                   >  1:347570/1‑149 (MQ=255)
                              tgcGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGcacca                                                                                                               >  2:314712/1‑143 (MQ=255)
                                cGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAg                                                                                                       >  1:226427/1‑149 (MQ=255)
                                  aaaTCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTc                                                                                                     >  1:369498/1‑149 (MQ=255)
                                        tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTc                                                                                                     >  2:435898/1‑143 (MQ=255)
                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGcaa                                                                                            <  1:144726/149‑1 (MQ=255)
                                                     gCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTATGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTAGACGATGGCCCCCCCGCCGGGTGTACCCGCAACAACGGGaaa                                                                                  >  1:145470/1‑149 (MQ=255)
                                                           actccgcTTCCGGCAAAAAACCTTAAACATAAGAGCAAACCGGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACGTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGgcg                                                                            <  1:259503/146‑1 (MQ=255)
                                                             gccgcATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGGCTACCCGCAACAACGGGAAACCGGCGCt                                                                          >  1:165921/1‑149 (MQ=255)
                                                              ccgcATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTg                                                                         >  1:300017/1‑149 (MQ=255)
                                                                       gCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGTACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTg                                                                <  1:66963/149‑1 (MQ=255)
                                                                       gCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGTACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTg                                                                <  1:67037/149‑1 (MQ=255)
                                                                       gCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTg                                                                >  1:527333/1‑149 (MQ=255)
                                                                       gCAAAAATCCCTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAACACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTg                                                                <  1:158433/149‑1 (MQ=255)
                                                                               ccTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGc                                                        <  1:399766/149‑1 (MQ=255)
                                                                               ccTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGc                                                        <  1:524089/149‑1 (MQ=255)
                                                                                  tAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTtc                                                     <  1:145028/149‑1 (MQ=255)
                                                                                    aaTATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATg                                                         <  2:203205/143‑1 (MQ=255)
                                                                                            agCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCtt                                                 >  2:458776/1‑143 (MQ=255)
                                                                                                 aCCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTGCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCt                                      >  1:22671/1‑149 (MQ=255)
                                                                                                    tGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGACgg                                         >  2:84703/1‑143 (MQ=255)
                                                                                                     gCAGGTCTGAATCTGTACAGAGCAATGGCGCCGTCCTGGTGCACTTCACGCTACAACTCGACGATTGCACCACCGCCGAGTCTACCCGCAACACCGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAggg                                        <  2:243283/143‑1 (MQ=255)
                                                                                                      cATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACATCGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGc                                       <  2:309085/143‑1 (MQ=255)
                                                                                                          tCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGCCGCTGTTCCGCCTGGGTGATGCTTCTGTTTCTGGAGGGCTGGa                                   >  2:395409/1‑143 (MQ=255)
                                                                                                           cTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGACGGGCTGGAg                                  <  2:124322/143‑1 (MQ=255)
                                                                                                                         agCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTggg                    <  2:472303/143‑1 (MQ=255)
                                                                                                                          gCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGaaa             >  1:216732/1‑149 (MQ=255)
                                                                                                                                   ccGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCTGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGAGGATTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGAAAGTg          >  2:11522/1‑143 (MQ=255)
                                                                                                                                   ccGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGGAAGTGGGCGAt    >  1:309545/1‑149 (MQ=255)
                                                                                                                                      tCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGAAAGTGGGc       <  2:160470/143‑1 (MQ=255)
                                                                                                                                           gTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGAAAGTGCGCGGTaa  >  2:414219/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTAAACCGGCGCTGTTCCGCCTGGGTGATGCTTCTCTTTCTGAAGGGCTGGAGCAACACCTGTTGGGGCTGAAAGTGGGCGATAA  >  CP009273/25723‑26004

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