Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,389,752 A→G N213D (AAT→GAT)  yfbP → TPR‑like repeats‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,389,7520AG97.3% 135.9 / NA 37N213D (AAT→GAT) yfbPTPR‑like repeats‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/16);  minor base T (0/1);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATAATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGATTTTAATATCGAT  >  NC_000913/2389605‑2389895
                                                                                                                                                   |                                                                                                                                               
aaTTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGAtt                                                                                                                                              >  2:181014/1‑151 (MQ=255)
        aTTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTg                                                                                                                                      >  2:210037/1‑151 (MQ=255)
           tGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGAt                                                                                                                                    <  1:204808/150‑1 (MQ=255)
                     gACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGAtgtgtg                                                                                                                         >  1:583021/1‑151 (MQ=255)
                     gACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGAtgtgtg                                                                                                                         <  1:56636/151‑1 (MQ=255)
                     gACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGAtgtgt                                                                                                                          <  1:572712/150‑1 (MQ=255)
                           aaaGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATa                                                                                                                   <  1:320131/151‑1 (MQ=255)
                                   tcattTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTAttt                                                                                                            <  2:407239/150‑1 (MQ=255)
                                       ttGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTgcgc                                                                                                       >  1:142399/1‑151 (MQ=255)
                                           tGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCag                                                                                                     >  2:581244/1‑149 (MQ=255)
                                            gTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCag                                                                                                     >  2:227146/1‑148 (MQ=255)
                                             tCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAgaga                                                                                                  >  2:306218/1‑150 (MQ=255)
                                                     cGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCtt                                                                                         <  1:537226/151‑1 (MQ=255)
                                                            gTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCAt                                                                                  <  1:284347/151‑1 (MQ=255)
                                                                  aaTTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAAccc                                                                            >  2:234495/1‑151 (MQ=255)
                                                                   aTTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCt                                                                           <  1:317055/151‑1 (MQ=255)
                                                                     ttCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTg                                                                          <  1:344453/150‑1 (MQ=255)
                                                                         aGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTAc                                                                      >  1:339607/1‑150 (MQ=255)
                                                                            atgatCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATaa                                                                  <  2:100485/151‑1 (MQ=255)
                                                                            atgatCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATa                                                                   >  2:591999/1‑150 (MQ=255)
                                                                                   tgtgATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGcc                                                           >  1:9084/1‑151 (MQ=255)
                                                                                       aTTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATaa                                                                  >  1:555682/1‑140 (MQ=255)
                                                                                       aTTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATaa                                                                  <  2:555687/140‑1 (MQ=255)
                                                                                             tGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAAGAACAAAGGAAAATAATTAGTACCCCTTGTTGCAGATGGGTGTTGATACTTATTTAGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATa                                                     <  1:51002/147‑1 (MQ=255)
                                                                                                   aaaGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTg                                           >  2:342093/1‑151 (MQ=255)
                                                                                                     aGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGAGGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGa                                          >  2:473508/1‑150 (MQ=255)
                                                                                                            ggATTAGATCTTAGCTTTGGAGGAACAAAGGAACATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTg                                  >  2:618598/1‑151 (MQ=255)
                                                                                                             gATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGc                                 >  2:553821/1‑151 (MQ=255)
                                                                                                              aTTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCa                                <  1:234495/151‑1 (MQ=255)
                                                                                                                 aGATCTTAGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCAttt                             <  1:280870/151‑1 (MQ=255)
                                                                                                                        aGCTTTGGAGGAACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTTTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGaa                       <  1:329117/150‑1 (MQ=255)
                                                                                                                              ggaggaACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATg                <  2:486568/151‑1 (MQ=255)
                                                                                                                              ggaggaACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATg                >  2:639825/1‑151 (MQ=255)
                                                                                                                                aggaACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGAt              >  1:490586/1‑151 (MQ=255)
                                                                                                                                  gaACAAAGGAAAATAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGAttt            <  2:281773/151‑1 (MQ=255)
                                                                                                                                           aaaaTAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGATTTTAATATCg    >  2:454504/1‑150 (MQ=255)
                                                                                                                                            aaaTAATGATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGATTTTAATATCGAt  >  1:278557/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AATTGGGCATTTGGTGTTCCTGACAATAAAGTCATTCATTTGCTGTCACAGGCCGCAGATGTCCCGAATTTCCAGGATGATCCTGTGATTCAGTGCATTAAAGAACAAGGATTAGATCTTAGCTTTGGAGGAACAAAGGAAAATAATAATTACCCCTTGATGCAGATGTGTGTTGATACTTATTTTGCGCAGAGAAAACCGCTTCAGGCATTAACCCTGTTACATAATTATGCCTGGATAATGTCCAGTGAAACGACTGCATTTCAGGAACGCTATGATTTTAATATCGAT  >  NC_000913/2389605‑2389895

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