Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,679,5250GT95.3% 122.5 / ‑6.1 43A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (22/19);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTTG  >  NC_000913_3_pae_tpiA/2679384‑2679670
                                                                                                                                             |                                                                                                                                                 
tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGAtt                                                                                                                                            >  2:410732/1‑149 (MQ=255)
tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCAGCTTCCGGCCAATCCGGCtt                                                                                                                                            >  1:121591/1‑149 (MQ=255)
   ttattaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                                                                                                         <  1:824218/149‑1 (MQ=255)
    tctcatcaacGGCTTCTTTTAGTGTTTTTGTGATGAGCGGCAGCGCGTGCTGCTAAAATATTACAGTTAATAGAATCTAATGGATACGGTTAAAAAAAAAATGCATCACGGTAGTGAATAATGCCGCTTTAGGCCAATCCGGATTAATc                                                                                                                                        <  2:56868/143‑1 (MQ=255)
     attaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGGCAATCCGGATTAATCg                                                                                                                                       >  1:325327/1‑149 (MQ=255)
           aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCa                                                                                                                                 >  2:630012/1‑149 (MQ=255)
           aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAAAAGGGTTTAGTGAACGCGGTAAAAAATAAAAGGCATCAAGGTGGTAAATACTGCTGCTATAAGCCAATCCGCATGAATCGCCGTCa                                                                                                                                 >  2:444707/1‑149 (MQ=255)
              ggCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAAtt                                                                                                                              <  1:157612/149‑1 (MQ=255)
               gCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCCTCACGGCAGTGAATAAGGCAGCTTCAACCCAAGCACGCTTAATCCCCGTCACTTc                                                                                                                             >  2:843880/1‑149 (MQ=255)
                        aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGCATAATGCCGCTTTAGACCAATCCGGATTAATTGCCGGCAATTCAgctccggc                                                                                                                    >  2:401631/1‑143 (MQ=255)
                          ttattTTGTGATGAGCGGCAGCGCCTGCACTGCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                                                                                                                  >  2:668338/1‑149 (MQ=255)
                              tttGTGATGGGCGGCAGCGCGTGCCGTACAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTc                                                                                                              <  2:668885/149‑1 (MQ=255)
                                    atgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTGAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTg                                                                                                        >  2:810186/1‑149 (MQ=255)
                                        gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                                                                                                    >  2:240532/1‑149 (MQ=255)
                                              gcgcGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGCGCAGACTCATATTGCGCAACCTcc                                                                                              >  2:496084/1‑149 (MQ=255)
                                                    gCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGAGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCGCTCAATTCAGCGCAAGACAACTCATATCGACCAAccccccgtatt                                                                                        >  2:417089/1‑144 (MQ=255)
                                                            aaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCggg                                                                                >  1:380501/1‑148 (MQ=255)
                                                            aaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGa                                                                                <  1:810186/149‑1 (MQ=255)
                                                                ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCAGCTTTCGACCAATCCGGATTAATCACCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACATCCCCGAGTCCCCCCGACGcc                                                                            >  2:13244/1‑149 (MQ=255)
                                                                             ggTTTAGGGAACACGGATAAAAATAAAATGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCTTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGccacca                                                               <  2:325327/149‑1 (MQ=255)
                                                                                   gggAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGc                                                         <  1:321200/147‑1 (MQ=255)
                                                                                         gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGACGACAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAttttt                                                   >  1:365099/1‑149 (MQ=255)
                                                                                                 aaaTAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgc                                           <  1:410732/149‑1 (MQ=255)
                                                                                                   aTAAAAGGCAACACGGTGGTGAATAACGCCGCTTTCGGCCAATCCGGATTAATCGCCGCCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTAtgctg                                         <  1:119881/149‑1 (MQ=255)
                                                                                                        aGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgtt                                    <  2:230022/149‑1 (MQ=255)
                                                                                                        aGGCATCACGGTGGTCAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgtt                                    <  1:536477/149‑1 (MQ=255)
                                                                                                          gCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCtgttgt                                  <  1:329591/149‑1 (MQ=255)
                                                                                                             tCACGGTGCTGAATAATGCGGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGTGGAACCTCCCCGAGTCCAGCGGACGCCAGCATCGCCACCATCAGTCAGTCTTCTGTATGCTGCTGTTGTaaa                               >  2:784623/1‑149 (MQ=255)
                                                                                                                cGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCg                            >  1:124448/1‑149 (MQ=255)
                                                                                                                  gtggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGcc                          <  1:630012/149‑1 (MQ=255)
                                                                                                                   tggtgAATAATGCCGTTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAATCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCa                         <  2:68682/149‑1 (MQ=255)
                                                                                                                   tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGa                                                                                          <  1:583333/84‑1 (MQ=255)
                                                                                                                   tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGa                                                                                          >  2:583333/1‑84 (MQ=255)
                                                                                                                        aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATgctgct                                       >  1:696513/1‑130 (MQ=255)
                                                                                                                        aataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATgctgct                                       <  2:696513/130‑1 (MQ=255)
                                                                                                                         ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGc                   <  1:240532/149‑1 (MQ=255)
                                                                                                                               ccGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGacacc             >  2:227208/1‑147 (MQ=255)
                                                                                                                                 gtgtgctgCCAAGCCGTATTAATCGTCGACAAATCAGCGCTAGCAGAATTATATTGCGCAACCTCCAAGAGTCCCTTGTATTCTAGCATCCTCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGTcaacaa           <  2:725273/142‑1 (MQ=255)
                                                                                                                                 gCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATCGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGGCacaacag           >  2:923984/1‑148 (MQ=255)
                                                                                                                                  cTTTCGGCCAATCCGGATTAATCGACCACAACGGCGCTGAAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAg          >  1:225444/1‑149 (MQ=255)
                                                                                                                                   tttCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGt         <  1:527988/149‑1 (MQ=255)
                                                                                                                                     tCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCa       <  1:139060/149‑1 (MQ=255)
                                                                                                                                          cAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCAAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACTGCGGCAACAAGCCAGGTTg  >  1:703511/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGTTG  >  NC_000913_3_pae_tpiA/2679384‑2679670

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