Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,538,785 C→T *239* (TAG→TAA)  nanC ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,7850CT100.0% 176.4 / NA 55*239* (TAG→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/32);  total (23/32)

TCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACA  >  NC_000913/4538650‑4538920
                                                                                                                                       |                                                                                                                                       
tCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCtattac                                                                                                                                      >  1:1155201/1‑139 (MQ=255)
   ttATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGt                                                                                                                                   <  1:1248998/139‑1 (MQ=255)
      tCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                >  1:381822/1‑139 (MQ=255)
        aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                              <  1:717104/139‑1 (MQ=255)
          ttGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                       >  2:640401/1‑128 (MQ=255)
          ttGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATg                                                                                                                            >  1:241499/1‑139 (MQ=255)
             cTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         >  2:703702/1‑139 (MQ=255)
                   tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   <  2:241499/139‑1 (MQ=255)
                   tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   <  1:703702/139‑1 (MQ=255)
                           ttAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATaa                                                                                                           <  2:282021/139‑1 (MQ=255)
                                        gAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         >  1:1056945/1‑112 (MQ=255)
                                        gAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         <  2:1056945/112‑1 (MQ=255)
                                          gTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATaaa                                                                                            <  2:1294686/139‑1 (MQ=255)
                                             tCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTa                                                                                         >  2:261630/1‑139 (MQ=255)
                                                   aaTACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                   <  2:123004/139‑1 (MQ=255)
                                                     tACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                                 <  1:526939/139‑1 (MQ=255)
                                                            gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                                           <  2:119446/90‑1 (MQ=255)
                                                            gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGACTATTACAGTTTAAATGa                                                                                                                           >  1:119446/1‑90 (MQ=255)
                                                                                  tcatGATAGCAAGCGCCGTTATTGTTTTCTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               >  2:847883/1‑112 (MQ=255)
                                                                                  tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               <  1:847883/112‑1 (MQ=255)
                                                                                  tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGt                                                    <  1:517949/139‑1 (MQ=255)
                                                                                   catGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  2:381822/139‑1 (MQ=255)
                                                                                   catGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  2:46893/139‑1 (MQ=255)
                                                                                         aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               >  2:1012688/1‑105 (MQ=255)
                                                                                         aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               <  1:1012688/105‑1 (MQ=255)
                                                                                            aaGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGAt                                                                                                             <  2:313478/72‑1 (MQ=255)
                                                                                            aaGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGAt                                                                                                             >  1:313478/1‑72 (MQ=255)
                                                                                            aaGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtat                                          >  2:5087/1‑139 (MQ=255)
                                                                                             aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatg                                         <  2:1155201/139‑1 (MQ=255)
                                                                                              gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   >  1:649425/1‑128 (MQ=255)
                                                                                              gcgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  2:649425/128‑1 (MQ=255)
                                                                                               cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       >  2:350341/1‑139 (MQ=255)
                                                                                               cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       <  2:313060/139‑1 (MQ=255)
                                                                                               cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       >  1:585371/1‑139 (MQ=255)
                                                                                                gcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  2:457390/126‑1 (MQ=255)
                                                                                                gcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   >  1:457390/1‑126 (MQ=255)
                                                                                                   gTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGc                                   >  2:1087871/1‑139 (MQ=255)
                                                                                                              aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                        <  1:5087/139‑1 (MQ=255)
                                                                                                              aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                        <  2:497378/139‑1 (MQ=255)
                                                                                                              aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTG‑CGTTTTAAAAACCCTGACGGTCAAGGTAGGCATATTCTTTGTATGGCGTTATTTTGGGc                       >  2:1089794/1‑139 (MQ=255)
                                                                                                               ttCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGc                       >  2:90672/1‑139 (MQ=255)
                                                                                                                tCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                      <  1:863263/139‑1 (MQ=255)
                                                                                                                tCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCg                      <  1:534984/139‑1 (MQ=255)
                                                                                                                    tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtt                                                                                                      >  2:485859/1‑55 (MQ=255)
                                                                                                                    tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtt                                                                                                      <  1:485859/55‑1 (MQ=255)
                                                                                                                    tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                                 <  2:508194/108‑1 (MQ=255)
                                                                                                                    tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                                 >  1:508194/1‑108 (MQ=255)
                                                                                                                    tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCat                  >  2:754175/1‑139 (MQ=255)
                                                                                                                     gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCata                 <  1:1087871/139‑1 (MQ=255)
                                                                                                                        aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATgg              <  1:1323296/139‑1 (MQ=255)
                                                                                                                        aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATgg              <  1:90672/139‑1 (MQ=255)
                                                                                                                         cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAAt                                                                                                                  >  1:133327/1‑38 (MQ=255)
                                                                                                                         cTGTCTCCTGTCTATTACAGTTTAAATGACACACCAAt                                                                                                                  <  2:133327/38‑1 (MQ=255)
                                                                                                                            tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTAtt          >  2:752455/1‑139 (MQ=255)
                                                                                                                                cTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGt        >  1:749008/1‑137 (MQ=255)
                                                                                                                                cTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGt        <  2:749008/137‑1 (MQ=255)
                                                                                                                                cTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAg      <  1:1229562/139‑1 (MQ=255)
                                                                                                                                    ctattaCAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACa  <  1:261630/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TCGTTATCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACA  >  NC_000913/4538650‑4538920

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