Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,418,840 +G coding (1070/1089 nt) b3269 → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,8371.G98.1% 274.1 / ‑1.4 105coding (1067/1089 nt)b3269putative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (56/47);  total (57/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

AAAACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGC  >  NC_000913/3418762‑3418907
                                                                             |                                                                       
aaaaCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGAc                                                                          <  1:3936196/76‑1 (MQ=255)
 aaaCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACg                                                                        <  1:3976969/76‑1 (MQ=255)
 aaaCTCCTCTCTTCCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACg                                                                        <  2:1433749/76‑1 (MQ=255)
   aCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggg                                                                      <  2:3059214/76‑1 (MQ=255)
   aCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggg                                                                      <  2:2894807/76‑1 (MQ=255)
     tcctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggggc                                                                    <  1:2716957/76‑3 (MQ=255)
      cctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACCGGGCa                                                                    >  2:3646570/1‑76 (MQ=255)
      cctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggggca                                                                   <  1:3330966/76‑4 (MQ=255)
      cctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggggca                                                                   <  2:3003414/76‑4 (MQ=255)
         ctctTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGc                                                                <  1:735530/76‑1 (MQ=255)
         ctctTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGc                                                                >  1:1976623/1‑76 (MQ=255)
           ctTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCa                                                              <  1:3327557/76‑1 (MQ=255)
            ttGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCAt                                                             >  2:2957993/1‑76 (MQ=255)
            ttGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCAt                                                             >  1:3327879/1‑76 (MQ=255)
            ttGCCCCCGCTATGGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCCCCGTGCTGAATCAGCCGGGGCAAGCCAt                                                             <  1:607020/76‑1 (MQ=255)
             tGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATc                                                            >  2:3250040/1‑76 (MQ=255)
              gccgccGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCg                                                           >  1:4538564/1‑76 (MQ=255)
               ccgccgCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCga                                                          <  1:4705581/76‑1 (MQ=255)
               ccgccgCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCga                                                          <  1:4115197/76‑1 (MQ=255)
               ccgccgCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCga                                                          <  2:2667094/76‑1 (MQ=255)
                cgcggcTATTGGATGTCCCGATGTGGTTTCGCTGTTTGCGGGCATCGTGCTGAATGAGACGGGGCACGCCATCggg                                                         <  2:842904/76‑3 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCgaga                                                        >  2:3369226/1‑76 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCACCgaga                                                        >  2:2482802/1‑76 (MQ=255)
                  ccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGAc                                                       >  1:3192357/1‑76 (MQ=255)
                  ccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGAc                                                       >  1:1113656/1‑76 (MQ=255)
                  ccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGAc                                                       <  2:4157532/76‑1 (MQ=255)
                  ccgcTATTGGCGATCCCGATATGGTTTCGCTGTTTGCCGGCACCGGGCTGAATCAGACGGGGCAAGCCATCGAGAc                                                       <  2:4452220/76‑1 (MQ=255)
                   cgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACg                                                      >  1:3141729/1‑76 (MQ=255)
                   cgcCATTGGCTACCCCGACATGGTTTCGCTGTTTGCCGGCCCCGAGCTGAATCAGCCGGGGCAAGCCACCGAGACg                                                      <  1:1153552/76‑1 (MQ=255)
                    gcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCCTCGAGACGa                                                     >  2:3418788/1‑76 (MQ=255)
                    gcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGa                                                     <  1:1301446/76‑1 (MQ=255)
                    gcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGa                                                     <  1:966795/76‑1 (MQ=255)
                      tATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATa                                                   >  1:1626827/1‑76 (MQ=255)
                      tATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATa                                                   <  1:3770515/76‑1 (MQ=255)
                        ttGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGc                                                 >  1:4411555/1‑76 (MQ=255)
                        ttGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGc                                                 >  1:1974530/1‑76 (MQ=255)
                         tGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGCCGATAccc                                                >  2:731352/1‑76 (MQ=255)
                         tGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGcc                                                <  1:933759/76‑1 (MQ=255)
                         tGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGcc                                                <  1:2429143/76‑1 (MQ=255)
                         tGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGcc                                                >  1:4163892/1‑76 (MQ=255)
                             tATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGa                                            <  1:2787706/76‑1 (MQ=255)
                             tATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGa                                            >  2:592337/1‑76 (MQ=255)
                              aTCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAc                                           >  1:2236640/1‑76 (MQ=255)
                               tCCCGATATGTTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAcc                                          <  1:1846333/76‑1 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCCTGAcc                                          >  2:4255969/1‑76 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAcc                                          <  1:3397802/76‑1 (MQ=255)
                                cccGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCa                                         >  1:4814623/1‑76 (MQ=255)
                                 ccGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCAt                                        >  1:88350/1‑76 (MQ=255)
                                 ccGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATACCCATGcccct                                        >  2:4536588/1‑74 (MQ=255)
                                  cGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGTTAGCCATGACCATg                                       >  1:4817051/1‑76 (MQ=255)
                                  cGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATg                                       >  2:4729720/1‑76 (MQ=255)
                                   gATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGt                                      >  1:2663102/1‑76 (MQ=255)
                                    atatGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTc                                     <  2:2437963/76‑1 (MQ=255)
                                    atatGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTc                                     >  1:2477877/1‑76 (MQ=255)
                                       tGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGt                                  >  2:1182274/1‑76 (MQ=255)
                                       tGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGt                                  >  2:3105218/1‑76 (MQ=255)
                                        ggTTGCGCTGTTTGCCGGCCCCGTGCGGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTc                                 <  1:3747643/76‑1 (MQ=255)
                                          tttCGTGGTTGGCGGGCACCGGGCTGAATCGGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTa                               <  2:568097/76‑1 (MQ=255)
                                          tttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTa                               >  1:4308768/1‑76 (MQ=255)
                                          tttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTa                               <  2:4369604/76‑1 (MQ=255)
                                           ttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTAt                              >  1:515195/1‑76 (MQ=255)
                                             cGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCt                            <  1:4708197/76‑1 (MQ=255)
                                              gCTGTTTGCCGGCACCGGGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTGCGTCTATCTg                           >  1:4104296/1‑76 (MQ=255)
                                              ccTGTTTGCCGGCACCGTGCTGACTCAGACGGGGCAAGCCATCGAGCCGATAGCCCAGACCATTTCTGGCTATCtt                           >  2:4872602/2‑75 (MQ=255)
                                               cTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGa                          >  1:3528333/1‑76 (MQ=255)
                                               cTGTTTGCCGGCACCGCGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCACGACCCTGTCGGTCTATCTGa                          >  2:705814/1‑76 (MQ=255)
                                                tGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         <  2:2203027/76‑1 (MQ=255)
                                                tGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         >  2:1384348/1‑76 (MQ=255)
                                                tGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         >  1:2253071/1‑76 (MQ=255)
                                                 ttttGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAtt                        <  1:4396006/75‑1 (MQ=255)
                                                 gTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAtt                        >  2:3176029/1‑76 (MQ=255)
                                                  tttGCCCACACCGTGCTGAAGCAGACGGGGCAAGCCACCGAGACGATAGCCACGCCCATGAAGGTCTATCTGATTa                       >  2:2988762/1‑76 (MQ=255)
                                                   ntGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGATACGATAGCCATGACCATGTCGGTCTATCTGATTAt                      >  1:3537651/2‑76 (MQ=255)
                                                    tGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATc                     >  1:2840213/1‑76 (MQ=255)
                                                    tGCCGGCACCGGGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATc                     <  1:3590275/76‑1 (MQ=255)
                                                     gCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCa                    <  1:1323925/76‑1 (MQ=255)
                                                     gCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCa                    >  1:1507830/1‑76 (MQ=255)
                                                      ccGCCACCGTGCTAAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAg                   <  1:4819027/76‑1 (MQ=255)
                                                       cGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTCTCTGATTATCAGc                  >  2:4852195/1‑76 (MQ=255)
                                                       cGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGc                  >  1:4417546/1‑76 (MQ=255)
                                                        ggCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGcc                 >  2:4265277/1‑76 (MQ=255)
                                                        ggCACCGTGCTGAATAAGACGGGGCAGGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGcc                 <  2:4167608/76‑1 (MQ=255)
                                                           cccGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGa              <  1:124193/75‑1 (MQ=255)
                                                           aCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTgc              >  2:4660020/1‑75 (MQ=255)
                                                           aCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGa              >  1:4441586/1‑76 (MQ=255)
                                                             cGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACt            <  2:1507829/76‑1 (MQ=255)
                                                              gTGCTGAATCAGAGGGGGAACGCCATCGAGACGATAGCCATGACCATGTCGGTCNATCTGATTATCAGCCTGACTa           <  1:2074386/76‑1 (MQ=255)
                                                               tGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAt          <  1:940613/76‑1 (MQ=255)
                                                               tGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAt          >  2:4581967/1‑76 (MQ=255)
                                                                gCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtc         <  2:3557662/76‑1 (MQ=255)
                                                                 cTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtct        <  1:3359639/76‑1 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGGCTATCTCg      <  2:2490915/76‑1 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCg      >  2:4619511/1‑76 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCg      >  1:3492534/1‑76 (MQ=255)
                                                                    aaTCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgc     >  2:2020757/1‑76 (MQ=255)
                                                                     agcAGACGGGGGAAGCCCTCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgct    <  2:3706887/74‑1 (MQ=255)
                                                                     aNCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgct    <  2:3277264/76‑1 (MQ=255)
                                                                      tCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgctg   <  1:4375713/76‑1 (MQ=255)
                                                                       cAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgctgc  >  2:1392466/1‑76 (MQ=255)
                                                                       cAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCgctgc  <  1:731352/76‑1 (MQ=255)
                                                                             |                                                                       
AAAACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGC  >  NC_000913/3418762‑3418907

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