Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,803,567 G→A pseudogene (50/1476 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,803,5670GA100.0% 101.8 / NA 32pseudogene (50/1476 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/17);  total (15/17)

AACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAAT  >  NC_000913/1803433‑1803695
                                                                                                                                      |                                                                                                                                
aaCCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCa                                                                                                                              >  2:106273/1‑139 (MQ=255)
         cGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTg                                                                                                                     <  1:495968/139‑1 (MQ=255)
                    ggATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTc                                                                                                          <  1:352926/139‑1 (MQ=255)
                          tatTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGc                                                                                                    >  1:249815/1‑139 (MQ=255)
                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAtctc                                                                                                              <  2:527171/124‑1 (MQ=255)
                               aacaggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAAtctc                                                                                                              >  1:527171/1‑124 (MQ=255)
                                   ggaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAATAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTg                                                                                           <  2:34829/139‑1 (MQ=255)
                                    gaacagTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGt                                                                                          <  1:256883/139‑1 (MQ=255)
                                        agTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGcc                                                                                      <  1:106273/139‑1 (MQ=255)
                                         gTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATCGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCACTCTCCTTCTGACGCTATTTTTTGTCGCCa                                                                                     >  1:1020954/1‑139 (MQ=255)
                                               tgctgAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTc                                                                               >  1:617043/1‑139 (MQ=255)
                                                                            tGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGaa                                                  >  2:966673/1‑139 (MQ=255)
                                                                              tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATa                                                                      >  1:1038803/1‑117 (MQ=255)
                                                                              tCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATa                                                                      <  2:1038803/117‑1 (MQ=255)
                                                                                tAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCAtt                                                       <  2:292594/130‑1 (MQ=255)
                                                                                tAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCAtt                                                       >  1:292594/1‑130 (MQ=255)
                                                                                                   aaaaaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCAt                                                                                    >  1:1345876/1‑82 (MQ=255)
                                                                                                   aaaaaCTATACACCCAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCAt                                                                                    <  2:1345876/82‑1 (MQ=255)
                                                                                                               cccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGa                                                                                                       <  2:26233/51‑1 (MQ=255)
                                                                                                               cccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGa                                                                                                       >  1:26233/1‑51 (MQ=255)
                                                                                                                ccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                     >  2:1440745/1‑84 (MQ=255)
                                                                                                                ccAACAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATAc                                                                     <  1:1440745/84‑1 (MQ=255)
                                                                                                                   aCAGTAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTActcct           >  2:180840/1‑139 (MQ=255)
                                                                                                                       tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt       <  1:180840/139‑1 (MQ=255)
                                                                                                                       tAAATCTCAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGt       <  2:617043/139‑1 (MQ=255)
                                                                                                                            ctcAATGGTTATGCATGTTTTTTGCAATCTCCTTGTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAAt  <  1:966673/139‑1 (MQ=255)
                                                                                                                            ctcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAgg                                        >  1:380445/1‑101 (MQ=255)
                                                                                                                            ctcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAgg                                        <  2:380445/101‑1 (MQ=255)
                                                                                                                             tcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTc                                 >  2:623429/1‑107 (MQ=255)
                                                                                                                             tcAATGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTc                                 <  1:623429/107‑1 (MQ=255)
                                                                                                                                 tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGa                          <  2:936784/110‑1 (MQ=255)
                                                                                                                                 tGGTAATGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGACCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGa                          >  1:936784/1‑110 (MQ=255)
                                                                                                                                      |                                                                                                                                
AACCATTTGCGAAGCATACTGGATTATATTGAACAGGAACAGTTAACTGCTGAAGATCGTAACCATTGTATGAAACTGTCTAAAAAAATCCATAGAGAAAAAAACTATACACCCAACAGTAAATCTCAATGGTAGTGCATGTTTTTTGCAATCTCCTTCTGACGCTATTTTTTGTCGCCATCTGTCATTGCAATACGCCCTTGATTCATTGAGAAATGGAAAAGGCAAAGTCAACCTGATTAAACATTACTCCTCCGTTGAAT  >  NC_000913/1803433‑1803695

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