Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,103,105 C→A 5.8% pseudogene (285/446 nt) wbbL ← pseudogene, lipopolysaccharide biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,103,1050CA 5.8% 244.4 / 5.3 69pseudogene (285/446 nt)wbbLpseudogene, lipopolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGG  >  NC_000913/2102965‑2103252
                                                                                                                                            |                                                                                                                                                   
acaATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                             >  3:145851/1‑149 (MQ=255)
         tCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                      >  4:430320/1‑147 (MQ=255)
         tCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                      <  3:430320/147‑1 (MQ=255)
             aaaGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTc                                                                                                                                >  3:144947/1‑149 (MQ=255)
             aaaGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTc                                                                                                                                <  4:144947/149‑1 (MQ=255)
                    aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                              >  2:203602/1‑128 (MQ=255)
                    aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                              <  1:203602/128‑1 (MQ=255)
                    aGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACa                                                                                                                         >  3:377405/1‑149 (MQ=255)
                     ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                              <  1:941144/127‑1 (MQ=255)
                     ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                              >  2:941144/1‑127 (MQ=255)
                      gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                           >  1:1062348/1‑145 (MQ=255)
                      gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                           <  2:1062348/145‑1 (MQ=255)
                            ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 <  1:1060649/149‑1 (MQ=255)
                            ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGACATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 <  4:145851/149‑1 (MQ=255)
                              cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                             >  2:833902/1‑119 (MQ=255)
                              cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                             <  1:833902/119‑1 (MQ=255)
                                     aataatCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                          >  1:1142402/1‑115 (MQ=255)
                                     aataatCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                          <  2:1142402/115‑1 (MQ=255)
                                         atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                    >  1:227608/1‑149 (MQ=255)
                                         atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                    <  4:377405/149‑1 (MQ=255)
                                            aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 >  1:443960/1‑133 (MQ=255)
                                            aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 <  2:443960/133‑1 (MQ=255)
                                            aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt                                                                                                                                         >  3:18397/1‑108 (MQ=255)
                                            aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt                                                                                                                                         <  4:18397/109‑2 (MQ=255)
                                                aaaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtata                                                                                             >  2:280100/1‑149 (MQ=255)
                                                 aaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                            >  1:96153/1‑149 (MQ=255)
                                                   gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt                                                                                                  <  4:299926/141‑1 (MQ=255)
                                                   gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt                                                                                                  >  3:299926/1‑141 (MQ=255)
                                                      ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                       <  2:227608/149‑1 (MQ=255)
                                                               tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTc                                                                                                                                                  <  1:39084/81‑1 (MQ=255)
                                                               tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTc                                                                                                                                                  >  2:39084/1‑81 (MQ=255)
                                                                        aGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                    <  2:112852/118‑1 (MQ=255)
                                                                        aGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                    >  1:112852/1‑118 (MQ=255)
                                                                           cATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAtttt                                                                                         <  2:1282153/126‑1 (MQ=255)
                                                                           cATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATAtttt                                                                                         >  1:1282153/1‑126 (MQ=255)
                                                                             tAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat                                                                <  2:96153/149‑1 (MQ=255)
                                                                             tAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATAATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat                                                                >  4:961368/1‑149 (MQ=255)
                                                                              aaTGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 <  3:357685/99‑1 (MQ=255)
                                                                              aaTGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                 >  4:357685/1‑99 (MQ=255)
                                                                                 gTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCt                                                                                     >  1:638012/1‑124 (MQ=255)
                                                                                 gTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCt                                                                                     <  2:638012/124‑1 (MQ=255)
                                                                                   aCTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattat                                                          <  1:280100/149‑1 (MQ=255)
                                                                                     tAAAAGCGTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa                                                                                                                       >  2:1200250/1‑86 (MQ=255)
                                                                                     tAAAAGCGTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa                                                                                                                       <  1:1200250/86‑1 (MQ=255)
                                                                                     tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                       >  2:532387/1‑118 (MQ=255)
                                                                                     tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                       <  1:532387/118‑1 (MQ=255)
                                                                                     tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                            >  2:629187/1‑145 (MQ=255)
                                                                                     tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                            <  1:629187/145‑1 (MQ=255)
                                                                                      aaaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattattat                                                       >  1:546639/1‑149 (MQ=255)
                                                                                        aaGCATAACGCTTACTTTCGACATCTTTAATATATGTCATCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGa                                                     <  3:961368/149‑1 (MQ=255)
                                                                                         aGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                    <  2:546639/149‑1 (MQ=255)
                                                                                            aTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                       <  4:839310/111‑1 (MQ=255)
                                                                                            aTAACGCTTACGTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                       >  3:839310/1‑111 (MQ=255)
                                                                                                       tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGc                                              <  2:1094146/141‑1 (MQ=255)
                                                                                                       tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGc                                              >  1:1094146/1‑141 (MQ=255)
                                                                                                                 tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc                                                                                      >  1:413058/1‑91 (MQ=255)
                                                                                                                 tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc                                                                                      <  2:413058/91‑1 (MQ=255)
                                                                                                                 tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTaa                            >  4:11393/1‑149 (MQ=255)
                                                                                                                     atatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata                        <  2:1986/149‑1 (MQ=255)
                                                                                                                        tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc                                                                                                               >  1:49685/1‑59 (MQ=255)
                                                                                                                        tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatc                                                                                                               <  2:49685/59‑1 (MQ=255)
                                                                                                                            tCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa                                                >  3:366813/1‑118 (MQ=255)
                                                                                                                            tCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa                                                <  4:366813/118‑1 (MQ=255)
                                                                                                                                     catcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                              <  1:83201/63‑1 (MQ=255)
                                                                                                                                     catcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                              >  2:83201/1‑63 (MQ=255)
                                                                                                                                      atcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg                                            >  1:183207/1‑112 (MQ=255)
                                                                                                                                      atcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg                                            <  2:183207/112‑1 (MQ=255)
                                                                                                                                         atGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCt    >  2:1053255/1‑149 (MQ=255)
                                                                                                                                           gCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTgg  >  3:812136/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
ACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGG  >  NC_000913/2102965‑2103252

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