Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 108,835 A→G 7.1% Y186C (TAC→TGC)  secA → preprotein translocase subunit, ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913108,8350AG 7.1% 164.6 / 3.3 56Y186C (TAC→TGC) secApreprotein translocase subunit, ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base G (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGAAA  >  NC_000913/108692‑108981
                                                                                                                                               |                                                                                                                                                  
tGACGCCGAAAACACCCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTgc                                                                                                                                               >  7:207441/1‑149 (MQ=255)
     ccGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGac                                                                                                                                          <  4:72285/149‑1 (MQ=255)
         aaaCAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGTGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGacaaca                                                                                                                                      <  4:101238/149‑1 (MQ=255)
             acccGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCTGGCATTCGAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGc                                                                                                                                  <  8:28372/147‑1 (MQ=255)
                      cTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGccc                                                                                                                         >  3:205041/1‑149 (MQ=255)
                          ttGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTg                                                                                                                       >  6:238398/1‑147 (MQ=255)
                            gAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTGCCTGCGCGACAACATGGCGTTCAGCCCTgaag                                                                                                                    >  7:40564/1‑148 (MQ=255)
                            gAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTGCCTGCGCGACAACATGGCGTTCAGCCCTgaag                                                                                                                    >  7:36538/1‑148 (MQ=255)
                                tCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgt                                                                                                               >  4:63868/1‑149 (MQ=255)
                                 ccTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgt                                                                                                               <  7:143245/148‑1 (MQ=255)
                                   ttGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTa                                                                                                            >  8:251957/1‑149 (MQ=255)
                                   ttGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTa                                                                                                            >  4:229386/1‑149 (MQ=255)
                                       ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              >  2:175752/1‑143 (MQ=255)
                                       ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              >  6:131190/1‑143 (MQ=255)
                                       ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              <  5:131190/143‑1 (MQ=255)
                                       ccTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACgtg                                                                                                              <  1:175752/143‑1 (MQ=255)
                                             tGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCCCTTACGGGACGAACGACGAATACGGCTTTGACTACCTGCGCGACGACCTGGCGTTCCGCCCTGAAGAACGGGTACAGCGTAAAc                                                                                                  >  6:155998/1‑149 (MQ=255)
                                                  gTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCAc                                                                                             >  6:164059/1‑149 (MQ=255)
                                                     tCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTAt                                                                                          <  1:141704/149‑1 (MQ=255)
                                                          cTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCt                                                                                     >  7:137134/1‑149 (MQ=255)
                                                             ccGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCtt                                                                                    <  8:36446/147‑2 (MQ=255)
                                                               gggCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCtggtg                                                                                 >  1:64273/1‑148 (MQ=255)
                                                                  cATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGAc                                                                              >  8:88179/1‑148 (MQ=255)
                                                                       cAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg                                                                        >  3:178830/1‑149 (MQ=255)
                                                                       cAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTg                                                                        <  5:164059/149‑1 (MQ=255)
                                                                         tgcccGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGa                                                                      <  4:104792/146‑1 (MQ=255)
                                                                            ccGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTc                                                                   >  4:65062/1‑149 (MQ=255)
                                                                             cGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTc                                                                   >  2:41422/1‑148 (MQ=255)
                                                                               gCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCAt                                                                >  3:51166/1‑149 (MQ=255)
                                                                                cAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATc                                                               <  4:51166/149‑1 (MQ=255)
                                                                                 aaaGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATcc                                                              >  2:167519/1‑149 (MQ=255)
                                                                                        gAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGa                                                       >  5:80382/1‑149 (MQ=255)
                                                                                         aaGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGAt                                                      <  3:63868/149‑1 (MQ=255)
                                                                                            cTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGaa                                                   <  2:64273/149‑1 (MQ=255)
                                                                                                       gACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACAcc                                        >  1:6648/1‑149 (MQ=255)
                                                                                                        aCATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCg                                       <  7:251957/149‑1 (MQ=255)
                                                                                                         cATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGc                                      <  1:41422/149‑1 (MQ=255)
                                                                                                               ttACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCa                                >  1:66413/1‑149 (MQ=255)
                                                                                                               ttACGGTACGAACAACGAATACGGATTTGACTACCTGCGCGACAACATGGCGTTCAGACCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATTAAGCGAGTACACCGCTGATCa                                >  7:217120/1‑149 (MQ=255)
                                                                                                                 aCGGTACGCACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGACCGTGTACAGCGTAAACTGCACTATGCGCGGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCa                                <  5:155998/147‑1 (MQ=255)
                                                                                                                    gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATc                                 <  5:167420/143‑1 (MQ=255)
                                                                                                                    gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATc                                 >  6:167420/1‑143 (MQ=255)
                                                                                                                    gTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTcc                           >  5:89392/1‑149 (MQ=255)
                                                                                                                       cGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTccgg                         <  2:66413/148‑1 (MQ=255)
                                                                                                                         aacaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGt                                             >  7:7212/1‑126 (MQ=255)
                                                                                                                         aacaacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGt                                             <  8:7212/126‑1 (MQ=255)
                                                                                                                            aacGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTccggcccgg                    >  3:82848/1‑148 (MQ=255)
                                                                                                                                aaTACGGCTTTTACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTTGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaa               <  8:217120/149‑1 (MQ=255)
                                                                                                                                aaTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaa               <  3:65062/149‑1 (MQ=255)
                                                                                                                                 aTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCagaag              <  4:169976/149‑1 (MQ=255)
                                                                                                                                 aTACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCag                 >  7:117267/1‑146 (MQ=255)
                                                                                                                                    cGGCTTTGACTACCTGCGCGACAACCTGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGAc            <  7:88179/148‑1 (MQ=255)
                                                                                                                                          tGACTGCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTc       <  8:36538/147‑1 (MQ=255)
                                                                                                                                          tGACTGCCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTc       <  8:40564/147‑1 (MQ=255)
                                                                                                                                             cTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGaaa  >  7:251262/1‑149 (MQ=255)
                                                                                                                                             cTACCTGCGCGACAACATGGCGGTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGaaa  <  5:130218/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCCGGGCATGCCAGCACCGGCAAAGCGCGAAGCTTACGCAGCTGACATCACTTACGGTACGAACAACGAATACGGCTTTGACTACCTGCGCGACAACATGGCGTTCAGCCCTGAAGAACGTGTACAGCGTAAACTGCACTATGCGCTGGTGGACGAAGTGGACTCCATCCTGATCGATGAAGCGCGTACACCGCTGATCATTTCCGGCCCGGCAGAAGACAGCTCGGAAA  >  NC_000913/108692‑108981

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