Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,986,136 T→C intergenic (‑8/‑789) araF ← / → ftnB arabinose ABC transporter periplasmic binding protein/putative ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,986,1360TC100.0% 77.1 / NA 22intergenic (‑8/‑789)araF/ftnBarabinose ABC transporter periplasmic binding protein/putative ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/16);  total (6/16)

CGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGAAAACAAAGCATTACC  >  NC_000913/1986060‑1986211
                                                                            |                                                                           
cGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCa                                                                          <  1:339057/80‑1 (MQ=255)
cGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCa                                                                          <  1:339056/80‑1 (MQ=255)
 gCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGc                                                                        <  1:339058/81‑1 (MQ=255)
 gCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGc                                                                        >  2:339087/1‑81 (MQ=255)
 gCCAAAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGc                                                                        >  1:339083/1‑81 (MQ=255)
   cATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCtt                                                                      <  1:339059/81‑1 (MQ=255)
    aTAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCttt                                                                     <  2:339060/81‑1 (MQ=255)
      aGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAg                                                                   >  2:339090/1‑81 (MQ=255)
               tgACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGtttt                                                          <  1:339062/81‑1 (MQ=255)
               tgACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGttt                                                           <  2:339061/80‑1 (MQ=255)
                gACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTtg                                                         >  1:339085/1‑81 (MQ=255)
                 aCATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTtgt                                                        <  2:339063/81‑1 (MQ=255)
                    tAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTa                                                     <  1:339064/81‑1 (MQ=255)
                       cGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAggg                                                  <  1:339065/81‑1 (MQ=255)
                              cAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAAc                                           <  1:339066/81‑1 (MQ=255)
                                             cAGGGCTTTAGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGtt                             >  1:339091/1‑80 (MQ=255)
                                                      aGTAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATc                    <  2:339067/80‑1 (MQ=255)
                                                        tAAATTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGa                 <  1:339068/81‑1 (MQ=255)
                                                           aTTTGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGaaaa              <  2:339069/81‑1 (MQ=255)
                                                            tttGTGCATGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGAAAAc             <  1:339070/81‑1 (MQ=255)
                                                                   aTGGTTCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGAAAACAAAGCAt      >  1:339088/1‑81 (MQ=255)
                                                                       ttCTCCCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGAAAACAAAGCATTAcc  <  1:339071/81‑1 (MQ=255)
                                                                            |                                                                           
CGCCATAGCGGATTGTGACATAACGGCTGCCAGACCAATGGCTGCCAGGGCTTTAGTAAATTTGTGCATGGTTCTCTCCAGCTTTAGTGTCGTTTTGTGTAGGGCAAAAACGAATGACATTCGTTAAATTAATCGGAAAACAAAGCATTACC  >  NC_000913/1986060‑1986211

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