Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 366,116 A→G F64L (TTT→CTT)  lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,1160AG100.0% 127.4 / NA 36F64L (TTT→CTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/18);  total (18/18)

GGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTT  >  NC_000913/365986‑366252
                                                                                                                                  |                                                                                                                                        
ggATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCaaa                                                                                                                                 >  1:95883/1‑140 (MQ=255)
      gTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGcgcca                                                                                                                           >  1:13035/1‑140 (MQ=255)
      gTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTACCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGcgcca                                                                                                                           <  1:177555/140‑1 (MQ=255)
           gTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGcgccattcgc                                                                                                                      >  2:169575/1‑140 (MQ=255)
                        ggCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcg                                                                                                         <  1:169575/140‑1 (MQ=255)
                         gcgcATCGTAACCGTGCATCTGGCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcgc                                                                                                        <  1:216605/140‑1 (MQ=255)
                         gcgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcgc                                                                                                        <  2:462453/140‑1 (MQ=255)
                         gcgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcgc                                                                                                        <  2:455596/140‑1 (MQ=255)
                          cgcATCGTAACCGTGCATCTGCCAGTTTGAGGTGACGACGACAGTATCGGCCTCAGGAAGATCGCTCTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcgc                                                                                                        <  1:215473/139‑1 (MQ=255)
                          cgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGcgccattc                                                                                                                        >  2:323869/1‑123 (MQ=255)
                          cgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGcgccattc                                                                                                                        <  1:323869/123‑1 (MQ=255)
                          cgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTgcgc                                                                                                        >  2:215473/1‑139 (MQ=255)
                          cgcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCa                                                                                                       >  1:392740/1‑140 (MQ=255)
                           gcATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCaa                                                                                                      <  2:392740/140‑1 (MQ=255)
                            cATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAAc                                                                                                     >  2:427465/1‑140 (MQ=255)
                              tCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTg                                                                                                   >  2:404716/1‑140 (MQ=255)
                               cGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGt                                                                                                  >  2:212180/1‑140 (MQ=255)
                                gTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGtt                                                                                                 <  1:212180/140‑1 (MQ=255)
                                       tGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAgg                                                                                          >  1:235294/1‑140 (MQ=255)
                                            cTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGAt                                                                                     <  2:235294/140‑1 (MQ=255)
                                                   tttGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCg                                                                              <  1:427465/140‑1 (MQ=255)
                                                            acgacgacAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTc                                                                     <  1:38414/140‑1 (MQ=255)
                                                              gacgacAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGc                                                                   >  2:235907/1‑140 (MQ=255)
                                                                      tATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGc                                                           >  2:532216/1‑140 (MQ=255)
                                                                         cGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAg                                                        <  1:404716/140‑1 (MQ=255)
                                                                                      aTCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTgg                                                    >  2:259870/1‑131 (MQ=255)
                                                                                      aTCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTgg                                                    <  1:259870/131‑1 (MQ=255)
                                                                                      aTCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAgggg                                           >  1:152337/1‑140 (MQ=255)
                                                                                       tCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAggggg                                          <  2:95883/140‑1 (MQ=255)
                                                                                       tCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAggggg                                          <  2:152337/140‑1 (MQ=255)
                                                                                          cACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGAtg                                       <  1:532216/140‑1 (MQ=255)
                                                                                               cagccagcTTTCCGGCACCGCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGtgctg                                  >  1:286938/1‑140 (MQ=255)
                                                                                                                  gCTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGt               <  2:286938/140‑1 (MQ=255)
                                                                                                                   cTTCTGGTGCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTa              >  2:199578/1‑140 (MQ=255)
                                                                                                                           gCCGGAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAgg      >  2:415847/1‑140 (MQ=255)
                                                                                                                               gAAGCCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGttt  >  1:381928/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GGATAGGTCACGTTGGTGTAGATGGGCGCATCGTAACCGTGCATCTGCCAGTTTGAGGGGACGACGACAGTATCGGCCTCAGGAAGATCGCACTCCAGCCAGCTTTCCGGCACCGCTTCTGGTGCCGGAAACCAGGCAAAGCGCCATTCGCCATTCAGGCTGCGCAACTGTTGGGAAGGGCGATCGGTGCGGGCCTCTTCGCTATTACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTT  >  NC_000913/365986‑366252

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