Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,184,181 (T)5→6 intergenic (+200/‑22) ygiL → / → yi215 putative fimbrial‑like protein/IS2 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,184,1761.T100.0% 113.7 / NA 37intergenic (+195/‑27)ygiL/yi215putative fimbrial‑like protein/IS2 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (20/17);  total (20/17)

AAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACA  >  NC_000913/3184085‑3184201
                                                                                            |                         
aaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                          <  1:3820260/93‑1 (MQ=255)
 aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         >  1:570232/1‑93 (MQ=255)
  atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                        >  2:1600389/1‑93 (MQ=255)
   tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  2:3438089/1‑93 (MQ=255)
   tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  1:4276779/1‑93 (MQ=255)
   tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  2:254515/1‑93 (MQ=255)
   tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       <  2:2503379/93‑1 (MQ=255)
   tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  2:1608367/1‑93 (MQ=255)
    atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGCCACCAAATATACCTTAAATAAGTGATAGTCTTAATACTAGttttt                      >  2:3723154/1‑93 (MQ=18)
    atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                      <  2:34865/93‑1 (MQ=255)
    atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                      <  1:1380605/93‑1 (MQ=255)
     tttattGATCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttt                     >  2:1933791/1‑92 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    >  1:482172/1‑91 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  1:3625368/93‑3 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  1:3625369/93‑3 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  1:3870169/93‑3 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  1:4276701/93‑3 (MQ=255)
      ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  2:3160452/93‑3 (MQ=255)
       tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttag                   >  1:778451/1‑90 (MQ=255)
       tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttag                   <  1:1409194/93‑4 (MQ=255)
       tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttag                   >  1:3787614/1‑90 (MQ=255)
        attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  2:340797/1‑89 (MQ=255)
        attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  1:2453108/93‑5 (MQ=255)
          tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 <  1:3550872/93‑1 (MQ=255)
          tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 >  1:3223475/1‑93 (MQ=255)
          tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 <  1:2247095/93‑1 (MQ=255)
          tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 <  1:1930611/93‑1 (MQ=255)
           gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACtg                >  1:3953361/1‑92 (MQ=255)
           gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                <  2:2050608/93‑1 (MQ=255)
           gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                >  2:1518492/1‑93 (MQ=255)
            ttaaTTATCAGTTATTAGACTGGCCCACTGAATCTACAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               >  2:479388/4‑93 (MQ=14)
             tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGt              >  2:864540/1‑93 (MQ=255)
              cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTc             <  2:1292293/93‑1 (MQ=255)
              cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTc             <  1:4509550/93‑1 (MQ=255)
               aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCa            >  2:1752667/1‑93 (MQ=255)
                ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           >  2:3893076/1‑93 (MQ=255)
                ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           <  1:280899/93‑1 (MQ=255)
                ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           >  1:294879/1‑93 (MQ=255)
                ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           <  2:2125346/93‑1 (MQ=255)
                 tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          >  2:511457/1‑93 (MQ=255)
                 tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  2:2398456/93‑1 (MQ=255)
                 tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  2:512767/93‑1 (MQ=255)
                 tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  2:3398851/93‑1 (MQ=255)
                 tATCCGTTATTAGACTGGCACCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  2:512768/93‑1 (MQ=255)
                  aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTg         <  2:3014926/93‑1 (MQ=255)
                   tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  1:935919/1‑93 (MQ=255)
                   tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        <  1:823415/93‑1 (MQ=255)
                   tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        <  1:823416/93‑1 (MQ=255)
                    ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTGGTCATTGga       <  1:4063417/93‑1 (MQ=255)
                    ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGga       >  1:506007/1‑93 (MQ=255)
                    ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGga       >  2:636706/1‑93 (MQ=255)
                    ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGga       >  1:256504/1‑93 (MQ=255)
                      gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgaga     >  1:2106903/1‑93 (MQ=255)
                       ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    >  1:3522980/1‑93 (MQ=16)
                       ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    <  1:766623/93‑1 (MQ=16)
                       ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    >  2:1886094/1‑93 (MQ=16)
                        tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgagagc   >  1:3683668/1‑91 (MQ=255)
                         attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  2:4733956/1‑93 (MQ=9)
                         attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  2:468237/1‑93 (MQ=9)
                         attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  2:4297008/1‑93 (MQ=9)
                                                                                            |                         
AAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACA  >  NC_000913/3184085‑3184201

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