Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,718,430 C→T G33G (GGG→GGA yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,718,4300CT100.0% 105.8 / NA 30G33G (GGG→GGAyfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/17);  total (13/17)

GGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTT  >  NC_000913/2718352‑2718506
                                                                              |                                                                            
gggCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTtcc                                                                            <  2:485853/81‑1 (MQ=255)
      aCTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTAcg                                                                      <  1:485854/81‑1 (MQ=255)
          tGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTAcgcg                                                                    <  2:485855/79‑1 (MQ=255)
            tttGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGtt                                                                <  2:485856/81‑1 (MQ=255)
              tGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTtg                                                               <  2:485857/80‑1 (MQ=255)
                tCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTgg                                                            >  1:485914/1‑81 (MQ=255)
                  aTCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTggg                                                           <  1:485858/80‑1 (MQ=255)
                  aTCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTggg                                                           >  1:485913/1‑80 (MQ=255)
                    cgcgGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTAc                                                        >  1:485918/1‑81 (MQ=255)
                       gggcggggcgACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGtt                                                     >  2:485904/1‑81 (MQ=255)
                         gcggggcgACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTa                                                   >  1:485907/1‑81 (MQ=255)
                           ggggcgACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTAt                                                  >  2:485921/1‑80 (MQ=255)
                           ggggcgACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATc                                                 <  2:485859/81‑1 (MQ=255)
                              gcgACGGCCACCGTCAGCACGAGATTTTCCTGGGCGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                              <  1:485860/81‑1 (MQ=255)
                                gACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAat                                            >  2:485903/1‑81 (MQ=255)
                                  cGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatca                                          <  1:485861/81‑1 (MQ=255)
                                   ggCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcat                                         >  2:485906/1‑81 (MQ=255)
                                    gCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcatc                                        <  2:485862/81‑1 (MQ=255)
                                      cACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatcatcat                                      >  1:485910/1‑81 (MQ=255)
                                        ccGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAAtcatcatca                                    <  1:485863/81‑1 (MQ=255)
                                          gTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACt                                  >  1:485912/1‑81 (MQ=255)
                                              gcgcGAGATTTTCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCgg                              <  2:485864/78‑1 (MQ=255)
                                                        ttCCTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATc                     <  1:485865/80‑1 (MQ=255)
                                                          ccTGGACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACt                  >  1:485902/1‑81 (MQ=255)
                                                             ggACGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACttt                >  1:485911/1‑80 (MQ=255)
                                                               aCGCCCGCCCCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAc             <  1:485866/81‑1 (MQ=255)
                                                                  cccgcccCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCtt          <  2:485867/81‑1 (MQ=255)
                                                                  cccgcccCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCtt          <  2:485868/81‑1 (MQ=255)
                                                                     gcccCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGcc       >  2:485916/1‑81 (MQ=255)
                                                                     gcccCCTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGcc       <  1:485869/81‑1 (MQ=255)
                                                                          cTTTTCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCttt  <  2:485870/81‑1 (MQ=255)
                                                                              |                                                                            
GGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTT  >  NC_000913/2718352‑2718506

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