Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,164 C→T 11.1% D453N (GAT→AAT)  cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,1640CT11.1% 98.1 / 3.6 36D453N (GAT→AAT) cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base T (2/2);  total (16/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 = 7.12e-01

AATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACC  >  NC_000913/2123026‑2123291
                                                                                                                                          |                                                                                                                               
aaTTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAAtt                                                                                                                       >  6:149417/1‑149 (MQ=255)
  ttttattttTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAAtttt                                                                                                                     >  6:214043/1‑149 (MQ=255)
           tGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGtt                                                                                                            >  2:211107/1‑149 (MQ=255)
             tCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGATTTTTCCCCGGGttt                                                                                                           <  6:8378/148‑1 (MQ=255)
                ttATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCa                                                                                                       <  4:90579/149‑1 (MQ=255)
                      acacCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCaa                                                                                                  <  5:149417/148‑1 (MQ=255)
                      acacCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt                                                                                                 >  3:146985/1‑149 (MQ=255)
                        acCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt                                                                                                 <  8:227538/147‑1 (MQ=255)
                        acCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt                                                                                                 >  7:227538/1‑147 (MQ=255)
                                  cgaattcattGCATGATGCGACGTTTACACCCTGCCTTAGCGATACGCGAAACTCACCACATAATAATCTTCGAGATAAGTTCCGGAGCGCACTTCAATTAAATCTATCTGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGcccca                                                                                      <  6:123326/145‑2 (MQ=255)
                                   gCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAg                                                                                    >  6:134488/1‑149 (MQ=255)
                                     attcattGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCg                                                                                  <  2:69080/149‑1 (MQ=255)
                                          ttGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAata                                                                             <  4:76352/149‑1 (MQ=255)
                                                 aTGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGATCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAAt                                                                        <  3:16974/147‑1 (MQ=255)
                                                    cGACGTTTACTCCCGTCCTTCGCGATCCGCGACACGCACCACATAATCCTCTTCGAGATAAGAGACGGAGCGCAGTTCAATTAAATCGTGCTTAATTTTCCACGGGTTTTCAATGCAATGCTTCGCCCGAAGCTGTTTATAAATGGAc                                                                    <  4:146985/148‑1 (MQ=255)
                                                     gACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTc                                                                  >  6:211612/1‑149 (MQ=255)
                                                        gTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGtt                                                               <  1:211107/149‑1 (MQ=255)
                                                             cacCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCAc                                                          <  5:211612/149‑1 (MQ=255)
                                                                ccgtccgtAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCAcc                                                         >  7:243802/1‑147 (MQ=255)
                                                                      gtAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAATAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGtt                                                 <  3:49519/149‑1 (MQ=255)
                                                                          cGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGat                                             <  5:134488/149‑1 (MQ=255)
                                                                               cgcgAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCAc                                        >  3:188553/1‑149 (MQ=255)
                                                                                gcgAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCAc                                        >  2:229734/1‑148 (MQ=255)
                                                                                       gCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGaa                                                                               <  2:140740/102‑1 (MQ=255)
                                                                                       gCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGaa                                                                               >  1:140740/1‑102 (MQ=255)
                                                                                             catcatCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt                                                                                                 <  1:17797/78‑1 (MQ=255)
                                                                                             catcatCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAAt                                                                                                 >  2:17797/1‑78 (MQ=255)
                                                                                                    ctctTCGATATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGAAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGGTTGATATCACCATCAATGGTGACTTTTGCCg                   <  5:214043/149‑1 (MQ=255)
                                                                                                      ctTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGtt                 >  3:163036/1‑149 (MQ=255)
                                                                                                           agaTATGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCTGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg                                <  4:29193/129‑1 (MQ=255)
                                                                                                           agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg                                <  4:32346/129‑1 (MQ=255)
                                                                                                           agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg                                >  3:32346/1‑129 (MQ=255)
                                                                                                           agaTAGGAGCCGGAGCGCACTTCAATTAAATTGAGCGGAATTTTCACCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATg                                >  3:29193/1‑129 (MQ=255)
                                                                                                                     cGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCTGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAaccacc  <  8:243802/149‑1 (MQ=255)
                                                                                                                             acatCAATCACATCGCGCGGACCTTTCCCCGGGGATTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGCTTTCACCAAGCAGGTTGATATCACCATCAATGGTGACTTTAGCCGTGCCCGCGa          <  8:181689/130‑1 (MQ=255)
                                                                                                                             aCTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGa          >  7:181689/1‑133 (MQ=255)
                                                                                                                                          |                                                                                                                               
AATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAGCCGGAGCGCACTTCAATTAAATCGAGCGGAATTTTCCCCGGGTTTTCCAGGCAATGCGTCGCCCCCAGCGGAATATAAATGGACTCGTTTTCACCAAGCAGTTTGATATCACCATCAATGGTGACTTTTGCCGTTCCCGCGACAACCACC  >  NC_000913/2123026‑2123291

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