Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,633 G→T 10.8% G296G (GGC→GGA) ‡ cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6330GT10.8% 117.3 / 5.6 37G296G (GGC→GGA) ‡cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base G (13/20);  new base T (2/2);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTT  >  NC_000913/2123485‑2123781
                                                                                                                                                    |                                                                                                                                                    
gTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTg                                                                                                                                                      <  3:429162/149‑1 (MQ=255)
                    aTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                  <  4:118159/149‑1 (MQ=255)
                     tCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTg                                                                                                                                 >  3:290387/1‑149 (MQ=255)
                        gCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtc                                                                                                                              <  1:114979/149‑1 (MQ=255)
                          aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtctc                                                                                                                            <  1:403465/149‑1 (MQ=255)
                              cATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCAt                                                                                                                        <  4:290387/149‑1 (MQ=255)
                                    tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                  <  4:93352/149‑1 (MQ=255)
                                    tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                  <  3:97429/149‑1 (MQ=255)
                                    tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                  <  1:295161/149‑1 (MQ=255)
                                    tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAAAGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                  <  4:276617/149‑1 (MQ=255)
                                       tCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAc                                                                                                               <  2:26262/149‑1 (MQ=255)
                                                 tACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCaa                                                                                                     <  4:39598/149‑1 (MQ=255)
                                                 tACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCAcggcgat                                                                                                                                           <  4:442261/111‑3 (MQ=37)
                                                   ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                          >  3:245586/1‑126 (MQ=255)
                                                   ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                          <  4:245586/126‑1 (MQ=255)
                                                   ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACa                                                                                                   >  3:282364/1‑149 (MQ=255)
                                                   ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCAcggcgat                                                                                                                                           >  4:46427/1‑107 (MQ=37)
                                                    cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACAt                                                                                                  >  2:219764/1‑149 (MQ=255)
                                                    cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACAt                                                                                                  <  1:190118/149‑1 (MQ=255)
                                                     aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCa                                                                                                      <  2:142385/144‑1 (MQ=255)
                                                     aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCa                                                                                                      >  1:142385/1‑144 (MQ=255)
                                                       aTCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGc                                                                                               >  3:319480/1‑149 (MQ=255)
                                                         cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             <  3:184544/149‑1 (MQ=255)
                                                                 ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaa                                                                                                                                   >  4:78586/1‑98 (MQ=21)
                                                                 ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaa                                                                                                                                   <  3:78586/103‑6 (MQ=18)
                                                                                    cAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                  <  1:219764/149‑1 (MQ=255)
                                                                                        ttCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCAt                                                              >  1:370387/1‑149 (MQ=255)
                                                                                              aTACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             >  3:348320/1‑112 (MQ=255)
                                                                                              aTACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                             <  4:348320/112‑1 (MQ=255)
                                                                                                  acaTAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACg                                                    >  3:363979/1‑149 (MQ=255)
                                                                                                           ttttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                  >  2:467121/1‑126 (MQ=255)
                                                                                                           ttttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                  <  1:467121/126‑1 (MQ=255)
                                                                                                                aGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGt                                      <  2:187178/149‑1 (MQ=255)
                                                                                                                       tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                          <  3:453187/138‑1 (MQ=255)
                                                                                                                       tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                          >  4:453187/1‑138 (MQ=255)
                                                                                                                       tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                               >  2:35118/1‑149 (MQ=255)
                                                                                                                          ttAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCGAGCCAGCATCCATCGGCACCACAACAGCATCTGCCGTAAGTTCCATGACCGCGTAATACACCGa                            >  4:240142/1‑149 (MQ=255)
                                                                                                                                         tGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGc             <  1:35118/149‑1 (MQ=255)
                                                                                                                                                    gCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGctt  <  2:370387/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTT  >  NC_000913/2123485‑2123781

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