Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,401,970 G→A 7.3% R187H (CGC→CAC)  spr → predicted peptidase, outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,401,9700GA 7.3% 276.4 / 8.2 96R187H (CGC→CAC) sprpredicted peptidase, outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base A (4/3);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCT  >  Exported/1401901‑1402034
                                                                     |                                                                
ttattaTTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCg                                                                  >  1:5171596/1‑70 (MQ=255)
  attattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCg                                                                  <  1:400803/68‑1 (MQ=255)
  attattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCg                                                                  <  1:5067055/68‑1 (MQ=255)
  attattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCa                                                                <  1:5522145/70‑1 (MQ=255)
   ttattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCg                                                                  <  1:2026893/67‑1 (MQ=255)
    tattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGc                                                              >  1:5627303/1‑70 (MQ=255)
     attTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCt                                                             >  1:4052654/1‑70 (MQ=255)
      tttCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGc                                                              >  1:2278391/1‑68 (MQ=255)
         ccAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCtaata                                                         <  1:4694354/70‑1 (MQ=255)
         ccAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCtaat                                                          <  1:4057848/69‑1 (MQ=255)
              atgaatgaACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGc                                                                 <  1:4080478/57‑1 (MQ=255)
              atgaatgaACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGt                                                   >  1:370094/1‑71 (MQ=255)
              atgaatgaACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGt                                                   <  1:4046752/71‑1 (MQ=255)
              atgaatgaACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGt                                                   <  1:1971114/71‑1 (MQ=255)
                 aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGc                                                                 >  1:3517906/1‑54 (MQ=255)
                 aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                >  1:1020844/1‑71 (MQ=255)
                 aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                >  1:4318380/1‑71 (MQ=255)
                 aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                <  1:1904570/71‑1 (MQ=255)
                  atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTgg                                               <  1:1897319/71‑1 (MQ=255)
                  atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                >  1:2653920/1‑70 (MQ=255)
                  atgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                <  1:3764505/70‑1 (MQ=255)
                    gaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCtaat                                                          >  1:3868052/1‑58 (MQ=255)
                    gaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGAt                                             >  1:5552251/1‑71 (MQ=255)
                     aaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAAc                                                      >  1:1946384/1‑61 (MQ=255)
                     aaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGAt                                             >  1:5403289/1‑70 (MQ=255)
                      aCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATg                                            <  1:1508703/70‑1 (MQ=255)
                        cGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGc                                           <  1:2619807/69‑1 (MQ=255)
                          tACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCTCAGCt                                                             <  1:4352983/49‑1 (MQ=255)
                           aCTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATc                                       <  1:5614551/70‑1 (MQ=255)
                           aCTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATc                                       <  1:4622786/70‑1 (MQ=255)
                            cTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATcc                                      >  1:2727684/1‑70 (MQ=255)
                              ggAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCtt                                   >  1:3644943/1‑71 (MQ=255)
                               gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGa                                              >  1:3101296/1‑59 (MQ=255)
                               gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCtt                                   <  1:893403/70‑1 (MQ=255)
                               gaagaaGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTg                                  <  1:3425478/71‑1 (MQ=255)
                               gaagaaGCGTTACAACGAAGCACGCCGGGTTCTAAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTg                                  >  1:4545274/1‑71 (MQ=255)
                                aagaagCGTTACACCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCtt                                   <  1:3739176/69‑1 (MQ=255)
                                aagaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGa                                              <  1:1021169/58‑1 (MQ=255)
                                aagaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGa                                              <  1:2463362/58‑1 (MQ=255)
                                 agaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCaa                                         >  1:4458021/1‑62 (MQ=255)
                                 agaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGa                                              <  1:4983359/57‑1 (MQ=255)
                                 agaagCGTTACAACGAAGCACGCCGGCTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGc                                >  1:4372010/1‑71 (MQ=255)
                                  gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGATg                                            >  1:3640984/1‑58 (MQ=255)
                                  gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGATg                                            >  1:1802759/1‑58 (MQ=255)
                                  gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGATg                                            >  1:2294106/1‑58 (MQ=255)
                                  gaagCGTTACAACGAAGCACGCCGGGTTCTCAGCCACAGCTAATAAACCGTTTGGATg                                            >  1:5461549/1‑58 (MQ=255)
                                   aagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTgg                                 <  1:3282401/68‑1 (MQ=255)
                                   aagCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTa                              <  1:3629167/71‑1 (MQ=255)
                                   aagCGTTACAACGAAGCACGCCGGGATCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGc                                <  1:498355/69‑1 (MQ=255)
                                    agCGTTACAACGAAGCATGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTAt                             >  1:295262/1‑71 (MQ=255)
                                      cGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATcc                           >  1:110898/1‑71 (MQ=255)
                                       gTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTa                              <  1:5545843/67‑1 (MQ=255)
                                       gTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTa                              <  1:5524170/67‑1 (MQ=255)
                                        ttACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATc                            >  1:2372435/1‑68 (MQ=255)
                                        ttACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCt                          >  1:4384731/1‑70 (MQ=255)
                                        ttACAACGAAGCACGACGGGTTCTCAGCCGCAGCtaa                                                           <  1:3004345/37‑1 (MQ=255)
                                         tACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATcc                           <  1:1125459/68‑1 (MQ=255)
                                         tACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATcc                           <  1:1736567/68‑1 (MQ=255)
                                          aCAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAAc                                                      >  1:1385809/1‑40 (MQ=255)
                                          aCAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAAc                                                      >  1:1667769/1‑40 (MQ=255)
                                          aCAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGa                                              >  1:1537548/1‑48 (MQ=255)
                                          aCAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGAc                       >  1:3015546/1‑71 (MQ=255)
                                           cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGAc                       <  1:4441657/70‑1 (MQ=255)
                                              cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCtt                                   <  1:2889378/55‑1 (MQ=255)
                                              cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCt                          >  1:4918430/1‑64 (MQ=255)
                                              cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGt                   >  1:933427/1‑71 (MQ=255)
                                              cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGt                   >  1:4370170/1‑71 (MQ=255)
                                              cGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGt                   >  1:2751229/1‑71 (MQ=255)
                                               gAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGtt                  <  1:3744221/71‑1 (MQ=255)
                                                aaGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTa                 <  1:371621/71‑1 (MQ=255)
                                                  gCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTg                                                <  1:2518243/38‑1 (MQ=255)
                                                   cACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACg                      <  1:2576564/63‑1 (MQ=255)
                                                   cACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGtt                  <  1:2588796/67‑1 (MQ=255)
                                                      gCCGGGTTCTCAGCCGCATCTAATAAACCGTTTGGATGCaa                                         >  1:4709240/1‑41 (MQ=255)
                                                      gCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCaa                                         >  1:3781743/1‑41 (MQ=255)
                                                      gCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATcc                           >  1:1952772/1‑55 (MQ=255)
                                                        cGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAg         >  1:2190010/1‑71 (MQ=255)
                                                        cGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAg         <  1:4386232/71‑1 (MQ=255)
                                                         gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCt                                    <  1:1352235/43‑1 (MQ=255)
                                                         gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATc                            <  1:1873743/51‑1 (MQ=255)
                                                         gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATc                            <  1:2792390/51‑1 (MQ=255)
                                                         gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATc                            <  1:512523/51‑1 (MQ=255)
                                                         gggTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAg         >  1:4233986/1‑70 (MQ=255)
                                                            ttCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACt              <  1:108118/62‑1 (MQ=255)
                                                             gatcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTg    <  1:3763804/69‑1 (MQ=255)
                                                             gatcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTg    <  1:11479/69‑1 (MQ=255)
                                                              ctcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGc   >  1:2702969/1‑71 (MQ=255)
                                                              ctcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGc   >  1:68852/1‑71 (MQ=255)
                                                               tcAGCCGCAGCTAATAAACCGTTTGGATGCAATccc                                     >  1:1589579/1‑36 (MQ=255)
                                                               tcAGCCGCAGCTAATAAACCGTTTGGATGCAATccc                                     >  1:3535552/1‑36 (MQ=255)
                                                                cAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGa                        >  1:897816/1‑48 (MQ=255)
                                                                cAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGa                        >  1:3957422/1‑48 (MQ=255)
                                                                 aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGa            >  1:2813579/1‑59 (MQ=255)
                                                                 aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt  <  1:4916333/69‑1 (MQ=255)
                                                                 aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt  <  1:817081/69‑1 (MQ=255)
                                                                 aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt  <  1:3431154/69‑1 (MQ=255)
                                                                 aGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCt  <  1:3895043/69‑1 (MQ=255)
                                                                  gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAAc               >  1:3516209/1‑55 (MQ=255)
                                                                  gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAAc               >  1:1923293/1‑55 (MQ=255)
                                                                  gCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAAc               >  1:4484611/1‑55 (MQ=255)
                                                                     |                                                                
TTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCT  >  Exported/1401901‑1402034

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