Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,804,180 G→A 7.5% pseudogene (663/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,804,1800GA 7.5% 165.5 / 3.5 53pseudogene (663/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base G (21/28);  new base A (2/2);  total (23/30)
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.35e-01

GTCCCGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATA  >  NC_000913/1804032‑1804328
                                                                                                                                                    |                                                                                                                                                    
gTCCCGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAAtgttg                                                                                                                                                      >  2:391743/1‑149 (MQ=255)
 tCCCGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAAtgttgt                                                                                                                                                     >  1:262392/1‑149 (MQ=255)
  cccGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATgttgtt                                                                                                                                                    <  2:262392/149‑1 (MQ=255)
    cGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTtaa                                                                                                                                                  >  3:83916/1‑149 (MQ=255)
       tGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTtaataa                                                                                                                                               <  1:12139/149‑1 (MQ=255)
       tGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTtaataa                                                                                                                                               <  1:23968/149‑1 (MQ=255)
       tGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTtaataa                                                                                                                                               <  4:83916/149‑1 (MQ=255)
            aGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGACTTCCAGGGGCTATTTTCGCAGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCt                                                                                                                                          >  3:293826/1‑149 (MQ=255)
                cTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAg                                                                                                                                      <  1:6287/149‑1 (MQ=255)
                  gggCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGc                                                                                                                                    <  2:419496/149‑1 (MQ=255)
                         tagataTGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGtctc                                                                                                                             >  4:346923/1‑149 (MQ=255)
                                 ttCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCtt                                                                                                                     >  3:302775/1‑149 (MQ=255)
                                       aCTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACa                                                                                                               <  4:316404/149‑1 (MQ=255)
                                           tgCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACattat                                                                                                           >  1:17492/1‑149 (MQ=255)
                                                    tATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGt                                                                                                  <  4:302775/149‑1 (MQ=255)
                                                     aTTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGtt                                                                                                 <  1:421999/149‑1 (MQ=255)
                                                      ttCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTg                                                                                                >  3:391697/1‑149 (MQ=255)
                                                       tCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGc                                                                                               <  2:158015/149‑1 (MQ=255)
                                                        cTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGcc                                                                                              <  3:346923/149‑1 (MQ=255)
                                                        cTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGAGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGTCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGcc                                                                                              <  2:347252/149‑1 (MQ=255)
                                                             ccTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATAc                                                                                         <  1:358704/149‑1 (MQ=255)
                                                              cTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACt                                                                                        >  3:186904/1‑149 (MQ=255)
                                                                       gtgGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAAc                                                                                                                          >  2:465843/1‑106 (MQ=255)
                                                                       gtgGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAAc                                                                                                                          <  1:465843/106‑1 (MQ=255)
                                                                          ggAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGa                                                                            >  2:72377/1‑149 (MQ=255)
                                                                           gAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGaa                                                                           <  3:207137/149‑1 (MQ=255)
                                                                                 cAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGTCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCAtt                                                                     <  2:439010/149‑1 (MQ=255)
                                                                                 cAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCAtt                                                                     <  4:314139/149‑1 (MQ=255)
                                                                                 cAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCAtt                                                                     <  1:391743/149‑1 (MQ=255)
                                                                                     ggCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGaaa                                                                 >  2:363873/1‑149 (MQ=255)
                                                                                        tATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCaaa                                                                                                      >  3:444290/1‑109 (MQ=255)
                                                                                        tATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCaaa                                                                                                      <  4:444290/109‑1 (MQ=255)
                                                                                        tATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAaataa                                                              >  1:15168/1‑149 (MQ=255)
                                                                                         aTTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAaataat                                                             <  1:312139/149‑1 (MQ=255)
                                                                                              cgcgGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTg                                                        <  1:363873/149‑1 (MQ=255)
                                                                                               gcgGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTgc                                                       <  1:72377/149‑1 (MQ=255)
                                                                                               gcgGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTgc                                                       <  2:210838/149‑1 (MQ=255)
                                                                                               gcgGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCGACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTgc                                                       <  4:271244/149‑1 (MQ=255)
                                                                                                      ttGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGtctc                                                                                                                             <  3:129228/72‑1 (MQ=255)
                                                                                                      ttGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGtctc                                                                                                                             >  4:129228/1‑72 (MQ=255)
                                                                                                           tCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGt                                           >  3:396840/1‑149 (MQ=255)
                                                                                                            cACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGtt                                          <  4:391697/149‑1 (MQ=255)
                                                                                                               gAATGCCATTTGCCCCTATACATCCATCGGCAATGTTATTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGc                                                                                               >  1:164407/1‑93 (MQ=255)
                                                                                                               gAATGCCATTTGCCCCTATACATCCATCGGCAATGTTATTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGc                                                                                               >  1:40733/1‑93 (MQ=255)
                                                                                                               gAATGCCATTTGCCCCTATACATCCATCGGCAATGTTATTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGc                                                                                               <  2:164407/93‑1 (MQ=255)
                                                                                                               gAATGCCATTTGCCCCTATACATCCATCGGCAATGTTATTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGc                                                                                               <  2:40733/93‑1 (MQ=255)
                                                                                                                      aTTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGaaa                                                                 <  3:332206/116‑1 (MQ=255)
                                                                                                                      aTTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGaaa                                                                 >  4:332206/1‑116 (MQ=255)
                                                                                                                      aTTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTg                                >  2:264694/1‑149 (MQ=255)
                                                                                                                                ataCATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTAttttt                      <  4:186904/149‑1 (MQ=255)
                                                                                                                                ataCATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTAttttt                      >  3:310431/1‑149 (MQ=255)
                                                                                                                                    atccatcGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACGCCTGCTTTATTTTTGGCt                  >  3:39761/1‑149 (MQ=255)
                                                                                                                                    atccatcGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCt                  >  4:436703/1‑149 (MQ=255)
                                                                                                                                       catcGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGaaa                                                                 >  1:392042/1‑99 (MQ=255)
                                                                                                                                       catcGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGaaa                                                                 <  2:392042/99‑1 (MQ=255)
                                                                                                                                              aaTGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACg        >  1:459259/1‑149 (MQ=255)
                                                                                                                                                 gttgttAATAATGGCTACAGGCCAGTCTCAACAGCTTATTAGATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGAc     >  3:310850/1‑149 (MQ=255)
                                                                                                                                                    gttAATAAAGGCTAGAGGCCAGTCTCAACAGCTTATTACATTATTCAAACGGTTGCCCATAGTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAAGACCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATa  <  4:293826/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCCCGATGCTCAGACCTGGGCTGATAGATATGTTCGCAACTGTGCAATTGCTATTCTGCCCCTATTACCTGTGGGAGTTCCAGGGGCTATTTTCGCGGGTATTGCATCACGAATGCCATTTGCCCCTATACATCCATCGGCAATGTTGTTAATAATGGCTACAGGCCAGTCTCAACAGCTTATTACATTATTCAAACAGTTGCCCATACTCCCTGAAAAAGAAATCATTGAAATAATAACTGCGCAGAATAGCGTTGGTACACCTGCTTTATTTTTGGCTATGATGAACGGACATA  >  NC_000913/1804032‑1804328

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