Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 9.6% intergenic (+32/+387) hdeD → / ← arrS acid‑resistance membrane protein/small regulatory RNA ArrS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,657,5990TA 9.6% 99.8 / 5.6 42intergenic (+32/+387)hdeD/arrSacid‑resistance membrane protein/small regulatory RNA ArrS
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-02

GCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGG  >  NC_000913/3657457‑3657735
                                                                                                                                              |                                                                                                                                        
gCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGCTCATAGTCACt                                                                                                                                          >  1:565035/1‑143 (MQ=255)
  aaCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTg                                                                                                                                        >  1:113999/1‑143 (MQ=255)
   aCGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATTGTCACTTGa                                                                                                                                       >  1:203400/1‑143 (MQ=255)
           ggTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                                                               <  1:429693/143‑1 (MQ=255)
                     aCGCTGGTTTCCACACTGGTGCGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggt                                                                                                                     <  1:451957/143‑4 (MQ=39)
                       gCTGGTTTCCACACTGGTGGGAATTGAGCTGGTATTTGGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                                                   <  2:125310/143‑1 (MQ=255)
                        cTGGTTTCCGCTGTGTTGGGCATTGAACGGTTATTTGTCGCCGCCAGTCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCGCACGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                                                                  <  2:313860/143‑1 (MQ=255)
                           gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACACCGGTTGTCACCCGGCACATAGTCACTTGATGTGACTGTGATCCGATTAATAcc                                                                                                               >  2:432306/1‑142 (MQ=255)
                           gTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCAAGCCTGTTCAGCTTCGCCAGGTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTAACCCGGATCATAGTGGCTTGATGTGACTATGATGAGATTGATAct                                                                                                               >  1:60662/1‑143 (MQ=255)
                                cacaCTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACtctctc                                                                                                          >  1:166424/1‑143 (MQ=255)
                                 acacTGGTGGGACTTGGACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGTCTATGATCCGATTAATACTCTCTcc                                                                                                         <  2:64953/143‑1 (MQ=255)
                                   acTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGc                                                                                                       >  2:198679/1‑143 (MQ=255)
                                        tggGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                                                                  >  1:215741/1‑143 (MQ=255)
                                         gggAATTGGACTGATGGTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                                                 <  1:458501/143‑1 (MQ=255)
                                              ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTTACTATGATCCGATTAATACTCTCTCCGCTACGCAGcgtt                                                                                            >  1:618738/1‑143 (MQ=255)
                                                 agctgggatTTGGCGCCGCCAGCCTGTGCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                                                         <  2:205188/136‑1 (MQ=255)
                                                     gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTAACACGGTTCGCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                                                                     <  1:305509/143‑1 (MQ=255)
                                                      atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATcg                                                                                    >  1:509974/1‑142 (MQ=255)
                                                         tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                                                 >  2:585142/1‑143 (MQ=255)
                                                          ttggcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTg                                                                                <  1:388425/140‑1 (MQ=255)
                                                                          gTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                                                                <  1:263159/143‑1 (MQ=255)
                                                                                 ttCGCCAGTTTATTCCTTAAGCAGCAATAATTACCCCGTTTGTCACCCTGATCATAGTCACTTGATGTGACTATTATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                                                         <  2:657122/143‑1 (MQ=255)
                                                                                 gtggCCGGTTTGTTCGTTACGTAGCAATAATTACCCCGGTTGTCACTCGGATCATAGTCATTTGATGTGACTATGATCCGATTAATCCTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                                                         <  1:138676/140‑1 (MQ=255)
                                                                                    gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGTTGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAGCAATTGCGCACTATCATTGAAATAATTACCTg                                                      >  1:593747/1‑143 (MQ=255)
                                                                                    gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCCGTGTTGTAGATCAATTGCGCACTATCATTTAAATCATCAcccg                                                      >  1:132831/1‑141 (MQ=255)
                                                                                     cctgTTTGTGCGGTACGCAGTGATAATTACCGCGGTTGTCACGCGGATCATAGTCACTTGATGTGACTATGATGCGATTAAGACGCTCTCCGCTACGCAGGGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGc                                                     <  2:213961/140‑1 (MQ=255)
                                                                                          ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGCCTATGATCAGCTCAAAACTGTTTCCGAAGTATAAGCGAGGGGAGCAATTGCGCACTATTATTGAAATAATTACCTGCTCGTg                                                >  2:139402/1‑143 (MQ=255)
                                                                                            gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat                                              <  1:398647/143‑1 (MQ=255)
                                                                                              tCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatta                                            <  1:599961/143‑1 (MQ=255)
                                                                                                gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAACTAATTACCCGCTACTGattatt                                          >  1:241792/1‑143 (MQ=255)
                                                                                                 ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTCTCATTGAAATAATTACCTGCTAGCTAGTAttt                                         >  2:561698/1‑143 (MQ=255)
                                                                                                 ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTCCCTGCTAGTTATGAttt                                         >  2:132999/1‑143 (MQ=255)
                                                                                                    agcagcAATACTTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa                                      <  2:360189/143‑1 (MQ=255)
                                                                                                      cagcaATAATTACCCCGGTGGTCACCCGGATCATAGGCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc                                    <  1:14552/143‑1 (MQ=255)
                                                                                                           ataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGCTTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCCGCTCGTGATTATTTCAACCTACTg                               >  2:417653/1‑143 (MQ=255)
                                                                                                                tACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAAttt                          <  2:565035/143‑1 (MQ=255)
                                                                                                                 gtcccGGTTGTCACCCGTATCATAGTCACTGTATGTGGCTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTc                         <  2:481859/141‑1 (MQ=255)
                                                                                                                           tCACCCGGAGCATAGGCACTTGAGTTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAtttt               <  1:290610/143‑1 (MQ=255)
                                                                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAAttttt              <  1:198679/143‑1 (MQ=255)
                                                                                                                                   aTCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCActc       >  1:498208/1‑143 (MQ=255)
                                                                                                                                   aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATAACACGCTAACGATTATTTCAACCCTCTCAGTTTCATCTACTTTTTTTCActc       >  1:90147/1‑143 (MQ=255)
                                                                                                                                      aTACTCATTTGTTTTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTAt    <  1:587659/143‑1 (MQ=255)
                                                                                                                                        aGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGGGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATgg  >  1:344645/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                        
GCAACGCCGATGGTGTCTGTAACGCTGGTTTCCACACTGGTGGGAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCTACTGAATTTCATCTAATTTTTTTCACTCTATGG  >  NC_000913/3657457‑3657735

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