Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,687,976 A→G 7.5% T302A (ACA→GCA)  U712_08695 → Zinc metalloprotease rasP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,687,9760AG 7.5% 159.1 / 2.2 54T302A (ACA→GCA) U712_08695Zinc metalloprotease rasP
Reads supporting (aligned to +/- strand):  ref base A (17/32);  new base G (2/2);  total (20/34)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGATAACAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCATATATGACATGACAGACCAAGTGGCGAAAACAGGGATCGT  >  CP006881/1687838‑1688106
                                                                                                                                          |                                                                                                                                  
acacaaaaCAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCaa                                                                                                                                    <  1:367242/133‑1 (MQ=255)
   gaTAACAAAACGCTTCATATTTCGGGGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAg                                                                                                                                 <  8:54730/139‑1 (MQ=255)
    atacaaaaaCGCTTCATATATCGGTGACTCCGGAAGCTTTTAAAGATGAGAACAAAAACACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGt                                                                                                                                <  6:207557/134‑1 (MQ=255)
     tAACAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGtt                                                                                                                               >  5:184946/1‑139 (MQ=255)
     tAACAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATAGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGtt                                                                                                                               >  7:229401/1‑139 (MQ=255)
        cAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGTGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGAt                                                                                                                            <  2:132226/139‑1 (MQ=255)
             cGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCAc                                                                                                                       <  5:81821/139‑1 (MQ=255)
                ttCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGAt                                                                                                                            <  2:253156/131‑1 (MQ=255)
                ttCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGAt                                                                                                                            >  1:253156/1‑131 (MQ=255)
                   atatTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAg                                                                                                                  <  1:264709/138‑1 (MQ=255)
                   atatTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAg                                                                                                                  <  7:22815/138‑1 (MQ=255)
                   atatTTCGGTGACTCCGGAAGCTGTTAAAGAGGAGAACAAAAAAACACTCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAg                                                                                                                  <  5:302214/138‑1 (MQ=255)
                       ttCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAg                                                                                                                  <  4:273446/134‑1 (MQ=255)
                       ttCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAg                                                                                                                  >  3:273446/1‑134 (MQ=255)
                           gTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTaa                                                                                                         <  7:169511/139‑1 (MQ=255)
                              aCTCCGGATGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGAAGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTAAGGCGCGACATCACCAGTTGATGTCACCAAAGCCATTTTAAca                                                                                                       >  4:144719/1‑137 (MQ=255)
                              aCTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAAcc                                                                                                       >  1:56857/1‑138 (MQ=255)
                              aCTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAAcc                                                                                                       >  2:204693/1‑138 (MQ=255)
                                  cGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAAt                                                                                                    <  2:56857/137‑1 (MQ=255)
                                      acaagTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGc                                                                                              <  5:254256/135‑1 (MQ=255)
                                          gTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGTCGTACTCTCAGCGGTTGCTTACGGCGCGACAGCAACAGTTGAGGTCACCAAAGCCAAGTTATCCAATCTGAGCAAAt                                                                                          >  6:306199/1‑139 (MQ=255)
                                          gTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAAt                                                                                          >  4:322825/1‑139 (MQ=255)
                                          gTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAAt                                                                                          >  4:151542/1‑139 (MQ=255)
                                          gTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTACAAAGCGCCGACTGAAAAAGGCGAACTCTCATCGGTTGCTTACGGCGCGACAACAACAGTTGATGTCACCATAGCCATTTTACACAATAAGAGCaaa                                                                                           >  6:53295/1‑138 (MQ=255)
                                           ttAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTGGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAAtt                                                                                         <  6:352864/139‑1 (MQ=255)
                                                gATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTaa                                                                                    <  3:123329/139‑1 (MQ=255)
                                                     gaACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               >  7:192261/1‑139 (MQ=255)
                                                     gaACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               >  7:67706/1‑139 (MQ=255)
                                                     gaACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               <  7:99315/139‑1 (MQ=255)
                                                     gaACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               <  7:139945/139‑1 (MQ=255)
                                                     gaACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGGTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               <  6:370761/139‑1 (MQ=255)
                                                           aaaaaCAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAAtt                                                                          >  1:368308/1‑138 (MQ=255)
                                                            aaaaCAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTa                                                                        <  7:171400/139‑1 (MQ=255)
                                                               aCAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAAc                                                                     <  3:151542/139‑1 (MQ=255)
                                                                cAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACt                                                                    >  6:367921/1‑139 (MQ=255)
                                                                     ggACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGGTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGat                                                                <  5:53295/138‑1 (MQ=255)
                                                                                      tGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGt                                               <  6:184946/138‑1 (MQ=255)
                                                                                        cgcCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCgg                                            <  3:352442/139‑1 (MQ=255)
                                                                                        cgcCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCg                                             <  2:368308/138‑1 (MQ=255)
                                                                                        cgcCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACAAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCgg                                            <  8:229401/139‑1 (MQ=255)
                                                                                         gcCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGc                                           <  8:67706/139‑1 (MQ=255)
                                                                                          ccGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAAtt                                                                                         <  1:230782/92‑1 (MQ=255)
                                                                                          ccGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAAtt                                                                                         >  2:230782/1‑92 (MQ=255)
                                                                                            gACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCata                                        <  5:132070/139‑1 (MQ=255)
                                                                                            gACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCata                                        <  3:324089/139‑1 (MQ=255)
                                                                                                      ggCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCATATatgacatga                              <  3:322825/139‑1 (MQ=255)
                                                                                                           aCTCTCAGCGGTTGCTTACGGCGCGACATCAGCAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATg                                                             >  3:104521/1‑103 (MQ=255)
                                                                                                           aCTCTCAGCGGTTGCTTACGGCGAGACATCAGCAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATg                                                             <  4:104521/103‑1 (MQ=255)
                                                                                                                 aGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCATATATGACATGACAGACCAAGTg                   >  1:342005/1‑139 (MQ=255)
                                                                                                                 aGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGGCGGCATATATGACATGACAGACCaa                      <  5:367921/136‑1 (MQ=255)
                                                                                                                 aGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAAGCTGCGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCTTATATGACATGACAGACCCAGTg                   >  1:314220/1‑139 (MQ=255)
                                                                                                                   cGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGACATTTTAACCTATCTGAGCAAATTAGTAACAGGCGAATTTAATCTCGATATGCTGTCAGGTCCTGTCGGCATATATGACATGACAGACCAAGTGGc                 <  7:35433/139‑1 (MQ=255)
                                                                                                                                  cgACATCAGCAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               >  6:297503/1‑62 (MQ=255)
                                                                                                                                  cgACATCAGCAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGc                                                                               <  5:297503/62‑1 (MQ=255)
                                                                                                                                  cgACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCATATATGACATGACAGACCAAGTGGCGAAAACAGGGATCGt  <  1:204693/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                  
AGAGATAACAAAACGCTTCATATTTCGGTGACTCCGGAAGCTGTTAAAGATGAGAACAAAAAAACAATCGGACGTTTCGGTTCCTATGCGCCGACTGAAAAAGGCGTACTCTCAGCGGTTGCTTACGGCGCGACATCAACAGTTGATGTCACCAAAGCCATTTTAACCAATCTGAGCAAATTAGTAACAGGCCAATTTAAACTCGATATGCTGTCAGGTCCTGTCGGCATATATGACATGACAGACCAAGTGGCGAAAACAGGGATCGT  >  CP006881/1687838‑1688106

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