Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 539,975 A→G D258D (GAT→GAC allC ← allantoate deiminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273539,9750AG100.0% 166.8 / NA 51D258D (GAT→GACallCallantoate deiminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/14);  total (37/14)

AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC  >  NZ_CP009273/539890‑540059
                                                                                     |                                                                                    
aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                  >  2:288272/1‑90 (MQ=255)
aatcaatGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGc                                                                                  >  2:665736/1‑90 (MQ=255)
     atGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGTCAAAGGTCAGAACCAGCGGGTCGCCCATc                                                                             >  2:297797/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCtttt                                                                        >  2:410795/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCtttt                                                                        >  1:238810/1‑90 (MQ=255)
          gAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCACCCtttc                                                                        >  1:1387293/1‑89 (MQ=255)
           aaCGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTc                                                                       >  2:411105/1‑90 (MQ=255)
               tGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCtt                                                                   >  1:56626/1‑90 (MQ=255)
                  tcttACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCttttt                                                                <  1:1341205/88‑1 (MQ=255)
                    ttACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcg                                                              >  1:1119717/1‑90 (MQ=255)
                      aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                            >  1:179752/1‑90 (MQ=255)
                      aCCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgac                                                            >  2:555816/1‑90 (MQ=255)
                       cccGGCCCCCCATTTTCCCTTTTCGGGCGGGGGTCCTCCTTGCCAAAGGTCAGAACCAGGGGGTCGGCCATCCTTTTCGCTTTTTcgacc                                                           <  2:1387293/90‑1 (MQ=255)
                         cGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTcgaccga                                                         >  1:1466272/1‑90 (MQ=255)
                           gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                       >  2:186418/1‑90 (MQ=255)
                           gCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACt                                                       >  2:1261748/1‑90 (MQ=255)
                                acaTTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATgg                                                  >  1:568250/1‑90 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCg                                                                      >  2:49255/1‑65 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCg                                                                      <  1:49255/65‑1 (MQ=255)
                                     tACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAAt                                             >  2:900478/1‑90 (MQ=255)
                                        cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          <  1:761770/90‑1 (MQ=255)
                                        cGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          >  2:76638/1‑90 (MQ=255)
                                        cGTATTCGGGCGCAGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCg                                          >  1:540222/1‑90 (MQ=255)
                                           aTTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACt                                       >  2:582046/1‑90 (MQ=255)
                                              cGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaa                                    <  1:158547/90‑1 (MQ=255)
                                               gggCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGaaa                                   <  1:751389/90‑1 (MQ=255)
                                                     ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                             <  2:1119717/90‑1 (MQ=255)
                                                     ggCTCTACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTaa                             >  1:709279/1‑90 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  2:568250/90‑1 (MQ=255)
                                                       ctctACTTTGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  1:1206924/90‑1 (MQ=255)
                                                              ttGCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCa                    >  1:398485/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1530439/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1419838/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1358670/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  >  2:1320332/1‑90 (MQ=255)
                                                                gCCAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacg                  <  1:900478/90‑1 (MQ=255)
                                                                  cAAAGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacgac                >  2:511001/1‑90 (MQ=255)
                                                                   aaaGGTCAGAACCAGCGGGACGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCacgacg               >  2:798674/1‑90 (MQ=255)
                                                                     aGGTCAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGAt             >  2:1155027/1‑90 (MQ=255)
                                                                         cAGAACCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAAcc         >  1:1374158/1‑90 (MQ=255)
                                                                             aCCAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATc     <  1:763141/90‑1 (MQ=255)
                                                                              ccAGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCg    <  2:56626/90‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                     >  1:175029/1‑55 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGa                                     <  2:175029/55‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           <  2:613908/65‑1 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAAc                           >  1:613908/1‑65 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:506326/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:416867/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  2:662696/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  1:699138/1‑90 (MQ=255)
                                                                                aGCGGGTCGCCCATCCTTTTCGCTTTTTCGACCGACGGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGc  >  1:438952/1‑90 (MQ=255)
                                                                                     |                                                                                    
AATCAATGGTGAACGTGGTTTTACCCGGCACCACATTTACCGTATTCGGGCGCGGCTCTACTTTGCCAAAGGTCAGAACCAGCGGATCGCCCATCCTTTTCGCTTTTTCGACCGACTGATGGCAAATGCGACTGAAAGCGTAAACTGTATCACGACGATAACCCATCGGC  >  NZ_CP009273/539890‑540059

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