Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 950,098 T→A 9.3% intergenic (+885/‑270) prsA1 → / → yqiG peptidylprolyl isomerase/flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730950,0980TA 9.3% 134.8 / 5.4 43intergenic (+885/‑270)prsA1/yqiGpeptidylprolyl isomerase/flavin reductase
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base A (2/2);  total (23/20)
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.81e-01

CTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTGTTTTT  >  CP000730/949958‑950237
                                                                                                                                            |                                                                                                                                           
cTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTttat                                                                                                                                             >  1:149169/1‑141 (MQ=255)
 ttACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTttatt                                                                                                                                            >  2:190342/1‑141 (MQ=255)
       agccggATAAAACAGTCAAAGCGAGCAATTAGACACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATAc                                                                                                                                      <  2:371272/137‑1 (MQ=39)
         cTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTa                                                                                                                                    >  2:227475/1‑141 (MQ=255)
          tGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                                                                   <  1:478565/141‑1 (MQ=255)
          tGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTaa                                                                                                                                   <  2:448323/141‑1 (MQ=255)
            gATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTAATAAAATTTAATAGATACTAAtc                                                                                                                                 <  1:63380/141‑1 (MQ=255)
                aaaCAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCtt                                                                                                                             <  2:133861/141‑1 (MQ=255)
                 aaCAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCttt                                                                                                                            >  1:31930/1‑141 (MQ=255)
                              aGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATaa                                                                                                               >  1:336917/1‑141 (MQ=255)
                               gCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTACTCTCTTTTTGTCTACGATAAc                                                                                                              >  1:3954/1‑141 (MQ=255)
                                 aaTTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGt                                                                                                            >  1:11792/1‑141 (MQ=255)
                                     cGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAata                                                                                                        >  1:84006/1‑141 (MQ=255)
                                         aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                                                                    >  1:239459/1‑141 (MQ=255)
                                         aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                                                                    <  1:159500/141‑1 (MQ=255)
                                               aaaGGAACAGATAGATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                                                                              <  2:421189/141‑1 (MQ=255)
                                               aaaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTc                                                                                              <  2:199173/141‑1 (MQ=255)
                                                aaGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCt                                                                                             >  1:210497/1‑141 (MQ=255)
                                                  ggAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGTTTACGTAATATGATTGACTCTAt                                                                                           >  2:408887/1‑141 (MQ=255)
                                                       aGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTaca                                                                                      >  1:399244/1‑141 (MQ=255)
                                                        gATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                                                                                     >  1:274964/1‑141 (MQ=255)
                                                        gATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTacac                                                                                     <  1:34794/141‑1 (MQ=255)
                                                           aaaTCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTa                                                                                  <  2:255181/141‑1 (MQ=255)
                                                              tCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaa                                                                               <  1:20895/141‑1 (MQ=255)
                                                               cGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACaaa                                                                              <  2:503262/141‑1 (MQ=255)
                                                                 tACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATg                                                                            <  2:23564/141‑1 (MQ=255)
                                                                       aaaTCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTTTCTACGTTAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTaa                                                                      >  2:190983/1‑141 (MQ=255)
                                                                             tACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGAc                                                                <  1:161312/141‑1 (MQ=255)
                                                                                      tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt                                                       <  2:185701/141‑1 (MQ=255)
                                                                                      tCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCAtt                                                       >  2:455201/1‑141 (MQ=255)
                                                                                         gaaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATc                                                    >  2:444097/1‑141 (MQ=255)
                                                                                          aaatgaaaGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCt                                                   <  2:101838/141‑1 (MQ=255)
                                                                                                       ttAGACTATAAAAGAGATTAGTATCTATTAGATTTTATTAGATACTAATCTCTTTTTGTCTACTATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTtagataggt                                      >  2:124787/1‑139 (MQ=255)
                                                                                                           aCTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGt                                  >  2:387193/1‑141 (MQ=255)
                                                                                                                aaaaGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGAtt                             >  1:378851/1‑141 (MQ=255)
                                                                                                                   agagATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTTTTTATGATACCGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATTTCCATTCTCTTTGTTATATTTAATCGTTGATTTgg                          >  2:412009/1‑140 (MQ=255)
                                                                                                                      gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAat                       <  2:274964/141‑1 (MQ=255)
                                                                                                                      gaTTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATAGTTGATTTGCAat                       >  2:69487/1‑141 (MQ=255)
                                                                                                                      gaTTAGTATCTATTAAAGTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAat                       <  1:27056/141‑1 (MQ=255)
                                                                                                                          aGTATCTAATAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTg                   >  1:498769/1‑141 (MQ=255)
                                                                                                                           gTATCTATTAAATTTTATTAGATACTAATCTCTTTTTTTTTATTATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGt                  >  2:314988/1‑141 (MQ=255)
                                                                                                                            tATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTCCAAGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTa                 <  1:249497/141‑1 (MQ=255)
                                                                                                                               ctaatAAAATTTAATAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTAtgt              >  2:410638/5‑141 (MQ=255)
                                                                                                                                tATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTAtgtg             <  1:155969/141‑1 (MQ=255)
                                                                                                                                      aTTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAAAATCGTTGATTTGCCATATTGTATGTGGATTTg       >  2:347284/1‑141 (MQ=255)
                                                                                                                                          tattaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTGtttt   >  1:451986/1‑141 (MQ=255)
                                                                                                                                           attaGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTGttttt  >  1:444470/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTGTTTTT  >  CP000730/949958‑950237

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