Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,860,370 T→C E28E (GAA→GAG yeaA ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,860,3700TC100.0% 166.0 / NA 63E28E (GAA→GAGyeaAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/38);  total (25/38)

ATATACGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGTTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAG  >  NC_000913/1860322‑1860440
                                                |                                                                      
cagtaCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTaa                                             <  2:2334189/73‑1 (MQ=255)
cagtaCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTaa                                             <  1:4656919/73‑1 (MQ=255)
 agtaCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTaaa                                            >  2:931482/3‑76 (MQ=255)
  gtaCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGCCTCGGTCCCATGAGTCTGCGTCACGTaaaa                                           >  2:4538582/2‑76 (MQ=255)
    aCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACt                                         <  1:1104533/76‑1 (MQ=255)
    aCGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACt                                         <  1:2195760/76‑1 (MQ=255)
     cGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTg                                        <  2:468947/76‑1 (MQ=255)
      gCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGc                                       >  2:2404462/1‑76 (MQ=255)
      gCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGc                                       <  1:2660393/76‑1 (MQ=255)
        cGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAt                                     >  2:394531/1‑76 (MQ=255)
         gTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtc                                    <  2:1019267/76‑1 (MQ=255)
           cACTCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtctc                                  <  1:634617/76‑1 (MQ=255)
           cACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtctc                                  <  1:2372430/76‑1 (MQ=255)
           cACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtctc                                  <  2:541991/76‑1 (MQ=255)
           cACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtctc                                  <  1:2812273/76‑1 (MQ=255)
           cACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCAtctc                                  <  2:1046976/76‑1 (MQ=255)
            aCGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCg                                 <  2:2111547/76‑1 (MQ=255)
            aCGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCg                                 >  1:4734130/1‑76 (MQ=255)
            aCGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCg                                 <  2:433361/76‑1 (MQ=255)
            aCGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCg                                 <  2:2619356/76‑1 (MQ=255)
             cGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCgg                                <  2:2604888/76‑1 (MQ=255)
              gCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTACAACTGCATCTCggc                               >  2:2223447/1‑75 (MQ=255)
              gCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGa                               >  2:2215273/1‑76 (MQ=255)
              gCTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGa                               >  1:1793529/1‑76 (MQ=255)
              gCTTGTTAGGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGa                               <  2:146610/76‑1 (MQ=255)
               cTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCggcc                              >  2:4606598/1‑74 (MQ=255)
               cTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGAc                              >  2:4544902/1‑76 (MQ=255)
               cTTGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGAc                              >  1:133864/1‑76 (MQ=255)
                ttGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACa                             <  1:765392/76‑1 (MQ=255)
                ttGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACa                             >  1:2826965/1‑76 (MQ=255)
                ttGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACa                             >  2:787423/1‑76 (MQ=255)
                 tGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACaa                            <  1:898018/76‑1 (MQ=255)
                 tGTTATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACaa                            >  2:4761912/1‑76 (MQ=255)
                    tATGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAAtt                         >  1:3227864/1‑76 (MQ=255)
                     aTGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAAttt                        <  2:2582429/76‑1 (MQ=255)
                     aTGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCGCGGACAAAttt                        <  2:2608900/76‑1 (MQ=255)
                      tGCAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAAtttt                       <  2:3308638/76‑1 (MQ=255)
                        cAGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAAtttttt                     <  1:3862545/76‑1 (MQ=255)
                         aGTAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAAttttttt                    <  1:3408015/76‑1 (MQ=255)
                           tAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTc                  >  1:3656575/1‑76 (MQ=255)
                           tAAACGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTc                  <  1:2419779/76‑1 (MQ=255)
                             aaCGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAg                <  1:1544617/76‑1 (MQ=255)
                               cGACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGtt              <  2:4787854/76‑1 (MQ=255)
                                gACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGttt             >  1:103528/1‑75 (MQ=255)
                                gACCCGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGttc             >  1:1241554/1‑76 (MQ=255)
                                 gcccGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAACTTTTTTTTCAGttct            <  1:3272667/75‑1 (MQ=255)
                                  cccGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGttctt           <  1:3064767/76‑1 (MQ=255)
                                   ccGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGttcttc          <  2:3043214/76‑1 (MQ=255)
                                    cGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTtcttct         <  1:2577422/76‑1 (MQ=255)
                                    cGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTtcttct         >  2:392639/1‑76 (MQ=255)
                                    cGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTtcttct         <  2:1735202/76‑1 (MQ=255)
                                    cGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTtcttct         <  1:4846968/76‑1 (MQ=255)
                                    cGTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTtcttct         <  2:4595908/76‑1 (MQ=255)
                                     ggaaaTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTg        <  2:240496/74‑1 (MQ=255)
                                     gTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTCCAACTGCATCTCGGACAAATTTTTTTTCAGTTTTTCTg        >  2:4328722/1‑76 (MQ=255)
                                     gTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTg        >  2:4282109/1‑76 (MQ=255)
                                     gTAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTg        >  1:2370583/1‑76 (MQ=255)
                                      tAAATGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGc       <  1:55868/76‑1 (MQ=255)
                                       aaaTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGcc      >  2:2756910/1‑76 (MQ=255)
                                       aaaTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGcc      >  2:2728582/1‑76 (MQ=255)
                                         aTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGa    >  2:1403881/1‑76 (MQ=255)
                                         aTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGa    <  2:4068875/76‑1 (MQ=255)
                                         aTGGCGGCTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGa    <  1:1245743/76‑1 (MQ=255)
                                          gggcggcTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCGCGGACAAATTTTTTTTCAGTTCTTCTGCCGaa   <  2:195356/75‑1 (MQ=255)
                                           ggcggcTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGaac  >  2:199771/1‑75 (MQ=255)
                                                |                                                                      
ATATACGCCGTCACGCTTGTTATGCAGTAAACGACCCGTAAATGGCGGTTCTGTCCCATGATTCTGCGTCACGTAAAACTGCATCTCGGACAAATTTTTTTTCAGTTCTTCTGCCGAAG  >  NC_000913/1860322‑1860440

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