Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,458,416 T→A 7.1% I26I (ATT→ATA gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,4160TA 7.1% 172.4 / 4.0 56I26I (ATT→ATAgspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base A (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TTTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGC  >  NC_000913/3458269‑3458558
                                                                                                                                                   |                                                                                                                                              
tttGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGatta                                                                                                                                               >  8:57260/1‑149 (MQ=255)
   gCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGattatt                                                                                                                                             <  5:3297/148‑1 (MQ=255)
    ccGGTGCTGGATGAAAACCCGTTCAACATTCAGGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGACTTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattt                                                                                                                                            <  6:148133/148‑1 (MQ=255)
    ccGGTGCTGGATGAAAACACGTTCACCAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCAAAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattt                                                                                                                                            <  6:151557/148‑1 (MQ=255)
             gATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGATAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACTACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGcc                                                                                                                                    <  8:41313/147‑1 (MQ=255)
             gATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGcc                                                                                                                                    <  7:243928/147‑1 (MQ=255)
              aTGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGa                                                                                                                                  >  4:181330/1‑148 (MQ=255)
               tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTTGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATa                                                                                                                                >  5:111927/1‑149 (MQ=255)
               tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGccctct                                                                                                                                <  1:12573/149‑5 (MQ=255)
                  aaaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGc                                                                                                                             >  2:231764/1‑149 (MQ=255)
                   aaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGc                                                                                                                             >  6:225898/1‑148 (MQ=255)
                         ttCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGa                                                                                                                      <  1:231764/149‑1 (MQ=255)
                            accAGTCCCGCCCCTGCGCCAGCGTAACGGGGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCg                                                                                                                   <  6:126155/147‑1 (MQ=255)
                                            cgcCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTc                                                                                                     <  6:111927/147‑1 (MQ=255)
                                                 tCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGata                                                                                                        >  7:158037/1‑139 (MQ=255)
                                                 tCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGata                                                                                                        <  8:158037/139‑1 (MQ=255)
                                                    tCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGt                                                                                           >  8:3446/1‑149 (MQ=255)
                                                        ggTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGc                                                                                       <  7:57260/149‑1 (MQ=255)
                                                            aGGCCTTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGACTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGGCAGGTTGGCCAg                                                                                    <  4:156045/148‑1 (MQ=255)
                                                            aGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGt                                                                                   <  6:9829/149‑1 (MQ=255)
                                                              gCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGc                                                                                 >  5:200664/1‑149 (MQ=255)
                                                                      atGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                         >  1:241367/1‑149 (MQ=255)
                                                                       tGAGAATTCACTCACCGTACCCCGCGTGTTGGGCGCTGGCACAACGAATTGGTTATCTGTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                         <  4:110637/148‑1 (MQ=255)
                                                                       tGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATg                                                                        >  2:237803/1‑149 (MQ=255)
                                                                        gagaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCa                                                                                                    >  3:9477/1‑120 (MQ=255)
                                                                        gagaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCa                                                                                                    <  4:9477/120‑1 (MQ=255)
                                                                        gagaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGa                                                                       >  6:235349/1‑149 (MQ=255)
                                                                            aTTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAg                                                                   >  8:202816/1‑149 (MQ=255)
                                                                                actcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                               >  7:95603/1‑149 (MQ=255)
                                                                                  tcacCGTACCCCGCCAGTTGGGCGCTGGGACAACGAATTGGTTATCTCTATTCAGAGCGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCATGTTGGCCAGTGCCAGACAATGACCAGCTTGt                                                              >  1:96608/1‑148 (MQ=255)
                                                                                  tcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGt                                                              >  8:44537/1‑148 (MQ=255)
                                                                                    acCGTACCCCTCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTca                                                            >  3:69902/1‑148 (MQ=255)
                                                                                              gCCAGTTGGGCGCTGGCACAACGAATTGTTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCAGATACGCCATTCGAGCGGTTACTCGATATTACACGTCATGTTGGCCAGTGCCAGACCATGACAAGCTTGTCACAGGCTGAttt                                                 >  5:86710/1‑149 (MQ=255)
                                                                                                  gTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                             <  8:95603/149‑1 (MQ=255)
                                                                                                        cgcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                                       <  5:225898/149‑1 (MQ=255)
                                                                                                        cgcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                                       >  8:57871/1‑149 (MQ=255)
                                                                                                          cTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCt                                     <  7:44537/149‑1 (MQ=255)
                                                                                                          cTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCt                                     >  3:24594/1‑149 (MQ=255)
                                                                                                               acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                                <  1:52721/117‑1 (MQ=255)
                                                                                                               acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                                >  2:52721/1‑117 (MQ=255)
                                                                                                                 aaCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                              <  4:69902/149‑1 (MQ=255)
                                                                                                                  aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                             <  3:181330/149‑1 (MQ=255)
                                                                                                                    gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                           >  8:243285/1‑149 (MQ=255)
                                                                                                                         ggTTATCTCTATTCAGCGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc                       <  2:45015/148‑1 (MQ=255)
                                                                                                                         ggTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc                       <  2:241367/148‑1 (MQ=255)
                                                                                                                               ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                >  7:95737/1‑149 (MQ=255)
                                                                                                                                tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                >  2:97832/1‑148 (MQ=255)
                                                                                                                                tctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                >  3:117719/1‑148 (MQ=255)
                                                                                                                                    ttCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTAc                                                                                                              >  5:62600/1‑50 (MQ=255)
                                                                                                                                    ttCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTAc                                                                                                              <  6:62600/50‑1 (MQ=255)
                                                                                                                                      cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGt         >  6:178795/1‑149 (MQ=255)
                                                                                                                                       agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                                >  5:8994/1‑109 (MQ=255)
                                                                                                                                       agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                                <  6:8994/109‑1 (MQ=255)
                                                                                                                                       agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGt         <  7:202816/148‑1 (MQ=255)
                                                                                                                                       agagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGt         >  8:77579/1‑148 (MQ=255)
                                                                                                                                           ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGt                            >  3:110994/1‑125 (MQ=255)
                                                                                                                                           ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGt                            <  4:110994/125‑1 (MQ=255)
                                                                                                                                           ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTc    >  5:139012/1‑149 (MQ=255)
                                                                                                                                            gCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCg   >  3:217300/1‑149 (MQ=255)
                                                                                                                                              gagaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGc  >  2:89186/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TTTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGC  >  NC_000913/3458269‑3458558

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: