Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,059,464 G→A 6.6% R105R (CGG→CGA)  U712_05205 → Uncharacterized protein yhfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,059,4640GA 6.6% 194.4 / 2.7 61R105R (CGG→CGA) U712_05205Uncharacterized protein yhfF
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

TATGATAAACATGACAATCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCTTTTCATT  >  CP006881/1059327‑1059595
                                                                                                                                         |                                                                                                                                   
tATGATAAACATGACAATCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCgg                                                                                                                                     <  4:82460/138‑1 (MQ=255)
 aTGATAAACATGACACTCCCAGATTCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGa                                                                                                                                    <  7:111412/138‑1 (MQ=255)
     tAAACATGACAATCCAAGATGCAATATAAAAGTAACAGATGGCTTTGAAACAACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGc                                                                                                                                <  7:271372/138‑1 (MQ=255)
         cATGACAATCGAAGATCGAATTTAAAAGTAACAGATGTCTTTGACACCAGCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGTATTGAGGAGTTTCAAGAGGATCATCGGATGGCATggg                                                                                                                           <  5:285376/139‑1 (MQ=255)
         cATGACAATCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATggg                                                                                                                           <  3:8474/139‑1 (MQ=255)
              cAATCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGaa                                                                                                                      <  2:3332/139‑1 (MQ=255)
                 tCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACt                                                                                                                    >  8:210417/1‑138 (MQ=255)
                          cAATATAAAAGTAACCGATGTCTTTGAATCCACATTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGCGGGCAATGGAGACCGAATGGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAg                                                                                                           <  3:112872/138‑1 (MQ=255)
                           aaTATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCg                                                                                                         >  1:99124/1‑139 (MQ=255)
                                      aaCAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGTTTATTCACATGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAg                                                                                              <  8:85139/139‑1 (MQ=255)
                                       aCAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTT‑‑‑GAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTg                                                                                           >  4:33501/1‑138 (MQ=255)
                                       aCAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGGCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAg                                                                                              <  8:140073/138‑1 (MQ=255)
                                           aTGTCTTTGAAACCACCTTTGGTGCGCCTGATTTGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGaa                                                                                         <  7:217265/139‑1 (MQ=255)
                                              tCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGTGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAAt                                                                                        >  3:234170/1‑137 (MQ=255)
                                                              tGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTT‑‑‑GAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTAc                                                                    <  3:33501/138‑1 (MQ=255)
                                                                  gcgcCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACag                                                                  <  2:99124/139‑1 (MQ=255)
                                                                              aGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGc                                                      <  2:293905/139‑1 (MQ=255)
                                                                                   aGTCAATGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGc                                                                                             <  5:42903/95‑1 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAg                                                                                              >  3:11632/1‑94 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAg                                                                                              <  4:11632/94‑1 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGc                                                                                             >  6:42903/1‑95 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGaa                                                                                         >  8:273426/1‑99 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGaa                                                                                         <  7:273426/99‑1 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTAc                                                                    >  8:166699/1‑120 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTAc                                                                    <  7:166699/120‑1 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGt                                                  >  7:210059/1‑138 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGt                                                  >  3:219264/1‑138 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGAATGGc                                                                                                                                >  6:56998/1‑60 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGAATGGc                                                                                                                                >  6:61485/1‑60 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGAATGGc                                                                                                                                <  5:61485/60‑1 (MQ=255)
                                                                                   aGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGAATGGc                                                                                                                                <  5:56998/60‑1 (MQ=255)
                                                                                     tCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTa                                                >  3:153492/1‑138 (MQ=255)
                                                                                              agGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATggg                                                                                                                           >  5:152077/1‑54 (MQ=255)
                                                                                              agGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATggg                                                                                                                           <  6:152077/54‑1 (MQ=255)
                                                                                              agGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaa                                       >  2:22213/1‑138 (MQ=255)
                                                                                              agGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaa                                       >  2:155885/1‑138 (MQ=255)
                                                                                              agGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaa                                       >  8:214435/1‑138 (MQ=255)
                                                                                               gggCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaaa                                      >  1:83084/1‑138 (MQ=255)
                                                                                               gggCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaaa                                      >  6:215472/1‑138 (MQ=255)
                                                                                                ggCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTTAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTaaaaa                                     >  6:120103/1‑138 (MQ=255)
                                                                                                      ggAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGc                              <  4:61045/139‑1 (MQ=255)
                                                                                                      ggAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGGGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGc                              <  4:153335/139‑1 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCgg                                                                                                        >  6:62871/1‑61 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCgg                                                                                                        <  5:62871/61‑1 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTg                                                                                                      <  7:290579/63‑1 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTg                                                                                                      >  8:290579/1‑63 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTg                                                                                                      >  8:286734/1‑63 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTg                                                                                                      <  7:286734/63‑1 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAgg                                                                                    <  4:66487/81‑1 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAgg                                                                                    >  3:66487/1‑81 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCgggg                          >  4:189669/1‑139 (MQ=255)
                                                                                                          aCCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCgggg                          >  8:212989/1‑139 (MQ=255)
                                                                                                               aTTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAg                      >  8:279549/1‑138 (MQ=255)
                                                                                                                ttGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGaa                    <  7:214435/139‑1 (MQ=255)
                                                                                                                             agagGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCttt       >  3:174453/1‑139 (MQ=255)
                                                                                                                                 gATCATCGTTTGGCATGGGCGGAACTTATCAAAAGCGGTGAACTGGAGCTTATTGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAATATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCTTTTCa    <  4:219264/138‑1 (MQ=255)
                                                                                                                                  atcatcGGATGGCATGGGCGGAACTGATCAAAAGGGGTg                                                                                                      <  6:284518/39‑1 (MQ=255)
                                                                                                                                  atcatcGGATGGCATGGGCGGAACTGATCAAAAGCGGTg                                                                                                      >  5:284518/1‑39 (MQ=255)
                                                                                                                                  atcatcGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCTTTTCAtt  <  7:212989/139‑1 (MQ=255)
                                                                                                                                  atcatcGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCCGGGCAAGAACCCCGCAGGCATTTCAt   >  2:68931/1‑138 (MQ=255)
                                                                                                                                  atcatcGGATGGCAGGGGCGGAACTGATCAAACGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCTTTTCAtt  <  7:134203/139‑1 (MQ=255)
                                                                                                                                  atattcGGATGGCAGTGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATAGTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCGTGGGGGGCAAGAACCCCGCAGGCTTTTCAtt  <  8:210059/135‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TATGATAAACATGACAATCCAAGATGCAATATAAAAGTAACAGATGTCTTTGAAACCACCTTTGGTGCGCCTGATATGAGGTTAGTCAAAGGAGAGGGCAATGGAGACCGAATTGAGGAGTTTCAAGAGGATCATCGGATGGCATGGGCGGAACTGATCAAAAGCGGTGAACTGGAGCTGAATGAGGATACTGTTTTGATTACAGAGTTGTTTGTGCTTGTTAAAGATTAAAAAAGCCTGCGGGGCAAGAACCCCGCAGGCTTTTCATT  >  CP006881/1059327‑1059595

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: ‑