Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,778,338 C→A 12.1% *323Y (TAG→TAT arcC3 ← carbamate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,778,3380CA12.1% 118.4 / 9.2 50*323Y (TAG→TATarcC3carbamate kinase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base A (3/3);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

CTAATTTATTATTTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATCCTTCCA  >  CP000730/2778200‑2778477
                                                                                                                                          |                                                                                                                                           
cTAATTTATTATTTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTa                                                                                                                                           >  2:302432/1‑141 (MQ=255)
 tAATTTATTATTTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTAc                                                                                                                                          >  2:121480/1‑141 (MQ=255)
 tAATTTATTATTTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTAc                                                                                                                                          >  1:70130/1‑141 (MQ=255)
    ttgattggttCTAGCCAAAATATAGTTTTCTGTCATGAGTTAAACTGCTTCAAGCCTTGCAAAGTGTCAGCAGACAACTTTTTCCTTCCCCCCCAAGGAGAAACTTATTTAATGTGTTGTCGCGTGATAGAATGCTaacgg                                                                                                                                       >  2:82083/9‑137 (MQ=255)
      tattatTTCCACCCAAGATAAAGTTTTCTGTCATAAGTTAAACTCCTTCACACCTTGCGAAGTGTAAGCAGACAACTTTTCCATTTGAATCATCGGAGAAACATATTTAATGGTTTGTCGCCTGATAAAGTGCTACAtgtg                                                                                                                                     <  2:134488/141‑1 (MQ=255)
        ttattTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGaa                                                                                                                                   >  2:360168/1‑141 (MQ=255)
          attTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAAtg                                                                                                                                 >  1:298712/1‑141 (MQ=255)
                   cAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTAAATGTGTGTACCt                                                                                                                        >  2:285366/1‑141 (MQ=255)
                    aaaaTAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAATGCTACATGTAAATGTGTGTACCtt                                                                                                                       >  2:208794/1‑141 (MQ=255)
                            gTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTAATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCa                                                                                                               >  1:378102/1‑141 (MQ=255)
                             ttttCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCaa                                                                                                              >  2:43468/1‑141 (MQ=255)
                             ttttCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCaa                                                                                                              >  2:357753/1‑141 (MQ=255)
                             ttttCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCaa                                                                                                              >  2:43481/1‑141 (MQ=255)
                                aacttatgacagaaaaCTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCa                                                                                                           >  1:270337/14‑141 (MQ=255)
                                 cTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCaa                                                                                                          >  1:16642/1‑141 (MQ=255)
                                      aTAAGTTAACCTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCtt                                                                                                     <  1:144238/141‑1 (MQ=255)
                                      aTAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCtt                                                                                                     <  1:381483/141‑1 (MQ=255)
                                           ttAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTAt                                                                                                <  2:298712/141‑1 (MQ=255)
                                            tAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTTTTTTTTTTTTACGGAGAAACATATTTAATGTTTTATCGCGTGATACAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATg                                                                                               >  2:393175/1‑141 (MQ=255)
                                             aaaCTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGc                                                                                              >  2:220844/1‑141 (MQ=255)
                                                   cTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTAGATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGctc                                                                                        <  1:261893/141‑1 (MQ=255)
                                                    ttCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGctct                                                                                       >  2:106610/1‑141 (MQ=255)
                                                       aaaCCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTaaa                                                                                    >  1:178378/1‑141 (MQ=255)
                                                       aaaCCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTCGTCTaaa                                                                                    >  1:158120/1‑141 (MQ=255)
                                                          ccTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTg                                                                                 >  2:50386/1‑141 (MQ=255)
                                                                aaaGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATg                                                                           >  2:326311/1‑141 (MQ=255)
                                                                aaaGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATg                                                                           >  2:307271/1‑141 (MQ=255)
                                                                 aaGTGTCATGAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTACTGTTTTGTCGTGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGa                                                                          <  2:287533/141‑1 (MQ=255)
                                                                 aaGTGTCAGCAGACAACTTTTTCATTTGAATCAGCGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGa                                                                          <  2:232811/141‑1 (MQ=255)
                                                                 aaGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTGTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGa                                                                          <  1:59227/141‑1 (MQ=255)
                                                                 aaGGGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGa                                                                          <  1:299945/141‑1 (MQ=255)
                                                                           agACAACTTTTTTTTTTTTTTTAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCTTTCTCTATATTTGTAATGATAACtttttt                                                                >  2:219844/1‑141 (MQ=255)
                                                                               aaCTTGTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGGGAAAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGttt                                                            <  2:25719/141‑1 (MQ=255)
                                                                                   ttttCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTccc                                                        <  1:245215/141‑1 (MQ=255)
                                                                                         tttGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACttt                                                  <  1:43468/141‑1 (MQ=255)
                                                                                         tttGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACttt                                                  <  1:357753/141‑1 (MQ=255)
                                                                                         tttGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGATACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACttt                                                  <  1:43481/141‑1 (MQ=255)
                                                                                          ttGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAAGCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTc                                                 <  2:161183/141‑1 (MQ=255)
                                                                                           tGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTca                                                >  2:268570/1‑141 (MQ=255)
                                                                                                 aaCGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAAt                                          <  1:360168/141‑1 (MQ=255)
                                                                                                     gagaAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTa                                      <  1:212378/141‑1 (MQ=255)
                                                                                                               tttCATGCTGTGTCGCGTGGTAAATTCCTACATGTGAATGTGTGTACCTTTGTGTCCAATCAAATCTTCGTTTGCCTTCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCTCCCTTTCAACAATTCGGATAGcgtcgtt                            <  2:47313/141‑7 (MQ=255)
                                                                                                               tttAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTc                            <  1:121480/141‑1 (MQ=255)
                                                                                                                            cgcgTGATAAAGTGTTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAAc               >  1:140607/1‑141 (MQ=255)
                                                                                                                              cgTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACAt             >  1:256635/1‑141 (MQ=255)
                                                                                                                               gTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTCTTTTTGGCACCATc            >  2:336300/1‑141 (MQ=255)
                                                                                                                                    aaaGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATcc       <  1:49234/141‑1 (MQ=255)
                                                                                                                                    aaaGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATcc       <  1:220844/141‑1 (MQ=255)
                                                                                                                                     aaGTGCCACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATCCt      <  2:21772/141‑1 (MQ=255)
                                                                                                                                         gCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATCCTTCca  <  1:207563/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
CTAATTTATTATTTCTACCCAAAATAAAGTTTTCTGTCATAAGTTAAACTCCTTCAAACCTTGCAAAGTGTCAGCAGACAACTTTTTCATTTGAATCAACGGAGAAACATATTTAATGTTTTGTCGCGTGATAAAGTGCTACATGTGAATGTGTGTACCTTTATTACCAATCAAAGCTTCGTATGCCTGCTCTAAATTGGTAATGATAACTTTTTTGTTTTCCCCACTTTCAACAAATCGTATCGCAGCTTCTATTTTTGGCAACATCGATCCTTCCA  >  CP000730/2778200‑2778477

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