Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,461,853:1 +C 8.0% coding (360/936 nt) pyrB ← aspartate carbamoyltransferase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,461,8531.C 8.0% 215.7 / 14.5 50coding (360/936 nt)pyrBaspartate carbamoyltransferase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base . (19/27);  new base C (2/2);  total (21/29)
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.17e-01

CATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAG  >  CP009273/4461714‑4461991
                                                                                                                                            |                                                                                                                                          
cATATGTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCg                                                                                                                                             <  1:839711/140‑1 (MQ=255)
    tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGaca                                                                                                                                        >  2:510565/1‑138 (MQ=255)
    tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGaaa                                                                                                                                        >  1:316164/1‑140 (MQ=255)
    tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGaaa                                                                                                                                        >  2:667754/1‑140 (MQ=255)
    tttCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTGGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGaaa                                                                                                                                        >  1:406399/1‑140 (MQ=255)
       cAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACt                                                                                                                                     <  1:737808/140‑1 (MQ=255)
       cAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACt                                                                                                                                     <  1:685942/140‑1 (MQ=255)
               caaccaTTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGcc                                                                                                                             >  1:303989/1‑140 (MQ=255)
                         gACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGcgccg                                                                                                                   >  2:279420/1‑140 (MQ=255)
                         gACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGcgccg                                                                                                                   >  2:278782/1‑140 (MQ=255)
                                     gTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCg                                                                                                       <  2:729706/140‑1 (MQ=255)
                                                ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATc                                                                                            >  1:590918/1‑140 (MQ=255)
                                                 cTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTATGCGTCGGATGTTGGTTGGAGCCATCACCGGCAGTCAGTACCGGTACATTGCCCGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCa                                                                                           >  1:612121/1‑140 (MQ=255)
                                                         ccTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCa                                                                                   >  1:22148/1‑140 (MQ=255)
                                                               tAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGCGAGCTATTGGTGCCAGTGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGt                                                                            >  2:129526/1‑140 (MQ=255)
                                                               tAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGCGAGCTATTGGTGCCAGTGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGt                                                                            >  2:607258/1‑140 (MQ=255)
                                                                  tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                     <  2:601151/113‑1 (MQ=255)
                                                                  tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGa                                                                                                     >  1:601151/1‑113 (MQ=255)
                                                                       atgTCCAGCACGGTGTGCGGCGGATGTTGGTTGGCGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     <  2:660385/138‑1 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     >  2:508725/1‑140 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     >  2:724131/1‑140 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     <  2:331285/140‑1 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     <  2:819943/140‑1 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     <  2:87002/140‑1 (MQ=255)
                                                                       aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                                                     >  2:91240/1‑140 (MQ=255)
                                                                            cAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATaa                                                                <  1:102151/140‑1 (MQ=255)
                                                                                 aGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGCGAGCTATTGGTGCCAGTGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGaaa                                                          <  1:607258/140‑1 (MQ=255)
                                                                                 aGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGCGAGCTATTGGTGCCAGTGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGaaa                                                          <  1:129526/140‑1 (MQ=255)
                                                                                  ggTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGaaa                                                          <  1:819084/140‑1 (MQ=255)
                                                                                  ggTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGaaa                                                          <  1:333256/140‑1 (MQ=255)
                                                                                       gCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTa                                                     <  1:91240/140‑1 (MQ=255)
                                                                                       gCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTa                                                     <  1:508725/140‑1 (MQ=255)
                                                                                             gATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGcc                                               >  1:510897/1‑140 (MQ=255)
                                                                                              aTGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCa                                              <  2:594650/140‑1 (MQ=255)
                                                                                              aTGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGGGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCa                                              <  2:829084/140‑1 (MQ=255)
                                                                                                            cATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCttt                                <  2:776878/140‑1 (MQ=255)
                                                                                                            cATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCttt                                <  2:979320/140‑1 (MQ=255)
                                                                                                                 ccGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                           <  2:406399/140‑1 (MQ=255)
                                                                                                                 ccGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                           <  1:661059/140‑1 (MQ=255)
                                                                                                                 ccGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                           <  1:486768/140‑1 (MQ=255)
                                                                                                                 ccGGCATTCAGGACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATCACGCTTCACTATCGCATCGACGTAAGTGCTGATAACCGAAAGGGTATCGGCCAGCGTTTCGCCCTTTTTAcc                           <  2:532630/140‑1 (MQ=255)
                                                                                                                 acGGCATTCAGTACCGTTATATTGCCTGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGAGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCTGCCAGCGTTTCGCCCTTTTTAcc                           <  2:316164/139‑1 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCt                                                                                                          <  2:41870/45‑1 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCt                                                                                                          >  1:41870/1‑45 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGTTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGgcg           >  2:613260/1‑140 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAccc                          >  2:517300/1‑125 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAccc                          >  2:21348/1‑125 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAccc                          <  1:517300/125‑1 (MQ=255)
                                                                                                                                 gTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTAccc                          <  1:21348/125‑1 (MQ=255)
                                                                                                                                  tACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGgcgc          <  2:510897/140‑1 (MQ=255)
                                                                                                                                          cGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGgag  <  1:613129/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAG  >  CP009273/4461714‑4461991

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: