Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,093,686 T→C 100% V130A (GTA→GCA)  ycdT → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,093,6860TC100.0% 271.3 / NA 77V130A (GTA→GCA) ycdTpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/31);  total (46/31)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTG  >  AP009048/1093619‑1093752
                                                                   |                                                                  
gtCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAtt                                                                 >  1:4109603/1‑71 (MQ=255)
   cAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                              >  1:550927/1‑71 (MQ=255)
    aaaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGt                                                              >  1:207670/1‑70 (MQ=255)
     aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                             >  1:4212671/1‑70 (MQ=255)
     aaaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTg                                                             >  1:4814382/1‑70 (MQ=255)
      aaTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                           >  1:2072727/1‑71 (MQ=255)
       aTGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCAt                                                                  <  1:2484347/63‑1 (MQ=255)
          aCATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaa                                                        <  1:3436679/70‑1 (MQ=255)
           cATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                      >  1:2163140/1‑71 (MQ=255)
           cATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                      <  1:1557459/71‑1 (MQ=255)
            aTTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATa                                                         <  1:1187973/67‑1 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                      >  1:2918751/1‑68 (MQ=255)
              tGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGt                                                     >  1:2168602/1‑69 (MQ=255)
                cAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                 <  1:4379949/71‑1 (MQ=255)
                 aaTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAg                                                      >  1:2495695/1‑65 (MQ=255)
                 aaTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGt                                                 >  1:1319076/1‑70 (MQ=255)
                 aaTCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTg                                                <  1:3823977/71‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGa                                                            <  1:1789639/57‑1 (MQ=255)
                   tCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGa                                               >  1:3351010/1‑70 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                             <  1:1526261/68‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                             <  1:2645627/68‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:5165647/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:3644436/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:4433199/67‑1 (MQ=255)
                       ttattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGATAGTGaa                                              <  1:4700445/67‑1 (MQ=255)
                        tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                             <  1:1383050/67‑1 (MQ=255)
                        tattTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                          >  1:5040219/1‑70 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                            >  1:5017723/1‑67 (MQ=255)
                         attTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAaca                                          >  1:3032077/1‑69 (MQ=255)
                           ttGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:430497/63‑1 (MQ=255)
                             gTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaa                                             >  1:954102/1‑62 (MQ=255)
                              tttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                   <  1:808128/71‑1 (MQ=255)
                               ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATaaa                                                       >  1:3348679/1‑50 (MQ=255)
                               ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAacac                                         >  1:3917418/1‑64 (MQ=255)
                               ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAAcaca                                        >  1:946997/1‑65 (MQ=255)
                               ttCGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                     >  1:1253466/1‑68 (MQ=255)
                                 cGTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                    <  1:1410476/51‑1 (MQ=255)
                                  gTCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                    >  1:4911814/1‑66 (MQ=255)
                                   tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:3848374/55‑1 (MQ=255)
                                   tCAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCa                              <  1:296648/71‑1 (MQ=255)
                                    cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGaa                                                           >  1:997009/1‑41 (MQ=255)
                                    cAGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAAc                                    <  1:3212457/64‑1 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                    >  1:1288027/1‑47 (MQ=255)
                                     aGATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGtt                                                    >  1:4800258/1‑47 (MQ=255)
                                       aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatt                           >  1:5199593/1‑70 (MQ=255)
                                        tGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                         >  1:3398243/1‑71 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                        <  1:4663161/70‑1 (MQ=255)
                                          aGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                        <  1:2663816/70‑1 (MQ=255)
                                           gTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttattt                      <  1:2606096/71‑1 (MQ=255)
                                            tttGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                     >  1:211701/1‑71 (MQ=255)
                                              tgtgCATATTAATATTTCTGGCATTGGTGAATAAAGTGAGTGAAAACACAAAACAGCGCAat                            >  1:3556327/1‑62 (MQ=255)
                                              tgtgCATATTAATATTACTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                    >  1:5175995/1‑70 (MQ=255)
                                                tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                  <  1:944927/70‑1 (MQ=255)
                                                tgCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCa                  >  1:663936/1‑70 (MQ=255)
                                                 gCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAgc                                 >  1:1404745/1‑54 (MQ=255)
                                                   atatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGaa                                              <  1:4737986/39‑1 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaa             <  1:4827688/71‑1 (MQ=255)
                                                    tatTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaa              >  1:1279066/1‑70 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAattta                         <  1:3014826/57‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCaaaaaaa            <  1:3723756/70‑1 (MQ=255)
                                                      ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAAt           >  1:4366049/1‑71 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGaaaa                                            >  1:3134821/1‑37 (MQ=255)
                                                       tAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACaaaa                                     >  1:2336974/1‑44 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACa                                   >  1:2483560/1‑45 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                    >  1:804064/1‑60 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa         >  1:2227304/1‑71 (MQ=255)
                                                        aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGa         >  1:4387792/1‑71 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                        >  1:915130/1‑54 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttat                        >  1:1071080/1‑54 (MQ=255)
                                                          tatTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTAtttt                     >  1:5009292/1‑57 (MQ=255)
                                                           atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAatttatt                       >  1:3640593/1‑54 (MQ=255)
                                                           atTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACtt      >  1:2198743/1‑71 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt     <  1:1857936/71‑1 (MQ=255)
                                                            tttCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACttt     >  1:5041190/1‑71 (MQ=255)
                                                               cTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgt   <  1:812121/70‑1 (MQ=255)
                                                                tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTc                    <  1:3867220/52‑1 (MQ=255)
                                                                tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTtgtg  >  1:544916/1‑70 (MQ=255)
                                                                   |                                                                  
GTCCAAAATGACATTGCAATCTATTATTTGTTTCGTCAGATGAGTTTGTGCATATTAATATTTCTGGTATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTG  >  AP009048/1093619‑1093752

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