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.G94.8% 297.6 / 10.8 117coding (1067/1089 nt)b3269putative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (58/52);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AAAACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGC  >  NC_000913/3418762‑3418907
                                                                              |                                                                        
aaaaCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTTCTGCATCAGAc                                                                           >  2:4073952/1‑76 (MQ=255)
 aaaCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTTCCGGCACCGTGCTGGATCAGACg                                                                         >  1:2341446/1‑76 (MQ=255)
 aaaCTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGGATCAGACg                                                                         <  1:5557365/76‑1 (MQ=255)
  acctcttcTCTCGCCGCCGCTATTTGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACgg                                                                        <  1:2491240/74‑1 (MQ=255)
  aaCTCCCCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACgg                                                                        <  1:4284034/76‑1 (MQ=255)
   acccctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggg                                                                      <  2:4827971/73‑1 (MQ=255)
   acccctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggg                                                                      <  2:3711870/73‑1 (MQ=255)
   aCTCCTCTCTTGCCGCCGCTATTGGCTATCCCTATATGGTTTCGCTGTTTGCCGGCCCCGTGCTGAATCAGACggg                                                                      <  1:4054285/76‑1 (MQ=255)
    ctcctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACgggg                                                                     <  1:2069480/76‑2 (MQ=255)
    ctcctcTCTTGCCGCCGCTATTGGCTATACCGATATGGTTTCGCTGTTTGCCCGCACCGTGCTCAATCAGACaggg                                                                     >  2:4277726/1‑75 (MQ=255)
    cccctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACgggg                                                                     <  2:3867461/74‑2 (MQ=255)
     tcctcTCTTGCCGCCGCTATTGGCTATCCCGATATTGGTTCTCTGTTTTCCGGCACCGTGCTTAATAAGACggggc                                                                    >  1:3080227/1‑74 (MQ=255)
     tcctcTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACggggc                                                                    >  1:933305/1‑74 (MQ=255)
            ttGTCGCCGCTTTGGGCACTCCGGTTTTGGTTTGGTTGTTAGCCGGGACCGTGCTGAATCATACGGGGCAAGCCAt                                                             <  1:625782/76‑1 (MQ=255)
            ttGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTCATTGGCGGCCACCTGCTGAATCAGACGGCGCAAGcacg                                                             >  2:2606880/1‑73 (MQ=255)
             tGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATc                                                            >  1:4680320/1‑76 (MQ=255)
             tGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATc                                                            <  1:3801483/76‑1 (MQ=255)
              gccgccGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCg                                                           >  1:4780028/1‑76 (MQ=255)
              gccgccGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCg                                                           <  2:2424167/76‑1 (MQ=255)
              gccgccGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGCGCTGAATCAGACGGGGCAAGCCAGCg                                                           >  1:5079850/1‑76 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTTCCGGCCCCGCGCTGAATCAGACGGGGCAAGCCACCgaga                                                        >  2:5019738/1‑76 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGGAAGCCATCgaga                                                        >  1:1971621/1‑76 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCCAGCCATCgcgc                                                        >  2:5387131/1‑73 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCgaga                                                        <  1:510065/76‑1 (MQ=255)
                 gccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCAGCgaga                                                        >  1:5292338/1‑76 (MQ=255)
                  ccgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGAc                                                       >  1:1650450/1‑76 (MQ=255)
                   cgcTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACg                                                      >  1:5271264/1‑76 (MQ=255)
                      tATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATa                                                   >  2:2120518/1‑76 (MQ=255)
                        ttGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCTATCGAGACGATAGc                                                 >  2:3958565/1‑76 (MQ=255)
                         tGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGcc                                                >  1:3641096/1‑76 (MQ=255)
                          ggCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCACCGCGACGATAAccc                                               >  2:1071121/1‑75 (MQ=255)
                           gatATCCCGATGTGGTTTCGCTGTTTGCCGCCTCCGTGCGGAATCAGACGGGGCAGGCCATCGAGACGATAGCCAt                                              <  2:2593065/74‑1 (MQ=255)
                           gCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCac                                              >  2:5060883/1‑75 (MQ=255)
                           gCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCAt                                              >  1:5165117/1‑76 (MQ=255)
                           gCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCAt                                              >  1:1070961/1‑76 (MQ=255)
                            ctgtCCCGATATGGTTTCGCTGTTGGCCGGCACGGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATg                                             <  2:1085883/73‑1 (MQ=255)
                             tATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGa                                            >  1:4265972/1‑76 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAcc                                          <  1:303156/76‑1 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCCTGAcc                                          >  2:5158915/1‑76 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAcc                                          >  2:4170172/1‑76 (MQ=255)
                               tCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGAcc                                          <  2:1070961/76‑1 (MQ=255)
                                cccGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCa                                         <  1:3175731/76‑1 (MQ=255)
                                  cGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATg                                       >  1:3996125/1‑76 (MQ=255)
                                  agATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGAAGGGGCAAGACATCGAGACGATAGCCACGCCCATg                                       >  2:4513627/2‑76 (MQ=255)
                                   gATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGt                                      <  1:4203766/76‑1 (MQ=255)
                                     tatGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCg                                    >  1:4615837/1‑76 (MQ=255)
                                     tatGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCg                                    >  1:1335641/1‑76 (MQ=255)
                                     tatGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCg                                    <  2:5165116/76‑1 (MQ=255)
                                      atGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCgg                                   >  1:2941039/1‑76 (MQ=255)
                                      atGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCgg                                   >  2:5374545/1‑76 (MQ=255)
                                       tGGTTTCGTTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGt                                  <  1:2802536/76‑1 (MQ=255)
                                        ggTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCCTGCCgggc                                 >  2:5747526/1‑74 (MQ=255)
                                          tttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTa                               <  1:2120518/76‑1 (MQ=255)
                                          tttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTa                               >  2:2058987/1‑76 (MQ=255)
                                           ttCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTAt                              <  1:4769262/76‑1 (MQ=255)
                                            tCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATc                             <  2:4643429/76‑1 (MQ=255)
                                            tCGCTGTTGGCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATc                             <  2:2041971/76‑1 (MQ=255)
                                             cGCTGTTTGCCGGCCCCGGGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCt                            <  1:443624/76‑1 (MQ=255)
                                             cGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCt                            <  1:4665102/76‑1 (MQ=255)
                                              gCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGCTAGCCATGACCATGTCGGTCTCTCTg                           >  2:1193206/1‑76 (MQ=255)
                                              gCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTg                           <  1:3413356/76‑1 (MQ=255)
                                               cTGTTTGCCGGCACCGTGCTGAATCGAACGGGGAAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGa                          <  1:2887596/76‑1 (MQ=255)
                                                tGTTTGCGGGCGCCGTGCTGAATCAGGGGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         <  2:613360/76‑1 (MQ=255)
                                                tGTTTGCCGGCACGGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         <  1:2338267/76‑1 (MQ=255)
                                                tGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         >  1:4850278/1‑76 (MQ=255)
                                                tGTTGGCGGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         <  2:1479310/76‑1 (MQ=255)
                                                ggTTTGCGGGGACGGTGGTGAATCAGACGGGGCACGCGATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAt                         <  2:865675/75‑1 (MQ=255)
                                                 gTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGAtt                        >  1:780100/1‑76 (MQ=255)
                                                  tttGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTa                       <  1:2824831/76‑1 (MQ=255)
                                                  tttGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTa                       <  1:5379399/76‑1 (MQ=255)
                                                   ttGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTAt                      >  1:2707873/1‑76 (MQ=255)
                                                   ttGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTAt                      <  1:2392438/76‑1 (MQ=255)
                                                   ttGCCGCCACCGGGCTGAATCAGCCGGGGAAAGCAATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTAt                      <  1:4073953/76‑1 (MQ=255)
                                                    tGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATc                     >  1:5602399/1‑76 (MQ=255)
                                                    tGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATc                     <  1:3312627/76‑1 (MQ=255)
                                                     gCCGGCCCCGTGGTGAACCAGAGGGGGCAAGCAATCGGGAAGATAGCCATGACCATGTCGGGCTATCTGATTATCa                    <  2:2491691/76‑1 (MQ=255)
                                                     gCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTAtaa                    <  2:3474288/76‑3 (MQ=255)
                                                     gCCGGCCCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCa                    <  2:2459130/76‑1 (MQ=255)
                                                      ccGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAg                   <  1:1596336/76‑1 (MQ=255)
                                                         gCACCGTGGTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCt                <  2:5677932/76‑1 (MQ=255)
                                                         gCACCGTGCTGAATCAGCCGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCt                <  1:3167163/76‑1 (MQ=255)
                                                         gCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCt                >  1:2553741/1‑76 (MQ=255)
                                                         gCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCt                <  1:3859855/76‑1 (MQ=255)
                                                          cACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTg               <  1:5447188/76‑1 (MQ=255)
                                                           aCCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTgc              >  2:661402/1‑75 (MQ=255)
                                                            ccGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCACCCTGAc             >  2:2965920/1‑76 (MQ=255)
                                                            cagTGCTGAACCAGCCAGGACAAGCCACCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGAc             <  1:362390/74‑1 (MQ=255)
                                                              gTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTa           <  1:4033039/76‑1 (MQ=255)
                                                              gTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCCATCTGATTATCAGCCTGAcct           >  2:2135414/1‑74 (MQ=255)
                                                               tGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAt          >  1:2766504/1‑76 (MQ=255)
                                                                gCTGAATCAGACGGGGCAAGCCATAGAGACGAGAGCCATGACCATGTCGGTCTATCTGATTATCAGCCCGACTAtc         >  1:835707/1‑76 (MQ=255)
                                                                 cTGAATCAGAGGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtct        <  1:657361/76‑1 (MQ=255)
                                                                  tGAATCAGGCGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtctc       <  2:4188355/76‑1 (MQ=255)
                                                                  tGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtctc       >  1:616927/1‑76 (MQ=255)
                                                                  tGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGGCTATCTGATTATCAGCCTGACTAtctc       >  1:5080429/1‑76 (MQ=255)
                                                                  tGAACCAGACGGGGAGAGCAATCAGGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTAtctc       <  1:1193206/76‑1 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGCTTATCAGCCTGACTATCTCg      >  2:1249806/1‑76 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTCTCTCg      >  1:3772463/1‑76 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCg      >  1:363888/1‑76 (MQ=255)
                                                                   gAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCg      >  2:1307477/1‑76 (MQ=255)
                                                                              |                                                                        
AAAACTCCTCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGC  >  NC_000913/3418762‑3418907

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