Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,460,356 T→C 11.6% pseudogene (708/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,460,3560TC11.6% 2266.3 / 179.4 799pseudogene (708/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (490/195);  new base C (9/81);  total (519/280)
Fisher's exact test for biased strand distribution p-value = 2.82e-30
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 831 total aligned reads displayed.

ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422
                                                                   |                                                                  
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:11370647/1‑68 (MQ=255)
aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                    >  1:21253840/1‑68 (MQ=255)
 cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTa                                                                   >  1:5132803/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:88581/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:25451850/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:22163532/1‑68 (MQ=255)
  ggTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTaa                                                                  >  1:21684931/1‑68 (MQ=255)
    tGAGCTTGTTATTACGGGTGATCATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:7501814/68‑3 (MQ=21)
     gAGCTTGTTATTACGGGTGATAATTCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                <  1:36854318/67‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAg                                                                <  1:29101968/67‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               <  1:14287276/68‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               <  1:4870126/68‑1 (MQ=255)
     gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCaagg                                                                <  1:38526496/67‑3 (MQ=255)
     gAGCTTGTTATTACGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                               <  1:6274584/68‑1 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAt                                                              >  1:598059/1‑67 (MQ=255)
       gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAATGAt                                                              >  1:16054215/1‑67 (MQ=255)
        cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                             >  1:37665845/1‑67 (MQ=255)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGAtt                                                             <  1:16983149/66‑1 (MQ=255)
         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTcc                                                           >  1:16571991/1‑68 (MQ=255)
          tgatattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCa                                                          <  1:9226216/65‑1 (MQ=18)
            ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        >  1:28404993/1‑68 (MQ=255)
             tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCAcc                                                        <  1:33998342/67‑1 (MQ=255)
              attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCgg                                                      >  1:35607961/1‑68 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:36080003/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:25207780/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:11743961/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:13902518/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:34753086/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:34462456/68‑1 (MQ=255)
               ttaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGt                                                     <  1:25714183/68‑1 (MQ=255)
               ttaCGGGTGAAAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGt                                                     <  1:12061203/68‑1 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:24431473/1‑68 (MQ=255)
                 aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAc                                                   >  1:33059063/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:4049268/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:31696021/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                  >  1:13530470/1‑68 (MQ=255)
                  cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACg                                                  <  1:34928597/68‑1 (MQ=255)
                   gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACgg                                                 >  1:19977166/1‑68 (MQ=255)
                   gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                 >  1:6620383/1‑68 (MQ=255)
                   gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACgg                                                 >  1:19151350/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGaa                                               >  1:23833218/1‑68 (MQ=255)
                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                <  1:4493309/67‑1 (MQ=255)
                        aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                            >  1:29933058/1‑68 (MQ=255)
                        aTAATGCCACAGTAAACAACAAAGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAAtt                                            >  1:30866537/1‑67 (MQ=17)
                          aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTGCGGAAATCaa                                          <  1:26456127/68‑1 (MQ=255)
                          aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCa                                           >  1:35924533/1‑67 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:10579620/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:18486318/1‑68 (MQ=255)
                           aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAc                                         >  1:14628706/1‑68 (MQ=255)
                             gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACgg                                       >  1:5670685/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:10006793/1‑68 (MQ=255)
                              ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      >  1:7064710/1‑68 (MQ=255)
                               cacaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:25970859/67‑1 (MQ=255)
                               cacaGTAAAAAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGt                                      <  1:22960593/67‑1 (MQ=255)
                                acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                    <  1:15431523/68‑1 (MQ=255)
                                 caGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                      <  1:20785429/49‑1 (MQ=11)
                                  aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:9431563/1‑68 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:15147683/1‑67 (MQ=255)
                                   gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  >  1:2722436/1‑67 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaata                                  <  1:11052906/65‑1 (MQ=255)
                                     aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACg                               >  1:22611941/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:24073540/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:24666451/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggg                             >  1:723932/1‑68 (MQ=255)
                                       acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACggc                             >  1:30546827/1‑67 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGa                            >  1:20185505/1‑68 (MQ=255)
                                        caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:19646065/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATGAACGGTAATAACGGGaaa                          >  1:37960930/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCGAGGGAAATAACGGGaaa                          >  1:3886703/1‑68 (MQ=17)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCACGGGAATAACGGGatt                          >  1:33656411/1‑66 (MQ=17)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          <  1:5839660/68‑1 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:17710150/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:38210006/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:27704410/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:3536326/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                          >  1:29695697/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGa                            >  1:23374409/1‑66 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATAAAGGGTAATAACGGGaaa                          >  1:18095634/1‑68 (MQ=255)
                                          acaaTGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTATAAACGGGaaa                          >  1:32486532/1‑68 (MQ=14)
                                            aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                        >  1:13187134/1‑68 (MQ=255)
                                            aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGt                        >  1:37448614/1‑68 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     >  1:3267158/1‑68 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:8187496/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     <  1:124576/68‑1 (MQ=255)
                                               ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                     >  1:37093995/1‑68 (MQ=255)
                                                  aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTc                   <  1:28833550/67‑1 (MQ=255)
                                                    aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                <  1:28115692/68‑1 (MQ=255)
                                                    aaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                <  1:17162807/68‑1 (MQ=255)
                                                     accaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGa               <  1:3034913/68‑1 (MQ=255)
                                                      ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc              >  1:26195660/1‑68 (MQ=255)
                                                        accGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg             >  1:31892428/1‑67 (MQ=255)
                                                          cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAGGGAAATAACGGGAAAGTGATTCAGGACGGCg          >  1:31192126/1‑68 (MQ=255)
                                                           gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGa         <  1:32178014/68‑1 (MQ=255)
                                                             tGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTAATTCAGGACGGCGAt        <  1:23262931/67‑1 (MQ=255)
                                                               aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTg     >  1:18640758/1‑68 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    <  1:14094527/68‑1 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:15010903/1‑68 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    >  1:33414646/1‑68 (MQ=255)
                                                                cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg    <  1:9678395/68‑1 (MQ=255)
                                                                  gTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGGAAGTGATTCAGGACGGCGATCTGGAt  >  1:26752885/1‑68 (MQ=255)
                                                                   |                                                                  
ACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAT  >  CP009273/1460289‑1460422

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: