Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,519,066 T→G N192H (AAT→CAT)  yjhA ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,519,0660TG90.7% 160.5 / ‑0.5 44N192H (AAT→CAT) yjhAN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base G (26/13);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATG  >  REL606/4518993‑4519137
                                                                         |                                                                       
tGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCAtgtt                                                                       <  1:683758/76‑1 (MQ=255)
tGACGGTCAAGGTAGTCATATTCTATGTATGGCGNAATATCGGGCGTCATANGACATTGTAGAACAAANGCAtgtt                                                                       >  2:953948/1‑76 (MQ=255)
  aCGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTcc                                                                     <  1:298262/76‑1 (MQ=255)
      tCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCg                                                                 >  1:714366/1‑76 (MQ=255)
       cAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGc                                                                >  1:377952/1‑76 (MQ=255)
        aaGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGcc                                                               >  1:807970/1‑76 (MQ=255)
         aGGTAGTCATATACTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGccc                                                              >  1:485787/1‑76 (MQ=255)
                cATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCtt                                                       >  1:842380/1‑76 (MQ=255)
                 atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTa                                                      >  2:64475/1‑76 (MQ=255)
                 atatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTa                                                      >  1:195669/1‑76 (MQ=255)
                  tatTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCCTTTCTTAt                                                     >  1:30707/1‑76 (MQ=255)
                   atTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATg                                                    >  2:121968/1‑76 (MQ=255)
                   atTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATg                                                    <  1:395831/76‑1 (MQ=255)
                   atTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATg                                                    <  1:395834/76‑1 (MQ=255)
                     tCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGt                                                  >  2:528808/1‑76 (MQ=255)
                     tCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGt                                                  >  2:454015/1‑76 (MQ=255)
                       tatgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGGCCATTTCTTATGGttt                                                <  2:1075735/76‑1 (MQ=255)
                       tatgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGttt                                                >  1:302414/1‑76 (MQ=255)
                       tatgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGttt                                                >  2:649232/1‑76 (MQ=255)
                        atgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTg                                               <  2:88731/76‑1 (MQ=255)
                        atgtatgACGTAATATCGGGCGTCATATGGTATTGTAGAACGAATGCNTGTTCCNTCGCCCATTTCTTATGCtttn                                               >  1:905437/1‑75 (MQ=255)
                         tgtatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGc                                              >  1:433926/1‑76 (MQ=255)
                           tatgGCGTAATATCGGTCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCAt                                            <  2:727243/76‑1 (MQ=255)
                           tatgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCAt                                            >  1:623374/1‑76 (MQ=255)
                             tgGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAg                                          <  1:1041000/76‑1 (MQ=255)
                              ggCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGCTTTGCATAGc                                         <  1:146192/76‑1 (MQ=255)
                               gcgcAATATCGGGCGTCATATGGTAGTGTAGAACAAATGCATCTTCCGTCGCCCATTTCTTATGGTTTGCATAGCg                                        <  2:924500/76‑1 (MQ=255)
                                 gTAATATCGGGCGTCATATGGTANTCTAGAANANATGCATCTTCCATCGCCCATTTCTGAAGGTTTGCATAGCGAt                                      >  1:997847/1‑76 (MQ=255)
                                  tAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATa                                     >  2:797532/1‑76 (MQ=255)
                                  tAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATa                                     <  1:643885/76‑1 (MQ=255)
                                  tAATATCGGGCGTCATATGGTATTGTAGAACAAAGGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATa                                     <  1:643881/76‑1 (MQ=255)
                                   aaTATCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATaa                                    <  2:151560/76‑1 (MQ=255)
                                    atatCGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAAt                                   >  2:801809/1‑76 (MQ=255)
                                        cGGGCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGtt                               >  1:525434/1‑76 (MQ=255)
                                         gggCGTCATATGGTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGtta                              >  2:1034466/1‑75 (MQ=255)
                                              tCATATGGTATTATAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGttt                         >  2:158220/1‑76 (MQ=255)
                                               cATATGGTATTGTAGAACTAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTg                        <  1:827786/76‑1 (MQ=255)
                                               cATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGCTTTGCATAGCGATAATCGTTCTGTTTg                        <  1:106364/76‑1 (MQ=255)
                                               cATATGGTATTGCAGAACCAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTg                        <  1:827773/76‑1 (MQ=255)
                                                     gTATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGtat                  <  2:673412/76‑1 (MQ=255)
                                                      tATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCTCTCTGTTTGCTGtata                 >  1:204808/1‑76 (MQ=255)
                                                      tATTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGtata                 <  1:125984/76‑1 (MQ=255)
                                                       aTTGTAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTGTAg                >  1:156175/1‑76 (MQ=255)
                                                          gTAGAACAAATGCATGTNCCGTCGCCCATTTCTTATGGTTTGCATAGCGATNATCGATGTGAGAGCTGTATAGcgt             >  2:893073/1‑76 (MQ=255)
                                                           tAGAACAAATGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTTTAGcgtc            <  1:689844/76‑1 (MQ=255)
                                                                cAAATGCATTTTCCGACGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGttca       <  2:139069/76‑3 (MQ=255)
                                                                  aaTGCATGTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGcc     >  1:206569/1‑76 (MQ=255)
                                                                     gCATGTTCCGTCGCCCATTTCTTATGGTTTGCATACCGATAATCGTTCTGTTTGCAGTATAGCGTCGTTTGCCATg  >  1:811587/1‑76 (MQ=255)
                                                                         |                                                                       
TGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACAAATGCATTTTCCGTCGCCCATTTCTTATGGTTTGCATAGCGATAATCGTTCTGTTTGCTGTATAGCGTCGTTTGCCATG  >  REL606/4518993‑4519137

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