Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 33,974 C→A A1053E (GCG→GAG)  carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi33,9740CA97.0% 96.3 / NA 33A1053E (GCG→GAG) carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (17/15);  minor base G (0/1);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAAGAATGGCGAATATACCTACATCATCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTT  >  NC_000913_3_bme_pgi/33829‑34118
                                                                                                                                                 |                                                                                                                                                
aTCAAGAATGGCGAATATACCTACATCATCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCt                                                                                                                                               <  1:637322/149‑1 (MQ=255)
    aGAATGCCGAATATACCTACATCATCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAAt                                                                                                                                           <  1:249722/149‑1 (MQ=255)
          gCGAATATACCTACATCATCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCACGATTGAGCTGAGTGCGGAt                                                                                                                                     >  1:108477/1‑149 (MQ=255)
              atatACCTACATCATCAAACCCACCTCAGGGCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGa                                                                                                                                 <  2:55884/149‑1 (MQ=255)
                       catcctCAACACCACCTCGGGCNCTCNTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAg                                                                                                                        <  1:634969/149‑1 (MQ=255)
                       catcatCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAg                                                                                                                        >  1:204152/1‑149 (MQ=255)
                           atcaACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAAt                                                                                                                    >  1:376575/1‑149 (MQ=255)
                           atcaACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAAAGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAAt                                                                                                                    >  1:657288/1‑149 (MQ=255)
                            tcaACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAAtt                                                                                                                   >  1:603801/1‑149 (MQ=255)
                             caACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGACTTACGACACCACCCTGAACGGCGGCTTTGCTACTGCGATGGACCTGAATCCCAACGCGACTGAAAAAGTAAttt                                                                                                                  >  2:250134/1‑149 (MQ=255)
                                  ccaccTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTg                                                                                                             >  1:686298/1‑149 (MQ=255)
                                   caccTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGGTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGATCTGAATGCGGATGCGAAAGAAAAAGTAATTTCGGTGc                                                                                                            <  2:561880/149‑1 (MQ=255)
                                      cTCAGGCCGTGGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAgg                                                                                                         >  1:814547/1‑149 (MQ=255)
                                         aGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGGGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGaaa                                                                                                      <  1:57687/149‑1 (MQ=255)
                                                   gCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACGACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGcaca                                                                                            <  2:742045/149‑1 (MQ=255)
                                                                 gcgcgTGATTCGTCGCAGTGCGCTGCAAGATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGc                                                                              <  2:70723/148‑1 (MQ=255)
                                                                 ccGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGc                                                                              <  1:428142/149‑1 (MQ=255)
                                                                   gcgTGATTCGCCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCATGAACGGCGGCTTTGCCACCGCGAGGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCgt                                                                            <  1:99699/149‑1 (MQ=255)
                                                                                        tGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATgca                                                                                                  >  2:750916/1‑106 (MQ=255)
                                                                                        tGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATgca                                                                                                  <  1:750916/106‑1 (MQ=255)
                                                                                        tGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATCTTTTTCAt                                                       >  1:212072/1‑149 (MQ=255)
                                                                                                 aaGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAAtt                                              >  1:287245/1‑149 (MQ=255)
                                                                                                  aGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAAttt                                             >  1:184389/1‑149 (MQ=255)
                                                                                                       aTTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCg                                        >  2:663571/1‑149 (MQ=255)
                                                                                                          aCGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCAGCCGCTAATTTGATCGAAt                                     <  2:603801/149‑1 (MQ=255)
                                                                                                              caccacCCTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACt                                 <  1:780474/149‑1 (MQ=255)
                                                                                                                    ccTGAACGGCGGCTTTGCCACCGCGATGGAGCTGTATGCCCATTCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACg                           >  1:411875/1‑149 (MQ=255)
                                                                                                                    ccTGAACGGCGGCTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACg                           <  2:686298/149‑1 (MQ=255)
                                                                                                                              ggcTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCtgat                 >  1:128966/1‑149 (MQ=255)
                                                                                                                               gcTTTGCCACCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATCGCGACATCCGGCGGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCtgatg                >  1:819922/1‑149 (MQ=255)
                                                                                                                                       aCCGCGATGGAGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt        >  2:528374/1‑149 (MQ=255)
                                                                                                                                           cgATGGAGGTGAATGCCGAAACGGGTGAAAAAGTAATTTCGGTGCATGAAAGGCACTCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGGTAATTAGATCTTAGAGCGAATACGGTTCTCTGATGAGGACCGttttt    <  2:680352/149‑1 (MQ=255)
                                                                                                                                             aTGGGGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt  <  2:575852/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATCAAGAATGGCGAATATACCTACATCATCAACACCACCTCAGGCCGTCGTGCGATTGAAGACTCCCGCGTGATTCGTCGCAGTGCGCTGCAATATAAAGTGCATTACGACACCACCCTGAACGGCGGCTTTGCCACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTT  >  NC_000913_3_bme_pgi/33829‑34118

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: