Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,003 G→T G131G (GGC→GGA pyrB ← aspartate carbamoyltransferase catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,0030GT97.9% 157.6 / ‑3.9 47G131G (GGC→GGApyrBaspartate carbamoyltransferase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (23/23);  minor base C (1/0);  total (24/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 = 9.56e-01

GAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTC  >  NC_000913/4471873‑4472137
                                                                                                                                  |                                                                                                                                      
gAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCg                                                                                                                               >  2:411647/1‑140 (MQ=255)
      gggAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTc                                                                                                                         <  2:918496/140‑1 (MQ=255)
      gggAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTc                                                                                                                         <  2:564168/140‑1 (MQ=255)
      gggAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTc                                                                                                                         <  2:236788/140‑1 (MQ=255)
        gAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAg                                                                                                                       <  1:1199877/140‑1 (MQ=255)
        gAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAg                                                                                                                       <  1:697476/140‑1 (MQ=255)
         aGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGt                                                                                                                      <  2:240226/140‑1 (MQ=255)
             aaCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCAc                                                                                                                                 <  1:237264/125‑1 (MQ=255)
             aaCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCAc                                                                                                                                 >  2:237264/1‑125 (MQ=255)
             aaCGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCg                                                                                                                  >  1:27915/1‑140 (MQ=255)
                 gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAg                                                                                                                       <  1:1101581/131‑1 (MQ=255)
                 gTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAg                                                                                                                       >  2:1101581/1‑131 (MQ=255)
                   gCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTa                                                                                                               <  1:1124300/137‑1 (MQ=255)
                   gCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTa                                                                                                               >  2:1124300/1‑137 (MQ=255)
                   gCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACAt                                                                                                            <  1:488352/140‑1 (MQ=255)
                               cAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACt                                                                                                >  1:945431/1‑140 (MQ=255)
                                 ggTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCg                                                                                              <  1:248744/140‑1 (MQ=255)
                                    cACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTg                                                                                           >  2:1138450/1‑140 (MQ=255)
                                          ccaTTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAgg                                                                                     <  2:126918/140‑1 (MQ=255)
                                                 gACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGCTGTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGcgccg                                                                              <  2:763219/140‑1 (MQ=255)
                                                     tGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCAcc                                                                          >  1:574306/1‑140 (MQ=255)
                                                            tGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGc                                                                   >  2:187251/1‑140 (MQ=255)
                                                                  gACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                             <  1:1017257/140‑1 (MQ=255)
                                                                  gACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                             >  1:649677/1‑140 (MQ=255)
                                                                       cccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCAt                                                        <  2:611622/140‑1 (MQ=255)
                                                                        ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                             <  2:355758/134‑1 (MQ=255)
                                                                        ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGa                                                             >  1:355758/1‑134 (MQ=255)
                                                                        ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATc                                                       >  1:765332/1‑140 (MQ=255)
                                                                        ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATc                                                       >  1:799058/1‑140 (MQ=255)
                                                                        ccTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGACCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATc                                                       >  2:369434/1‑140 (MQ=255)
                                                                            ggTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTa                                                   <  1:135453/140‑1 (MQ=255)
                                                                                tCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGc                                               <  1:1138450/140‑1 (MQ=255)
                                                                                 ccTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCa                                              >  1:437663/1‑140 (MQ=255)
                                                                                        aGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGt                                       <  2:799058/140‑1 (MQ=255)
                                                                                          tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGTCGTaa                                     >  2:305846/1‑140 (MQ=255)
                                                                                          tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTaa                                     >  1:305056/1‑140 (MQ=255)
                                                                                          tGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTaa                                     >  1:783874/1‑140 (MQ=255)
                                                                                            aaTAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGt                                   >  1:576731/1‑140 (MQ=255)
                                                                                               aaGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCt                                >  1:26492/1‑140 (MQ=255)
                                                                                                aGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTg                               >  1:802508/1‑140 (MQ=255)
                                                                                                 gTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTgg                              >  1:827830/1‑139 (MQ=255)
                                                                                                  tCCAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGAt                             >  1:760400/1‑140 (MQ=255)
                                                                                                    cAGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATaa                           <  2:783874/140‑1 (MQ=255)
                                                                                                     aGCAAGGTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAAc                          <  1:573273/140‑1 (MQ=255)
                                                                                                          ggTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGaaa                     <  2:656045/140‑1 (MQ=255)
                                                                                                           gTTTGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAAt                    >  2:1016950/1‑140 (MQ=255)
                                                                                                            tttGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGAAGTAAGTGCTGATAACCGaa                      <  2:697736/137‑1 (MQ=255)
                                                                                                            tttGCGTCGGATGTTGGTTGGATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGAAGTAAGTGCTGATAACCGaa                      >  1:697736/1‑137 (MQ=255)
                                                                                                                             ttggATCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTc  <  2:802508/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
GAGTCAGGGAGTGAACGGTGCGGCCATATTTCAGGTCACCAACCATTGCGACGTGGAGATTGTCCAGACGCCCCTGGGTTTCCTGAATAGTGAATAAGTCCAGCAAGGTTTGCGTCGGATGTTGGTTGGAGCCATCACCGGCATTCAGTACCGGTACATTGCCGGAAAACTCGGTGGCCAGGCGCGCCGCACCTTCCTGCGGATGACGCATCACTATCGCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTC  >  NC_000913/4471873‑4472137

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