Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 950,093 T→A 8.8% intergenic (+880/‑275) prsA1 → / → yqiG peptidylprolyl isomerase/flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730950,0930TA 8.8% 145.7 / 5.3 46intergenic (+880/‑275)prsA1/yqiGpeptidylprolyl isomerase/flavin reductase
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base A (2/2);  total (26/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.96e-01

GAACGCTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTG  >  CP000730/949953‑950232
                                                                                                                                            |                                                                                                                                           
gAACGCTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAAt                                                                                                                                             >  2:452068/1‑141 (MQ=255)
    gCTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTtta                                                                                                                                         >  2:399665/1‑141 (MQ=255)
     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:159500/141‑1 (MQ=255)
                                              aCGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATAtgat                                                                                               >  1:239459/1‑141 (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)
                                                                                                                                            |                                                                                                                                           
GAACGCTTACCAATCTGGATAAAACAGTCAAAGCGAGCAATTCGAAACGACTAAAGGAACAGATAAATCGTACCGAAAATCATACAAAGGGTCTGAAATGAAAGTTTCTTAGACTATAAAAGAGATTAGTATCTATTAAATTTTATTAGATACTAATCTCTTTTTGTCTACGATAACGTAATATGATTGATTCTATTTACACGTACAAATGGTTTAAGGTGACATATCCATTATCTTTGTTAGATAGAATCGTTGATTTGCAATATTGTATGTGGATTTG  >  CP000730/949953‑950232

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: ‑