Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,695,997 C→T *63* (TAG→TAA yhjR ← DUF2629 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,695,9970CT98.3% 187.4 / NA 60*63* (TAG→TAAyhjRDUF2629 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/25);  minor base A (0/1);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTGG  >  NC_000913/3695865‑3696132
                                                                                                                                    |                                                                                                                                       
ccGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCATTACttt                                                                                                                                   >  1:1313930/1‑139 (MQ=255)
           tCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaaa                                                                                                                        <  2:1169090/139‑1 (MQ=255)
             gACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAAct                                                                                                                      <  2:481762/139‑1 (MQ=255)
                  ccaGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCTCAAACTCTGcc                                                                                                                 >  1:1090431/1‑139 (MQ=255)
                           tttCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                        <  1:863680/139‑1 (MQ=255)
                             tctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCaaa                                                                                                                        >  2:397415/1‑121 (MQ=255)
                             tctcCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCATACTTTTGTTGCGCaaa                                                                                                                        <  1:397415/121‑1 (MQ=255)
                                         tGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGc                                                                                          >  1:614146/1‑139 (MQ=255)
                                          gTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGcc                                                                                         <  1:1292964/139‑1 (MQ=255)
                                           tAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                        <  1:801138/123‑1 (MQ=255)
                                           tAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                        >  2:801138/1‑123 (MQ=255)
                                               gACCAGGCTAATGCGGCGGTGATGGTTGTAGTCCCCACGCCTCCCGGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGc                                                                                    >  1:930751/1‑139 (MQ=255)
                                                aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCa                                                                                                             <  2:37293/113‑1 (MQ=255)
                                                aCCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCa                                                                                                             >  1:37293/1‑113 (MQ=255)
                                                        aaTGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTc                                                                           <  1:1429907/139‑1 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                            >  1:1219949/1‑136 (MQ=255)
                                                          tGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                            <  2:1219949/136‑1 (MQ=255)
                                                           gcggcgGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTcgcg                                                                        <  2:1208117/139‑1 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                        <  2:1417901/103‑1 (MQ=255)
                                                               cggTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGc                                                                                                        >  1:1417901/1‑103 (MQ=255)
                                                                 gTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                            >  1:1579692/1‑129 (MQ=255)
                                                                 gTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGtt                                                                            <  2:1579692/129‑1 (MQ=255)
                                                                   gATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTTCTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTcgcg                                                                        <  1:592157/131‑1 (MQ=255)
                                                                   gATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTcgcg                                                                        >  2:592157/1‑131 (MQ=255)
                                                                           ttgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCa                                                        <  1:741062/139‑1 (MQ=255)
                                                                            tgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCAt                                                       >  2:772779/1‑139 (MQ=255)
                                                                            tgtCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGACAACTGTTCGCGTTGGGAAATATCGGCAt                                                       >  2:1090999/1‑139 (MQ=255)
                                                                                       cTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGAGAAATATCGGCATAATCAATATCa                                            <  1:1090999/139‑1 (MQ=255)
                                                                                               aCCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGa                                    >  1:1515190/1‑139 (MQ=255)
                                                                                                 ccccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGaaa                                  >  2:1185277/1‑139 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                      >  1:293123/1‑69 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCa                                                                                                      <  2:293123/69‑1 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAAt                                >  2:1247430/1‑139 (MQ=255)
                                                                                                   ccTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAAt                                >  2:357228/1‑139 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTgcg                                                                                                                            >  1:327197/1‑46 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTgcg                                                                                                                            <  2:327197/46‑1 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         <  1:497903/113‑1 (MQ=255)
                                                                                                    cTGCTATCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGc                                                         >  2:497903/1‑113 (MQ=255)
                                                                                                         aTCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGc                          >  2:737603/1‑139 (MQ=255)
                                                                                                         aTCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCACAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGc                          >  1:722550/1‑139 (MQ=255)
                                                                                                          tCCCAGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCt                         >  1:67780/1‑139 (MQ=255)
                                                                                                              aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  1:357228/139‑1 (MQ=255)
                                                                                                              aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  2:614146/139‑1 (MQ=255)
                                                                                                              aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  1:1247430/139‑1 (MQ=255)
                                                                                                              aGTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTaa                     <  1:737603/139‑1 (MQ=255)
                                                                                                               gTACGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAAc                    >  2:818774/1‑139 (MQ=255)
                                                                                                                 aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATg                                                                                           <  2:369267/66‑1 (MQ=255)
                                                                                                                 aCGGCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATg                                                                                           >  1:369267/1‑66 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAata                                                              >  1:1226819/1‑92 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAata                                                              <  2:1226819/92‑1 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCcac               >  2:212258/1‑139 (MQ=255)
                                                                                                                    gCCATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCcac               >  1:1146734/1‑139 (MQ=255)
                                                                                                                      cATTCACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACaa             >  2:277578/1‑139 (MQ=255)
                                                                                                                          cACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAgg                                          >  1:771535/1‑106 (MQ=255)
                                                                                                                          cACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAgg                                          <  2:771535/106‑1 (MQ=255)
                                                                                                                          cACCAATCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATc         >  1:1529873/1‑139 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGcc                                                                                         <  1:147767/54‑1 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGcc                                                                                         >  2:147767/1‑54 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCt    >  2:1347653/1‑139 (MQ=255)
                                                                                                                               aTCCCTTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCt    >  2:381248/1‑139 (MQ=255)
                                                                                                                                 ccctTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTgg  <  1:1347653/139‑1 (MQ=255)
                                                                                                                                 ccctTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTgg  <  1:277578/139‑1 (MQ=255)
                                                                                                                                 ccctTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTgg  <  1:1570082/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CCGGGCAGGCATCGACCACCAGGACATTTTCTCCCAACATTTGTAATGACCAGGCTAATGCGGCGGTGATGGTTGTTGTCCCCACGCCTCCCCGCACCCCCTGCTATCCCAGTACGGCCATTCACCAATCCCCTACTTTTGTTGCGCAAACTCTGCCAGCAACGGCCAGCGTTTTAATGCCGCGGCCAACTGTTCGCGTTGGGAAATATCGGCATAATCAATATCAGGCAGTGAAAATGCCTGCTTTAACGCCACAATATCATTCTGG  >  NC_000913/3695865‑3696132

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