Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,329,709 C→G 9.2% G30A (GGT→GCT) ‡ yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,329,7090CG 9.2% 222.7 / 9.3 65G30A (GGT→GCT) ‡yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base G (4/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCT  >  NC_000913/2329575‑2329847
                                                                                                                                      |                                                                                                                                          
gATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGt                                                                                                                              <  1:503198/149‑1 (MQ=255)
      tATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAAc                                                                                                                        >  1:885642/1‑149 (MQ=255)
      tATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTAATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAAc                                                                                                                        >  4:766303/1‑149 (MQ=255)
           tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgcg                                                                                                                                        >  1:1166163/1‑128 (MQ=255)
           tccAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgcg                                                                                                                                        <  2:1166163/128‑1 (MQ=255)
            ccagccagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGCCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGtc                                                                                                                  <  1:1463906/149‑1 (MQ=255)
                 cagGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTAATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcctc                                                                                                             <  3:766303/149‑1 (MQ=255)
                            aTACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAc                                                                                                                                  <  4:968928/117‑1 (MQ=255)
                            aTACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCGCAGCAc                                                                                                                                  >  3:968928/1‑117 (MQ=255)
                            aTACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                  <  2:885642/149‑1 (MQ=255)
                              aCAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGTGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                  <  1:776388/147‑1 (MQ=255)
                              aCAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGTGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTg                                                                                                  >  2:776388/1‑147 (MQ=255)
                                  tttCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTAGCagag                                                                                            >  3:625365/1‑149 (MQ=255)
                                       acaACCTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCagag                                                                                            <  2:1359000/144‑1 (MQ=255)
                                       acaAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCagag                                                                                            >  1:1359000/1‑144 (MQ=255)
                                            cTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGc                                                                                  <  1:1398692/149‑1 (MQ=255)
                                             tCCCTAACGGCGTGGTGAGCGAGTCAGGTAATTCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgc                                                                                                                                         <  2:1179593/93‑1 (MQ=255)
                                             tCCCTAACGGCGTGGTGAGCGAGTCAGGTAATTCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcgc                                                                                                                                         >  1:1179593/1‑93 (MQ=255)
                                              cccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcg                                                                                                                                          >  3:1181202/1‑91 (MQ=255)
                                              cccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACcg                                                                                                                                          <  4:1181202/91‑1 (MQ=255)
                                               ccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTaa                                                                                         >  1:1423865/1‑139 (MQ=255)
                                               ccTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTaa                                                                                         <  2:1423865/139‑1 (MQ=255)
                                                     ggCGTGGTGAGGGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCGGTGCgcagc                                                                                                                                    >  3:1177660/1‑90 (MQ=255)
                                                                    tCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGc                                                                                                                     >  3:1236241/1‑90 (MQ=255)
                                                                    tCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGc                                                                                                                     <  4:1236241/90‑1 (MQ=255)
                                                                    tCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGa                                                                    <  1:267322/139‑1 (MQ=255)
                                                                    tCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGa                                                                    >  2:267322/1‑139 (MQ=255)
                                                                    tCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCa                                                          >  3:84417/1‑149 (MQ=255)
                                                                             gCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcct                                                                                                              >  1:568166/1‑88 (MQ=255)
                                                                             gCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTctcct                                                                                                              <  2:568166/88‑1 (MQ=255)
                                                                              cTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAt                                                >  1:1002966/1‑149 (MQ=255)
                                                                              cTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAt                                                <  2:1002966/149‑1 (MQ=255)
                                                                                gagaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAttt                                              >  1:1450397/1‑149 (MQ=255)
                                                                                 agaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                        <  1:635288/90‑1 (MQ=255)
                                                                                 agaCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                        >  2:635288/1‑90 (MQ=255)
                                                                                   aCTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTAGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                           <  4:625365/149‑1 (MQ=255)
                                                                                     tttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                        <  1:241191/86‑1 (MQ=255)
                                                                                     tttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCaa                                                                                                        >  2:241191/1‑86 (MQ=255)
                                                                                     tttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGAc                                                                      >  3:1009968/1‑120 (MQ=255)
                                                                                     tttGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGAc                                                                      <  4:1009968/120‑1 (MQ=255)
                                                                                                        gacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCa                                            >  4:14082/1‑127 (MQ=255)
                                                                                                        gacAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCa                                            <  3:14082/127‑1 (MQ=255)
                                                                                                         acAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACTGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                           >  2:1111989/1‑127 (MQ=255)
                                                                                                         acAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACTGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                           <  1:1111989/127‑1 (MQ=255)
                                                                                                            gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTc                                                                 <  2:915886/102‑1 (MQ=255)
                                                                                                            gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTc                                                                 >  1:915886/1‑102 (MQ=255)
                                                                                                            gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCgg                                <  4:152398/135‑1 (MQ=255)
                                                                                                            gCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCgg                                >  3:152398/1‑135 (MQ=255)
                                                                                                                     gTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGAt                                                                   >  2:458351/1‑91 (MQ=255)
                                                                                                                     gTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGAt                                                                   <  1:458351/91‑1 (MQ=255)
                                                                                                                      ttGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAacac        <  4:84417/149‑1 (MQ=255)
                                                                                                                        actggTGCTGCGCGGTGCGCAGCACGAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTAtt                                               >  4:1115439/5‑108 (MQ=38)
                                                                                                                         aTCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAgtgt                                                                                                    >  3:1410916/1‑54 (MQ=255)
                                                                                                                         aTCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAgtgt                                                                                                    <  4:1410916/54‑1 (MQ=255)
                                                                                                                            gTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTtgct  <  2:1450397/149‑1 (MQ=255)
                                                                                                                             tgctgcGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGAc                                                                      <  3:35093/80‑1 (MQ=255)
                                                                                                                             tgctgcGCGGTGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGAc                                                                      >  4:35093/1‑80 (MQ=255)
                                                                                                                               ctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                           >  3:1028485/1‑105 (MQ=255)
                                                                                                                               ctgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCAc                                           <  4:1028485/105‑1 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAg                                                         <  1:1090593/90‑1 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAg                                                         >  2:1090593/1‑90 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCt                                         >  2:630729/1‑106 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCt                                         <  1:630729/106‑1 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCaa           <  3:708053/136‑1 (MQ=255)
                                                                                                                                tgcGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCaa           >  4:708053/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                          
GATGGGTATCCTCCAGCCAGGCGTAAACATACAGTTTCCACAAACTCCCTAACGGCGTGGTGAGCGAGTCAGGTAATGCTGAGACTTTGGTCACGCCTGATGACGACAGCTGATACAGTTGATCGTGCTGCGCACCGCGCAGCACCAGTTGTAACGGCGTCTCCTCCGCAAGTGTTGGCAGAGTAACCAACGCCAATAGCCAGACGATTCTGCGCCAGTTCACTTATTTCACCTGCATCCCGGTCCATTCACTCCCTGCCGCAACACTTTGCT  >  NC_000913/2329575‑2329847

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