Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007779 4,371,274 (A)5→3 intergenic (‑9/+296) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,371,2710A.100.0% 185.6 / NA 36intergenic (‑6/+300)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (15/21);  total (15/21)
*NC_0077794,371,2720A.100.0% 185.6 / NA 36intergenic (‑7/+299)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (15/21);  total (15/21)

CGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCC  >  NC_007779/4371172‑4371363
                                                                                                   ||                                                                                           
cGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                                             <  4:34555/101‑1 (MQ=255)
 gATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                                                            >  3:740729/1‑101 (MQ=255)
 gATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                                                            >  1:128848/1‑101 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           >  3:64958/1‑100 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           <  1:1186043/101‑2 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           >  4:144875/1‑100 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           >  1:492190/1‑100 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           <  3:114589/101‑2 (MQ=255)
  aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                           >  1:847060/1‑100 (MQ=255)
     ggTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagc                                                                                        <  1:936350/101‑5 (MQ=255)
        gCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagcggc                                                                                     <  3:1377194/101‑8 (MQ=255)
         cGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                                                  <  1:965370/101‑1 (MQ=255)
           aGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                <  3:666139/101‑1 (MQ=255)
           aGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                <  3:897916/101‑1 (MQ=255)
           aGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                >  2:131821/1‑101 (MQ=255)
           aGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                >  2:132080/1‑101 (MQ=255)
           aGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                <  3:261939/101‑1 (MQ=255)
            ggATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                               <  3:1050612/101‑1 (MQ=255)
                cGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCg                                                                           >  1:554666/1‑101 (MQ=255)
                  tGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGcc                                                                         >  3:711155/1‑101 (MQ=255)
                   gATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCa                                                                        <  1:504607/101‑1 (MQ=255)
                       tATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑TAATAGCGGCGATAACCGCCATGAt                                                                    <  3:1318382/101‑1 (MQ=255)
                       tATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGAt                                                                    >  1:1161767/1‑101 (MQ=255)
                           tCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                                                >  3:1194753/1‑101 (MQ=255)
                           tCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                                                >  3:905941/1‑101 (MQ=255)
                              ggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                                             >  1:1252714/1‑101 (MQ=255)
                                   cGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGAc                                                        <  3:496578/101‑1 (MQ=255)
                                    ggCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACa                                                       <  1:969364/101‑1 (MQ=255)
                                    ggCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACa                                                       <  4:112049/101‑1 (MQ=255)
                                      cgcAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATc                                                     >  3:380968/1‑101 (MQ=255)
                                      cgcAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATc                                                     >  3:1514602/1‑101 (MQ=255)
                                           cagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCTAACg                                                <  4:131959/101‑1 (MQ=255)
                                           cagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACg                                                >  3:955572/1‑101 (MQ=255)
                                           cagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCAGGATAATGGATATGACAGCGAACg                                                <  4:131700/101‑1 (MQ=255)
                                            agcTTTTCTGTCTGCTGACCNCGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACgg                                               >  4:52765/1‑101 (MQ=255)
                                                 ttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt                                          >  3:1469471/1‑101 (MQ=255)
                                                      tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCg                                     <  3:725988/101‑1 (MQ=255)
                                                      tctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCg                                     <  3:571557/101‑1 (MQ=255)
                                                       ctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCgg                                    <  3:532597/101‑1 (MQ=255)
                                                          atgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGtt                                 <  3:855812/101‑1 (MQ=255)
                                                              aCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAAc                             <  3:180696/101‑1 (MQ=255)
                                                              aCCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAAc                             <  1:1445624/101‑1 (MQ=255)
                                                               ccAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACg                            <  3:1283312/101‑1 (MQ=255)
                                                                    tAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAAt                       <  3:116840/101‑1 (MQ=255)
                                                                         cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCgg                  >  3:1570284/1‑101 (MQ=255)
                                                                                 cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAAcc          <  1:5213/101‑1 (MQ=255)
                                                                                         tGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAcccc  >  1:191074/1‑101 (MQ=255)
                                                                                                   ||                                                                                           
CGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCC  >  NC_007779/4371172‑4371363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: