Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,803,104 A→G 100% pseudogene (11/426 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,803,1040AG95.5% 274.0 / 1.5 88pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (33/50);  total (36/52)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199
                                                                                                 |                                                                                               
taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCg                                                                                               >  2:2305163/1‑100 (MQ=255)
 aaCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGa                                                                                              <  2:1506105/100‑1 (MQ=255)
  aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                             >  1:1579352/1‑100 (MQ=255)
   cAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGata                                                                                            >  1:1441320/1‑100 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCTAACAGATTATGAGTCAAAGCGATATc                                                                                          >  1:2443125/1‑100 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                          <  1:855836/100‑1 (MQ=255)
     aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATc                                                                                          <  1:1261225/100‑1 (MQ=255)
       ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAt                                                                                        <  1:1587734/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  1:376223/100‑1 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      >  2:1865856/1‑100 (MQ=255)
         ttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTa                                                                                      <  2:763229/100‑1 (MQ=255)
          tCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTAt                                                                                     <  2:1606993/100‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    >  2:250869/1‑100 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    <  1:943140/100‑1 (MQ=255)
           cGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATc                                                                                    <  1:611635/100‑1 (MQ=255)
            gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCa                                                                                   <  2:24930/100‑1 (MQ=255)
              gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                 <  2:224799/100‑1 (MQ=255)
                acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATAGAAGGCGAACATATTATGAGTCAAAGCGATATCATTATCAGAAc                                                                               >  2:1585259/1‑100 (MQ=255)
                 caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                              <  2:579883/100‑1 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                         <  2:2267943/100‑1 (MQ=255)
                      cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATta                                                                         <  1:1179614/100‑1 (MQ=255)
                        ttCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtata                                                                       <  2:818384/100‑1 (MQ=255)
                         tCTGGCGATTATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATaa                                                                      <  1:694470/100‑1 (MQ=255)
                          cttacGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     <  1:2753773/96‑1 (MQ=255)
                          ctggttttaattattaattatCATCTTCTAAACCTCCAGTTTTATATATGGCAAAAATATTATGATTCAATTCGATCTCATTATCAGAACTCa                                                                            <  1:1297045/86‑1 (MQ=255)
                          cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAg                                                                     >  2:2842647/1‑100 (MQ=255)
                           tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGt                                                                    >  2:1929266/1‑100 (MQ=255)
                            ngCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                   >  2:2319208/2‑100 (MQ=255)
                               gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAGGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                <  1:2865740/100‑1 (MQ=255)
                                ataataATTAATCATCATCCTCAAAACCTCCAGATAGATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtct                                                                  <  1:2079618/97‑1 (MQ=255)
                                 taataaTTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcctc                                                              >  1:215390/1‑100 (MQ=255)
                                     aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                          <  2:797219/100‑1 (MQ=255)
                                      attaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                         >  2:2687924/1‑100 (MQ=255)
                                      attaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                         >  2:51006/1‑100 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                      <  2:342289/100‑1 (MQ=255)
                                         aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                      <  2:1921111/100‑1 (MQ=255)
                                         aatCATCACCCACAAAAAATCAAGATATATATAAGGAGAACTGATTATGAATAAAAGAGATATCATTATCAGAACTCCTTATCAGTAACCAAATAGAttn                                                      <  1:2612003/100‑2 (MQ=255)
                                          atcatcaacCTCAAGACCTCCAGATATAGATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                     <  1:2336681/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  1:1080829/100‑1 (MQ=255)
                                           tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCAATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                    <  1:918103/100‑1 (MQ=255)
                                            catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                   >  1:2305740/1‑100 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                >  1:1796392/1‑100 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  2:245511/100‑1 (MQ=255)
                                               catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                <  2:1579356/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                               <  1:840188/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                               <  2:1654554/100‑1 (MQ=255)
                                                atcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCg                                               >  1:1279704/1‑100 (MQ=255)
                                                 tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGa                                              <  1:63759/100‑1 (MQ=255)
                                                  ccTCAAAACCTCCAGATATATATAAGGCGACCAGATTATGAGTCAAAGCCATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACAt                                                 >  1:2214998/1‑96 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            >  2:815007/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            >  2:1147848/1‑100 (MQ=255)
                                                   cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                            <  1:905136/100‑1 (MQ=255)
                                                    tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                           <  2:350517/100‑1 (MQ=255)
                                                     cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                          <  2:1878373/100‑1 (MQ=255)
                                                      aaaaCCTGCAGCTATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCAGTATAAGTATCCTCATAGATTGCACATCGATAGCg                                         >  1:1524214/1‑100 (MQ=255)
                                                      aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                         >  2:1076883/1‑100 (MQ=255)
                                                        aaacTCAAGGTATATATAAGGCGAAAAAATTATAAGTCAAAGAGATATTATTATCAGAACTAATTAAAAATATACTCATAGATTGCACATCAATAGCAAc                                       >  1:739897/4‑100 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  1:1237120/100‑1 (MQ=255)
                                                             ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                  <  1:866150/100‑1 (MQ=255)
                                                              cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                 <  2:2177611/100‑1 (MQ=255)
                                                               aGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATTGATAGCGACATACCaa                                >  2:740673/1‑100 (MQ=255)
                                                                gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                               >  2:740900/1‑100 (MQ=255)
                                                                 atatatatAAGCCGAACAGAGTATGAGTCAAAACGATATCATTATCAGAACGCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCCac                               <  1:283287/99‑1 (MQ=255)
                                                                  tatatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAacac                             <  1:1757228/100‑1 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                           >  2:459029/1‑100 (MQ=255)
                                                                    tatataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                           >  2:2589024/1‑100 (MQ=255)
                                                                      tataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                         >  2:1833234/1‑100 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        <  2:1847000/100‑1 (MQ=255)
                                                                       ataAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                        <  1:459028/100‑1 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  2:1700001/1‑100 (MQ=255)
                                                                        taAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                       >  2:2528873/1‑100 (MQ=255)
                                                                         aaGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                      <  2:269714/100‑1 (MQ=255)
                                                                            gCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                   <  2:357683/100‑1 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                <  2:2441596/100‑1 (MQ=255)
                                                                               aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                <  2:2083267/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATATATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  1:591605/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATATCGACATACCAACACCTTCATCAGAGCCTAtt            <  2:475012/100‑1 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            >  1:248141/1‑100 (MQ=255)
                                                                                   gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt            <  1:2250179/100‑1 (MQ=255)
                                                                                     ttATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTan          <  1:2319197/100‑2 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc        <  2:1753237/100‑1 (MQ=255)
                                                                                       aTGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc        <  1:300942/100‑1 (MQ=255)
                                                                                        tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa       >  1:50612/1‑100 (MQ=255)
                                                                                        tGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa       >  2:148376/1‑100 (MQ=255)
                                                                                        ggAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa       >  1:503850/2‑100 (MQ=255)
                                                                                           gTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt    >  2:1484702/1‑100 (MQ=255)
                                                                                             cAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt   >  1:2355026/1‑99 (MQ=255)
                                                                                             cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg  >  2:1220842/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTG  >  NC_000913/1803007‑1803199

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