Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,093,020 T→G 100% F74L (TTT→TTG yqhH → predicted outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,093,0200TG99.2% 468.2 / ‑4.2 122F74L (TTT→TTGyqhHpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (55/66);  total (56/66)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 127 total aligned reads displayed.

ATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCC  >  REL606/3092972‑3093067
                                                |                                               
aTCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGa                                               <  1:4781639/51‑1 (MQ=255)
aTCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGa                                               <  1:214205/51‑1 (MQ=255)
 tCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGAt                                              >  1:8920043/1‑51 (MQ=255)
 tCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAt                                              >  1:2829619/1‑51 (MQ=255)
 tCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAt                                              <  1:3141419/51‑1 (MQ=255)
 tCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAt                                              >  1:9205291/1‑51 (MQ=255)
 tCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAt                                              >  1:5281505/1‑51 (MQ=255)
  ccGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAtt                                             >  1:5147010/1‑51 (MQ=255)
  ccGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAtt                                             <  1:28562/51‑1 (MQ=255)
  ccGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGAtt                                             <  1:7810031/51‑1 (MQ=255)
   cGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTg                                            >  1:10948644/1‑51 (MQ=255)
   cGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTg                                            >  1:9836803/1‑51 (MQ=255)
   cGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTg                                            >  1:8818079/1‑51 (MQ=255)
    gNAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGc                                           >  1:372941/1‑51 (MQ=255)
    gAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGc                                           >  1:9225251/1‑51 (MQ=255)
    gAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGc                                           >  1:1174344/1‑51 (MQ=255)
       gCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTg                                        <  1:4332470/51‑1 (MQ=255)
        cTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgc                                       <  1:5674380/51‑1 (MQ=255)
        cTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgc                                       <  1:2710408/51‑1 (MQ=255)
        cTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgc                                       <  1:9656834/51‑1 (MQ=255)
        cTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgc                                       <  1:5680638/51‑1 (MQ=255)
         tAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgcg                                      <  1:5132348/51‑1 (MQ=255)
         tAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgcg                                      >  1:5191330/1‑51 (MQ=255)
          aaCAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgcgc                                     >  1:3764373/1‑51 (MQ=255)
          aaCAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTgcgc                                     >  1:4268431/1‑51 (MQ=255)
           aCAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgct                                    <  1:4563704/51‑1 (MQ=255)
            cAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctg                                   <  1:2087464/51‑1 (MQ=255)
            cAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctg                                   >  1:7181799/1‑51 (MQ=255)
             agagCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctgc                                  >  1:3302440/1‑51 (MQ=255)
             agagCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctgc                                  >  1:1360910/1‑51 (MQ=255)
               agCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCtt                                <  1:2980476/51‑1 (MQ=255)
               agCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCtt                                <  1:6043871/51‑1 (MQ=255)
               agCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCtt                                <  1:8441869/51‑1 (MQ=255)
               agCCAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCtt                                <  1:1612729/51‑1 (MQ=255)
                 ccAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctg                                   <  1:3542305/46‑1 (MQ=255)
                  cAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCgctgc                                  <  1:984474/46‑1 (MQ=255)
                  cAATACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCg                             >  1:8404151/1‑51 (MQ=255)
                   aaTACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGt                            <  1:6539575/51‑1 (MQ=255)
                   aaTACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGt                            <  1:8520411/51‑1 (MQ=255)
                    aTACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTa                           >  1:7589747/1‑51 (MQ=255)
                    aTACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTa                           <  1:1226080/51‑1 (MQ=255)
                     tACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTAt                          >  1:2245228/1‑51 (MQ=255)
                     tACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTAt                          >  1:10604528/1‑51 (MQ=255)
                     tACGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTAt                          >  1:6018684/1‑51 (MQ=255)
                      aCGCGTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATg                         >  1:1171976/1‑51 (MQ=255)
                       cgcgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGt                        <  1:242897/51‑1 (MQ=255)
                       cgcgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGt                        <  1:1387442/51‑1 (MQ=255)
                        gcgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTa                       <  1:6201018/51‑1 (MQ=255)
                        gcgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTa                       <  1:8275325/51‑1 (MQ=255)
                         cgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTAc                      >  1:4958781/1‑51 (MQ=255)
                         cgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTAc                      >  1:8317242/1‑51 (MQ=255)
                         cgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTAc                      <  1:7693265/51‑1 (MQ=255)
                         cgTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTAc                      <  1:8622572/51‑1 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     <  1:7579584/51‑1 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     <  1:3405725/51‑1 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     <  1:2541426/51‑1 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     <  1:1688276/51‑1 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     >  1:2287318/1‑51 (MQ=255)
                          gTCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACg                     >  1:378898/1‑51 (MQ=255)
                           tCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGc                    <  1:6414082/51‑1 (MQ=255)
                           tCTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGc                    >  1:1424741/1‑51 (MQ=255)
                            cTTGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCa                   <  1:9462755/51‑1 (MQ=255)
                             ttGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAg                  <  1:5812728/51‑1 (MQ=255)
                             ttGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAg                  <  1:4285717/51‑1 (MQ=255)
                             ttGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAg                  >  1:9272077/1‑51 (MQ=255)
                             ttGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAg                  <  1:4871262/51‑1 (MQ=255)
                              tGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGa                 <  1:7818832/51‑1 (MQ=255)
                              tGATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGa                 >  1:5301112/1‑51 (MQ=255)
                               gATGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGaa                <  1:6225784/51‑1 (MQ=255)
                                aTGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAAt               >  1:9307399/1‑51 (MQ=255)
                                 tGCTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATg              >  1:9721803/1‑51 (MQ=255)
                                   cTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGAt            >  1:1349972/1‑51 (MQ=255)
                                   cTCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGAt            >  1:6473406/1‑51 (MQ=255)
                                    tCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATa           <  1:7818285/51‑1 (MQ=255)
                                    tCAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATa           <  1:6992329/51‑1 (MQ=255)
                                     cAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaa          <  1:5560756/51‑1 (MQ=255)
                                     cAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaa          <  1:5927145/51‑1 (MQ=255)
                                     cAGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaa          <  1:10191379/51‑1 (MQ=255)
                                      aGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaa         <  1:2429255/51‑1 (MQ=255)
                                      aGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaa         >  1:3258130/1‑51 (MQ=255)
                                      aGGACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaa         >  1:1259762/1‑51 (MQ=255)
                                       ggACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaa        <  1:6816402/51‑1 (MQ=255)
                                       ggACTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaa        <  1:5747644/51‑1 (MQ=255)
                                         aCTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaa      <  1:5549920/51‑1 (MQ=255)
                                         aCTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaa      <  1:1912600/51‑1 (MQ=255)
                                         aCTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaa      <  1:8521265/51‑1 (MQ=255)
                                         aCTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaa      <  1:10640272/51‑1 (MQ=255)
                                          cTATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaaa     >  1:4259799/1‑51 (MQ=255)
                                           tATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAAt    <  1:1844397/51‑1 (MQ=255)
                                           tATTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAAt    >  1:3574724/1‑51 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   <  1:828430/51‑1 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   <  1:9768239/51‑1 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   <  1:8992500/51‑1 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   >  1:3510995/1‑51 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   <  1:9725856/51‑1 (MQ=255)
                                            aTTTGGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATc   >  1:7259781/1‑51 (MQ=255)
                                             tttgGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATcc  <  1:2255949/51‑1 (MQ=255)
                                             tttgGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATcc  <  1:1290852/51‑1 (MQ=255)
                                             tttgGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATcc  >  1:1405410/1‑51 (MQ=255)
                                             tttgGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATcc  >  1:10188490/1‑51 (MQ=255)
                                                |                                               
ATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCC  >  REL606/3092972‑3093067

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