Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 11.8% intergenic (+42/‑206) ygaM → / → nrdH DUF883 domain‑containing protein YgaM/glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5170TA11.8% 114.3 / 4.5 48intergenic (+42/‑206)ygaM/nrdHDUF883 domain‑containing protein YgaM/glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (24/18);  new base A (2/4);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCA  >  NC_000913/2800379‑2800659
                                                                                                                                          |                                                                                                                                              
ttactCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCt                                                                                                                                            <  1:415044/139‑1 (MQ=255)
     cTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                                                                       <  1:691141/143‑1 (MQ=255)
      ttttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCgggc                                                                                                                                      >  2:494567/1‑142 (MQ=255)
         ttttggTGAAAGGCCCTGGTGGAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTAACCCCGAAGACAGGCTTCGGGGtca                                                                                                                                   <  2:537843/138‑3 (MQ=255)
           ttCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCt                                                                                                                                 >  1:386130/1‑143 (MQ=255)
           ttCGTGAAAGA‑CCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAACTAGTAGTGCAAAAATAATAATAAATACGCGTCTTTGACCCTGAAGCCTGTCTACTGGGATTCtt                                                                                                                                >  2:146729/1‑143 (MQ=255)
            tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAACTGATAATAAATACGCGCATGTGAACCCGACGCGTGGCTTCGGGGTTTGtt                                                                                                                                >  1:270064/1‑143 (MQ=255)
             cGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCttt                                                                                                                               >  1:313758/1‑143 (MQ=255)
               tGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATAAGCGCACTGTTAAGCATGCGCCAATCTATTTACAAAAATGATAATAAATACGAGTATTTGACACCGAAGCAAGTCATCGAGGTTTAttttt                                                                                                                             >  2:226697/1‑143 (MQ=255)
                 aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGc                                                                                                                           <  1:655642/143‑1 (MQ=255)
                  aaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTCACCCCGAAGCCTGTCTTCGGCGTTTCTTTTCGcc                                                                                                                          >  2:343612/1‑143 (MQ=255)
                   aGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATGATAAATACGCGTCTTTGACCACGAAGCCTGTTGTCGGGGTTTCTTTTTGCCt                                                                                                                         <  2:329268/143‑1 (MQ=255)
                          ggTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAAt                                                                                                                  <  2:614103/143‑1 (MQ=255)
                            tgtAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                                                                                >  1:441399/1‑143 (MQ=255)
                                   tGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCCAGTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGCCCCCGAAGCCTGTCTTCGGCGTTTCCTTTTCCCCGGTGAATcaccaccct                                                                                                         >  2:236745/1‑137 (MQ=255)
                                      gTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                                                                      >  1:155560/1‑143 (MQ=255)
                                          agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTa                                                                                                  <  1:598712/143‑1 (MQ=255)
                                                  cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATCATCACTGACAAACGCAGACCCCGGCGCCACACTGCACTGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                          <  1:213140/143‑1 (MQ=255)
                                                    gtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAAGCGGAGGGCGAAAATGATAATAAATACGCGTCTATTACCCCAAATCCTGTCTTCGTGGTTTCTTTATGCCTGGTGAATCACAAACATCCCCCTACCCCGacacc                                                                                        >  2:155711/1‑142 (MQ=255)
                                                      gggtaTTTGTATCGGCGCACTGTTAAGCATGGGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCTCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                                                                      <  2:414764/142‑1 (MQ=255)
                                                        gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCCCa                                                                                    >  2:817313/1‑143 (MQ=255)
                                                        gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGCGAATCACACAAACCCGCCAACCCCGTCCCGCTCa                                                                                    >  2:358016/1‑143 (MQ=255)
                                                         taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAGTCACAAAACTCCCCCTCCCCCGGCACGCTCat                                                                                   >  2:534873/1‑143 (MQ=255)
                                                          aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCCACCCCGCCACGccctca                                                                                  >  2:639883/1‑138 (MQ=255)
                                                          aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCGTTGCCCCCGCAGCCTATCTTCGGGGGTTCTTTTTGCCTGCTGACTCACAAAGATACCCCTCCCCCGTCACgactcta                                                                                  >  2:621164/1‑137 (MQ=255)
                                                              ttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                                                              >  2:177816/1‑143 (MQ=255)
                                                               tATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                             <  2:437272/143‑1 (MQ=255)
                                                                      gcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTGTTCGGGGTTTCTTTTTGCCTGGTGAATCACACAAATCCCCCTACCCCGTCACGCTCATATCCACGGTAAtc                                                                      >  2:483556/1‑142 (MQ=255)
                                                                      gcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGAGAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGACTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGCGTGAtt                                                                      >  1:149737/1‑143 (MQ=255)
                                                                             tAAGCATGCGCAAATCGTCGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGCTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATCTCCCGCGCAATTTCGcccc                                                               >  2:427297/1‑142 (MQ=255)
                                                                                gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCCCTa                                                            >  2:224467/1‑143 (MQ=255)
                                                                                 cATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                           >  1:756159/1‑143 (MQ=255)
                                                                                     cgaaaaTCGTAGTGCAAAAACGATAATAAATAAGCGTGTTTGACCCCGCCTCCTGTGAATCCAGTTACCTTCGGACTGGTGAATCACAAAAATCCCCCTACCAAACAACGCACATATCCAGGGTAATTTCGACCACTATTTGc                                                       <  2:696953/140‑1 (MQ=255)
                                                                                      gcAAATCGGAGTACAAAAATGATAATAAAGATGGGTCTTTTACCCCGAAGCGTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                      <  2:48921/143‑1 (MQ=255)
                                                                                       cAAATCGTTGTGCTAAGAGGATTATAAATACGCGTCTTTGATCCCGAAGCCAGTCTTCGGGGTGTCTCTTGCCCTGCTGAATCACAAAAATCCCCCTACCACGCCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCta                                                     <  2:474216/143‑1 (MQ=255)
                                                                                                 tGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                           >  1:211611/1‑143 (MQ=255)
                                                                                                    aaaaaTGATAATAAATATGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGt                                        <  2:625942/143‑1 (MQ=255)
                                                                                                     aaacTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGCGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                       <  1:77154/143‑1 (MQ=255)
                                                                                                     aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                       <  1:514676/143‑1 (MQ=255)
                                                                                                             aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAGCACTATTTGCTATATATTGTGTGGTTGAATCttt                               <  1:847947/143‑1 (MQ=255)
                                                                                                             aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                               <  2:313758/143‑1 (MQ=255)
                                                                                                              ataaataCGCGTCTTTGACCCCAAAGCCTGTCTTCCGGTCTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTTAATCtttt                              <  2:331264/143‑1 (MQ=255)
                                                                                                               taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTGTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTGGAATCttttt                             >  1:310124/1‑143 (MQ=255)
                                                                                                                    aCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATGTTTGGTGGTTGAATCTTTTTTCAAc                        >  2:902713/1‑143 (MQ=255)
                                                                                                                       cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTAc                     <  2:211611/143‑1 (MQ=255)
                                                                                                                            ttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTa                <  1:500424/143‑1 (MQ=255)
                                                                                                                               aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACCTCTAGTa             >  1:658583/1‑143 (MQ=255)
                                                                                                                                      aGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt      <  1:494567/143‑1 (MQ=255)
                                                                                                                                          tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGAGTAATTTCGACCACTATTTGCTATCTATTGTATGGTTGAATCTTTTTTCCACTACATCTAGTATCTCTGTCTCa  >  2:477495/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCA  >  NC_000913/2800379‑2800659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑