Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,724,478 (G)6→7 coding (178/522 nt) sodC ← superoxide dismutase (Cu‑Zn)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,724,4721.G97.6% 137.2 / ‑3.2 41coding (184/522 nt)sodCsuperoxide dismutase (Cu‑Zn)
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (19/21);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTG  >  NC_000913/1724333‑1724604
                                                                                                                                            |                                                                                                                                    
ccGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCAcc                                                                                                                                       >  1:961942/1‑140 (MQ=255)
 cGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑g                                                                                                                                     <  2:996804/140‑1 (MQ=255)
 cGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑g                                                                                                                                     >  2:1363358/1‑140 (MQ=255)
 cGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑g                                                                                                                                     <  2:1043626/140‑1 (MQ=255)
 cGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑g                                                                                                                                     <  1:17373/140‑1 (MQ=255)
    cACCTTCTGGCCTTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggg                                                                                                                                  <  2:549636/140‑1 (MQ=255)
      ccTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggg                                                                                                                                <  1:1363358/140‑1 (MQ=255)
      ccTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCACTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggg                                                                                                                                <  2:1367447/140‑1 (MQ=255)
        ttCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggggt                                                                                                                              <  2:695944/140‑3 (MQ=255)
         tCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggggta                                                                                                                             <  2:785749/140‑4 (MQ=255)
          cTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggggtaa                                                                                                                            <  2:491271/140‑5 (MQ=255)
          cTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggggtaa                                                                                                                            <  2:1211375/140‑5 (MQ=255)
                       tttACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGGAATGCTTTCAGATCgg                                                                                                               >  1:112685/1‑140 (MQ=255)
                            cGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGa                                                                                                           <  2:1454099/140‑1 (MQ=255)
                             ggTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGaa                                                                                                          <  2:961942/140‑1 (MQ=255)
                                  tttGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTc                                                                                                     >  1:1168962/1‑140 (MQ=255)
                                    tgtgGATCAAGATGCCCGCCTGCGGATTCCGCGGCGTTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCa                                                                                                   <  1:587615/140‑1 (MQ=255)
                                    tgtgGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCa                                                                                                   <  2:485070/140‑1 (MQ=255)
                                    tgtgGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCa                                                                                                   >  1:858205/1‑140 (MQ=255)
                                    tgtgGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCa                                                                                                   <  1:882514/140‑1 (MQ=255)
                                    tgtgGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCa                                                                                                   <  1:1140763/140‑1 (MQ=255)
                                             aGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTAt                                                                                          <  1:429988/140‑1 (MQ=255)
                                                 gCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                      <  2:858205/140‑1 (MQ=255)
                                                   ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                                                                                           <  2:928923/133‑1 (MQ=255)
                                                   ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTa                                                                                           >  1:928923/1‑133 (MQ=255)
                                                             aTTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCCGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGt                                                                          <  2:1256226/140‑1 (MQ=255)
                                                                 cgcgGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                      <  2:423350/124‑1 (MQ=255)
                                                                 cgcgGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                      >  1:423350/1‑124 (MQ=255)
                                                                  gcggcgCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGc                                                                     <  1:1340612/140‑1 (MQ=255)
                                                                            cTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                           <  2:1168962/140‑1 (MQ=255)
                                                                                gCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                           <  2:1180370/136‑1 (MQ=255)
                                                                                gCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                           >  1:1180370/1‑136 (MQ=255)
                                                                                 ccATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGTAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGcc                                                      <  2:877089/140‑1 (MQ=255)
                                                                                   aTCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTa                                                   <  2:112685/140‑1 (MQ=255)
                                                                                         ggtggCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                           <  1:329860/127‑1 (MQ=255)
                                                                                         ggtggCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                           >  2:329860/1‑127 (MQ=255)
                                                                                                  ctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTAcccc                                               >  2:851183/1‑130 (MQ=255)
                                                                                                  ctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTAcccc                                               <  1:851183/130‑1 (MQ=255)
                                                                                                  ctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTg                                     >  1:1160342/1‑140 (MQ=255)
                                                                                                            ccTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCa                           >  1:1037571/1‑140 (MQ=255)
                                                                                                                   cATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGAc                    >  1:1316627/1‑140 (MQ=255)
                                                                                                                      gAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACttt                 >  1:9279/1‑140 (MQ=255)
                                                                                                                       aaTATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACtttt                >  1:1232381/1‑140 (MQ=255)
                                                                                                                        atatGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACttttt               <  2:1232381/140‑1 (MQ=255)
                                                                                                                              aaGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCa             >  2:1461953/1‑136 (MQ=255)
                                                                                                                              aaGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCa             <  1:1461953/136‑1 (MQ=255)
                                                                                                                                gCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCa       >  1:1214493/1‑140 (MQ=255)
                                                                                                                                 ccATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAg      >  2:417651/1‑140 (MQ=255)
                                                                                                                                 ccATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAg      >  1:878228/1‑140 (MQ=255)
                                                                                                                                  cATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGc     <  2:1037571/140‑1 (MQ=255)
                                                                                                                                   aTGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCt    >  2:1373562/1‑140 (MQ=255)
                                                                                                                                    tGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCtt   >  1:1191403/1‑140 (MQ=255)
                                                                                                                                     gTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTtg  >  2:419257/1‑140 (MQ=255)
                                                                                                                                     gTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTtg  <  2:9279/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                    
CCGGCACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTG  >  NC_000913/1724333‑1724604

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