Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 3,725,237 A→G intergenic (‑69/‑460) U712_19205 ← / → U712_19210 Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,2370AG84.6% 124.4 / ‑1.3 52intergenic (‑69/‑460)U712_19205/U712_19210Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (34/10);  minor base C (1/7);  total (35/17)
Fisher's exact test for biased strand distribution p-value = 9.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAC  >  CP006881/3725109‑3725353
                                                                                                                                       |                                                                                                                     
cGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGaaaa                                                                                                               >  6:28077/1‑138 (MQ=255)
         gTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAt                                                                                                    >  4:63699/1‑139 (MQ=255)
          ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTGTTTTTGGCAGGCCGCCTGTTTAGCAGCCCGTTTTAGCAGCCGCAGCCCGCAAGAAAAGACGGCCGCtt                                                                                                   >  2:383935/1‑139 (MQ=255)
          ttCATCAAAAAGACAATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCGTTTTGGCCGGCAGCCTATTTAGCAGCCCGTATTCTCCCCCGCAGCCCCCAATAAAAGACGGcctcag                                                                                                   >  4:112230/1‑135 (MQ=255)
               caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCATCCCCTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                              >  3:340290/1‑139 (MQ=255)
                  acacaCCATCCACACGCTCAAATATACAAGCAAAGGTTTAAAAGAAAAACATCCCTCCGTTTCTTTTTGGCAGGCAGCATTTTTAGCAGCCCGTTTTCTCAGCCGCCGCCCCAAATAAAACACGGCCGAtttcttttct                                                                                           <  6:252234/135‑1 (MQ=255)
                  aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                                           >  7:216906/1‑139 (MQ=255)
                  aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcgttta                                                                                            >  5:347805/1‑137 (MQ=255)
                      aCCATCCACACCCGTAAAGATACACCAACACTTTTCTTTGAAAAACATCCGTCCTCTTTTTTTTGGAAGGCAGCCTTTTACGCGGCGCGTTTGCTCAGCCCCCGCCCCCACGAAAACTCGGCCGATTTCTTTTCTCGCa                                                                                       <  8:49583/139‑1 (MQ=255)
                      aCCATCCACACCCGAAACGATACAACCAACGGTTTCTTACAAAAACACCCCTCCGCTTCTTTTGGGCAGCCAGCCTTTTTCCCCCCCCGTTTTCTCCCCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                       <  8:49577/139‑1 (MQ=255)
                            cacaCCCTTAAATATAAAAGCAATGGTTTTATAAAAAAACATCCGTCGGTTTTTTTTTGGCGGGCTGGTTTTTTAGCCGCCGGTTTTCTCCCCCGCCGCGGCAATAAAAACACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                 <  4:340290/139‑1 (MQ=255)
                            acagcccGAAAATAAACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCCGCCCGCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                 <  3:89543/135‑1 (MQ=255)
                                   gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTATCAGCCGCAGCCCGCAAGAAAAGACGGCCTATTTATTTTGTCGCAAACTAGCggagat                                                                          >  5:244016/1‑135 (MQ=255)
                                    tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACTAGCGGctc                                                                           >  4:378060/1‑137 (MQ=255)
                                      aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCATCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTat                                                                        >  3:345310/1‑137 (MQ=255)
                                           aCAAGCAAACGTTTCATAGAAAAACGTCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTGTGGCAAACAACCGGCGATAGCcaaa                                                                   >  1:203756/1‑138 (MQ=255)
                                           aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCactgt                                                                  >  4:89543/1‑135 (MQ=255)
                                           aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                   >  3:400520/1‑138 (MQ=255)
                                           aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACAAGCGCCTCGAGCaaaa                                                                   >  6:132349/1‑138 (MQ=255)
                                           aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACAAGCCGCTCTAGCaaaa                                                                   >  4:230857/1‑138 (MQ=255)
                                             aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTGTCCGAAACTAGCGGCGAGAGCAAAAca                                                                 >  1:122685/1‑137 (MQ=255)
                                             aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGATTTCGTTTCTCGCAAACGAGCGGCTCTAGCAAAACTc                                                                >  8:47171/1‑139 (MQ=255)
                                             aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                                >  4:403585/1‑123 (MQ=255)
                                             aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                                <  3:403585/123‑1 (MQ=255)
                                             aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTACCAAACTACCGGCTCAAGCAAAACt                                                                 >  2:41973/1‑138 (MQ=255)
                                               gaaaaCGTATCAAACAAAACCCTCCCTCCGTTTATTTCGGCCACGCAGCCTTTCTATCCGCCGGTTCTCTCAGCCGCAGCCCCCATCAACAGACGCCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAACCTCAt                                                              <  8:323905/137‑1 (MQ=255)
                                               gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACTCAt                                                              >  6:399406/1‑139 (MQ=255)
                                               gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                              <  7:68126/139‑1 (MQ=255)
                                               gCAAACGTTTCATAGAAAAACATCCCTCCGCTTATTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGCCTCTAGCAAAACTCAt                                                              >  5:172020/1‑139 (MQ=255)
                                                 aaaCTTTCCAAATAAAAACATCCCTCCGCTTCTTTTTGGCAGCCAGCCTTTTAAGCAGCCCGTTTTCTCAGCCGCCGCCCGAAACAAAAGACGGCCGATTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                            <  8:216906/139‑1 (MQ=255)
                                                 aaaCGTTTCTTATAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAGCCTTTTTAGCAGCCCTTTTTCTCCGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                            <  6:347805/139‑1 (MQ=255)
                                                   ccGTTTCATACAAAAACATCCCTCCGCTTTTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCGCAGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  1:41973/138‑1 (MQ=255)
                                                   ccGTCTCATTGAAAAACAACCCTCGGCTGCTTTTTGGCCGGCAGCCTTTTTTGCCGCCCTTTTGCACCCCCGCCCCCCGAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  4:77369/138‑1 (MQ=255)
                                                   acttttCATACAAAAACACCCCTCCGCTTCTTTTTGGCTGGCAGCATTTTTAGCAGCCCTTTTCCTCACCCGCCCCCCCAAATAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  5:377912/136‑1 (MQ=255)
                                                   aCGTTTCATATAAAAACACCCCTCCGTGTATTTTTGGCAGGCTGCCTTTTATGCCGCCCGTTTGCTCCGCCGCTCCCCGCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  4:93119/139‑1 (MQ=255)
                                                   aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGCTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAcg                                                          >  5:115785/1‑138 (MQ=255)
                                                   aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACGCATGAct                                                          >  7:368816/1‑139 (MQ=255)
                                                   aCGTTTATTACAAAAACACCCCCCCGCTTGTTTTTGCCCGCCCGCCTTTTTTGCCCCCGGTTTTCCCCGCCCCACCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  6:82942/139‑1 (MQ=255)
                                                    cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGATTTCTCAGCCGCAGCCCGCAAGAAAAAACGGCCCAGTACTATAAACGCAAACTAACCCCTCTAGCAAAACTAAAGAct                                                          >  8:393387/1‑138 (MQ=255)
                                                    cGTGTCATACAAAAACACCCCCCCGCGTCTTTTTGGCAGGCCCCCTTTTTCCCCGCCCGGTTGCACCCCCGCCGCCCGCAATAAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  2:35259/139‑1 (MQ=255)
                                                               aaaaCACCCCCCGGCTTCTTTTGGCCCGGCCGGTTTTTAAGCCCCCCGTTTACACCGCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                              <  1:121942/139‑1 (MQ=255)
                                                                        cTCCGCTTCTTTTTGGCAGGCCGCCTTTTTAGCAGCCCGTTTTCTAAGCCGCAGCCCGCCAGAAAAGACGGCCGAGTTGTTTTGTCGCAAACATGCGGCGCGAGCAAAACTCATTCATCTATGAGATACCGGGGaaa                                       >  6:262288/1‑137 (MQ=255)
                                                                        cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGAGTTCGTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCAGGACTCTATGATAGACCGTGGaaa                                       >  7:165254/1‑137 (MQ=255)
                                                                        cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                                       >  1:358022/1‑137 (MQ=255)
                                                                        cTCCGCTTCTTTTTGGCAGGCAGCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTAATTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTATATGATAGACCGTGGaaa                                       >  3:318368/1‑137 (MQ=255)
                                                                                           gcagCCTTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat                  >  3:57848/1‑139 (MQ=255)
                                                                                                    ttttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGCCAAACACAGGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACCATATGattca            >  2:269257/1‑134 (MQ=255)
                                                                                                               agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATTTGAGTGAAAAAAcc     >  5:380102/1‑138 (MQ=255)
                                                                                                               agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc     >  8:312038/1‑138 (MQ=255)
                                                                                                               agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCGATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc     >  1:337899/1‑138 (MQ=255)
                                                                                                               agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTAGCAAACTAGCGGCTCTAGCAAAACTCCTGACTCTATGTTAGACAGTGGAAAGAGAATTCTTTTTGACAAAATGAGTGTAAAAAcc     >  5:40408/1‑138 (MQ=255)
                                                                                                                 cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc  >  4:215313/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                     
CGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTTTTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAC  >  CP006881/3725109‑3725353

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