Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 346,248 (A)7→8 coding (90/672 nt) yahN ← amino acid exporter for proline, lysine, glutamate, homoserine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913346,2411.A97.9% 180.3 / ‑3.4 48coding (97/672 nt)yahNamino acid exporter for proline, lysine, glutamate, homoserine
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (25/22);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAGAGATCACAAAGTGTCTT  >  NC_000913/346091‑346384
                                                                                                                                                       |                                                                                                                                               
cacaCTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGAtt                                                                                                                                                  >  2:157609/1‑151 (MQ=255)
      gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaaaaa                                                                                                                                           <  1:720705/151‑1 (MQ=255)
      gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaaaaa                                                                                                                                           <  1:302102/151‑1 (MQ=255)
      gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaaaa                                                                                                                                            >  1:697788/1‑150 (MQ=255)
      gCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaaa                                                                                                                                             <  2:238251/149‑1 (MQ=255)
          aaTTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTaaaaaaaag                                                                                                                                        <  2:941622/150‑3 (MQ=255)
             tAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAAt                                                                                                                                     >  1:422416/1‑151 (MQ=255)
                   tGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGa                                                                                                                                >  1:659045/1‑150 (MQ=255)
                    gCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGaa                                                                                                                               <  2:442650/150‑1 (MQ=255)
                       aGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGt                                                                                                                           >  1:248003/1‑151 (MQ=255)
                        gACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTc                                                                                                                          <  1:645241/151‑1 (MQ=255)
                            aaacaaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTAc                                                                                                                      >  2:740000/1‑151 (MQ=255)
                                aaacCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTc                                                                                                                  >  1:281838/1‑151 (MQ=255)
                                  acccaacccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAg                                                                                                                <  1:698565/151‑1 (MQ=255)
                                     caacccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGt                                                                                                              >  1:393985/1‑150 (MQ=255)
                                      aacccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTa                                                                                                             >  1:923897/1‑150 (MQ=255)
                                       acccGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTaaa                                                                                                           <  2:494956/151‑1 (MQ=255)
                                          cGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACgg                                                                                                        >  2:167236/1‑151 (MQ=255)
                                           ggAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGc                                                                                                       <  2:360684/151‑1 (MQ=255)
                                             aaTAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCAt                                                                                                     <  1:56283/151‑1 (MQ=255)
                                                 aaaTGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCa                                                                                                  >  2:370716/1‑150 (MQ=255)
                                                      cATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGa                                                                                             >  2:336401/1‑150 (MQ=255)
                                                        tCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATcc                                                                                          >  2:830909/1‑151 (MQ=255)
                                                          gCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCAt                                                                                        >  1:780376/1‑151 (MQ=255)
                                                            ccAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCAt                                                                                        >  1:884878/1‑149 (MQ=255)
                                                            ccAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAg                                                                                      <  1:109788/151‑1 (MQ=255)
                                                             cAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGt                                                                                     >  1:186978/1‑151 (MQ=255)
                                                               gcgcCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGa                                                                                   <  2:165718/151‑1 (MQ=255)
                                                                   cACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTc                                                                               <  1:548054/151‑1 (MQ=255)
                                                                      gCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATc                                                                            <  1:44746/151‑1 (MQ=255)
                                                                             ccGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGATAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaat                                                                     >  1:946795/1‑151 (MQ=255)
                                                                              cGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaata                                                                    <  2:840329/151‑1 (MQ=255)
                                                                               ggTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCataaataa                                                                   <  2:407913/151‑1 (MQ=255)
                                                                                  cAGCACCCCTACGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATACGt                                                                 >  1:906399/1‑150 (MQ=255)
                                                                                    gCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGaa                                                              <  1:730948/151‑1 (MQ=255)
                                                                                       ccccTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGaacta                                                           <  1:610852/151‑1 (MQ=255)
                                                                                               cgTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGcatc                                                   <  1:196930/151‑1 (MQ=255)
                                                                                                gTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGcatc                                                   >  2:768848/1‑150 (MQ=255)
                                                                                                   gACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCatat                                               >  2:257338/1‑151 (MQ=255)
                                                                                                     ccGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATAttt                                             <  2:274336/151‑1 (MQ=255)
                                                                                                          aGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTAcc                                        >  1:994813/1‑151 (MQ=255)
                                                                                                           gCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTAcct                                       <  2:794526/151‑2 (MQ=255)
                                                                                                                   ggTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAaataat                               <  2:42597/151‑1 (MQ=255)
                                                                                                                   ggTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAaataat                               >  2:620874/1‑151 (MQ=255)
                                                                                                                   ggTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAaata                                 <  2:740219/149‑1 (MQ=255)
                                                                                                                   ggTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGGCCTACTGTCAGGTAAATTTCATTCAAAGGATCGCTAGCGCTTTCATCCATAAATAAGTGAACTACCTTCCCCATATGGCAACCataaataagc                              <  2:119390/151‑3 (MQ=255)
                                                                                                                       tGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCggg                           >  1:648456/1‑151 (MQ=255)
                                                                                                                        gTACTACCACAACGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAATCGGCATGCAAAGGATCCATAGAGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAaata                                 >  2:76751/1‑144 (MQ=255)
                                                                                                                        gTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGt                          <  1:737879/151‑1 (MQ=255)
                                                                                                                           ctacCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCg                       >  2:293016/1‑151 (MQ=255)
                                                                                                                            tacCACAAAGAGATTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCg                       >  1:239763/1‑150 (MQ=255)
                                                                                                                                     gagaTTGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAGAGATca             <  2:810316/151‑1 (MQ=255)
                                                                                                                                          tGGCTCCCGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAGAGATCACaaa         <  2:697707/150‑1 (MQ=255)
                                                                                                                                                 cGGATTAAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAGAGATCACAAAGTGTCtt  >  2:881274/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
CACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTTGTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAAAGGATCCATAGTGATTTCATCCATAAATAAGTGAACTAACTGCATCATATTTCTACCAAAAATAATCGGGTGCGAGAGAGATCACAAAGTGTCTT  >  NC_000913/346091‑346384

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