Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,988,339 C→T Q40* (CAG→TAG)  yecR → lipoprotein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,3390CT98.2% 186.1 / ‑4.4 55Q40* (CAG→TAG) yecRlipoprotein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/30);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CATTAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACCAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCACGCGGTCCCC  >  NC_000913/1988247‑1988434
                                                                                            |                                                                                               
cATTAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCttt                                                                                          >  2:1219667/1‑100 (MQ=255)
 actaataaCGCTGACGGGGTGTACGGTTACTAGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCtttt                                                                                         >  1:879698/3‑100 (MQ=255)
  ttAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCttttt                                                                                        <  2:1251660/100‑1 (MQ=255)
  ttAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCttttt                                                                                        >  2:489143/1‑100 (MQ=255)
  ttAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCttttt                                                                                        <  1:948298/100‑1 (MQ=255)
   taataaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCtttttt                                                                                       >  2:1231721/1‑100 (MQ=255)
   taataaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCtttttt                                                                                       <  1:243256/100‑1 (MQ=255)
    aataaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTg                                                                                      >  2:597049/1‑100 (MQ=255)
     ataaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgc                                                                                     <  2:932583/100‑1 (MQ=255)
     ataaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgc                                                                                     >  1:1277044/1‑100 (MQ=255)
      taaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgca                                                                                    >  1:842407/1‑100 (MQ=255)
       aaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgcag                                                                                   <  2:582326/100‑1 (MQ=255)
       aaCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgcag                                                                                   <  1:414367/100‑1 (MQ=255)
        aCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgcagc                                                                                  >  1:26992/1‑100 (MQ=255)
        aCGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgcagc                                                                                  <  2:588531/100‑1 (MQ=255)
         cGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTgcagca                                                                                 <  1:450592/100‑1 (MQ=255)
             ggCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGcc                                                                             >  1:825002/1‑100 (MQ=255)
             ggCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGcc                                                                             <  1:565531/100‑1 (MQ=255)
               cGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCa                                                                           >  1:550337/1‑100 (MQ=255)
               cGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCa                                                                           >  2:740832/1‑100 (MQ=255)
               cGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCa                                                                           <  1:385911/100‑1 (MQ=255)
                ggggTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCAt                                                                          <  2:201316/100‑1 (MQ=255)
                ggggTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCAt                                                                          >  1:428294/1‑100 (MQ=255)
                   gtgtACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATAcc                                                                       <  1:876137/100‑1 (MQ=255)
                        cGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCg                                                                  <  1:702200/100‑1 (MQ=255)
                           ttACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTa                                                               <  1:454905/100‑1 (MQ=255)
                           ttACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTa                                                               <  1:268991/100‑1 (MQ=255)
                                cGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCa                                                          >  1:1162468/1‑100 (MQ=255)
                                  ccAAGCCCATGTGAGCGAAGTTGACGCGGCAACATGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAg                                                         <  2:599994/99‑1 (MQ=255)
                                      gCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTcgcg                                                    <  2:845749/100‑1 (MQ=255)
                                       cccATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCgg                                                   >  1:165273/1‑100 (MQ=255)
                                       cccATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCgg                                                   <  2:386979/100‑1 (MQ=255)
                                           tgtgAGCGAAGTTGACGCGGCATCAGGTATTCTACGGCTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATc                                               <  2:879691/100‑1 (MQ=255)
                                                gCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGa                                          >  2:670185/1‑100 (MQ=255)
                                                 cGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGAt                                         >  1:299840/1‑100 (MQ=255)
                                                 cGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTGTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCAGGATCGCCGAt                                         <  2:54727/100‑1 (MQ=255)
                                                    aGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATcgc                                      <  1:740837/100‑1 (MQ=255)
                                                     gTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATcgcg                                     <  1:1152145/100‑1 (MQ=255)
                                                     gTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATcgcg                                     >  1:1062739/1‑100 (MQ=255)
                                                      ttGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCgcgc                                    >  2:658270/1‑100 (MQ=255)
                                                       tGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGCGATCGGCGATCccgcg                                   <  1:1231732/100‑2 (MQ=255)
                                                       tGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACCAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATcgcg                                     >  2:289455/1‑98 (MQ=255)
                                                         aCGCGGCAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATg                                 <  1:662445/100‑1 (MQ=255)
                                                               cAACAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCcagca                           <  1:1219678/100‑1 (MQ=255)
                                                                aaCAGGTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCcagcag                          <  1:40559/100‑1 (MQ=255)
                                                                    ggTATTGTACGGTTAGTTTATGACTAGGCTTTTTTTTTGTATTTTTATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATTCTAGCAGGAAg                      >  2:947020/1‑100 (MQ=255)
                                                                    ggTATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAg                      >  1:298869/1‑100 (MQ=255)
                                                                      tATTGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGc                    >  2:1086840/1‑100 (MQ=255)
                                                                        ttGTACGGTTAGTTTATGACTAGGCTTTTTTGCAGCATGCGCATCCCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCta                  <  1:715322/100‑1 (MQ=255)
                                                                            aCGGTTAGTTTATGACTAGGCTTTTTTGCAGCATACCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGc                    >  2:467591/1‑94 (MQ=255)
                                                                              ggTTAGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATAGAAGTAGCAGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTCTAcaca            >  2:1159438/1‑100 (MQ=255)
                                                                                  aGTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCAcgcg        <  2:842401/100‑1 (MQ=255)
                                                                                   gTTTATGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCACGCgg       >  1:452304/1‑100 (MQ=255)
                                                                                      tATGACTAGGCTTTTTTGGAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCACGCGGTc     <  2:482207/99‑1 (MQ=255)
                                                                                        tGACTAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCACGCGGTcccc  <  1:597053/100‑1 (MQ=255)
                                                                                            |                                                                                               
CATTAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGTATTGTACGGTTAGTTTATGACCAGGCTTTTTTGCAGCATGCCCATACCGATCGTTATGTCAGTCGCGGGATCGCCGATCGCGCATGCCAGCAGGAAGGCTATACCCACGCGGTCCCC  >  NC_000913/1988247‑1988434

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