Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,091,145 C→T G136G (GGC→GGT)  sodA → superoxide dismutase [Mn]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,091,1450CT97.2% 127.6 / ‑3.2 36G136G (GGC→GGT) sodAsuperoxide dismutase [Mn]
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/18);  minor base A (0/1);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

GGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGCGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCG  >  NZ_CP009273/4091011‑4091279
                                                                                                                                      |                                                                                                                                      
ggCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGAt                                                                                                                                     >  2:656099/1‑138 (MQ=255)
        ccTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            <  2:594514/139‑1 (MQ=255)
        ccTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            <  1:395383/139‑1 (MQ=255)
               ggTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCggtg                                                                                                                         <  2:654526/135‑1 (MQ=255)
               ggTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCggtg                                                                                                                         >  1:654526/1‑135 (MQ=255)
                 tGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            >  1:313373/1‑130 (MQ=255)
                 tGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            <  2:313373/130‑1 (MQ=255)
                 tGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTc                                                                                                                    >  2:525892/1‑138 (MQ=255)
                 tGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTc                                                                                                                    >  2:277597/1‑138 (MQ=255)
                  gACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTa                                                                                                                  >  1:372365/1‑139 (MQ=255)
                    ccTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGGTtcttct                                                                                                                >  2:36216/1‑139 (MQ=255)
                          aGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAAc                                                                                                          <  2:715298/139‑1 (MQ=255)
                             ggCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAg                                                                                                       >  1:531322/1‑139 (MQ=255)
                               cTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTg                                                                                                               <  1:535959/129‑1 (MQ=255)
                               cTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTg                                                                                                               >  2:535959/1‑129 (MQ=255)
                                         tGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTg                                                                                           <  2:239118/139‑1 (MQ=255)
                                                ggCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTg                                                                                    <  2:229331/139‑1 (MQ=255)
                                                       ttGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATtct                                                                                                 <  1:325892/119‑1 (MQ=255)
                                                       ttGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATtct                                                                                                 >  2:325892/1‑119 (MQ=255)
                                                               ttCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTTTTCTGCTAACCAGGATTCTCCGCTGATGGGTGAAGGTATTTCTGgcg                                                                     >  2:745272/1‑139 (MQ=255)
                                                                 cacaGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            >  2:347175/4‑82 (MQ=255)
                                                                 cAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCg                                                                                                                            <  1:347175/82‑1 (MQ=255)
                                                                    aGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGCTTTTACCTTTAACCAGGATACTCCTCTGGTGGGGGGAGCTATTTCTGGCGCTTcc                                                                >  1:379319/1‑139 (MQ=255)
                                                                          aTTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTc                                                          <  1:392646/139‑1 (MQ=255)
                                                                             tGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGAGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTccc                                                        <  1:282856/138‑1 (MQ=255)
                                                                                 aaaGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTa                                                   <  1:656099/139‑1 (MQ=255)
                                                                                                     tGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTggg                                 <  2:594885/137‑1 (MQ=255)
                                                                                                      ggTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGaa                               <  1:743228/138‑1 (MQ=255)
                                                                                                      ggTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGaa                               >  1:107756/1‑138 (MQ=255)
                                                                                                       gTTCCGGCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACa                             <  1:749960/139‑1 (MQ=255)
                                                                                                             gCTGGGCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTta                       <  2:107756/139‑1 (MQ=255)
                                                                                                                 ggCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACa                             >  2:513281/1‑129 (MQ=255)
                                                                                                                 ggCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACa                             <  1:513281/129‑1 (MQ=255)
                                                                                                                  gCATGGCTGGTGCTGAAAGGTGATAAACTGGCGGTGGGTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTAcc                  >  2:156342/1‑139 (MQ=255)
                                                                                                                                gAAAGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAAc    >  2:439405/1‑139 (MQ=255)
                                                                                                                                  aaGGTGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAAccg  <  1:439405/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GGCACCACCCTGCAGGGTGACCTGAAAGCGGCTATCGAACGTGACTTCGGCTCCGTTGATAACTTCAAAGCAGAATTTGAAAAAGCGGCAGCTTCCCGCTTTGGTTCCGGCTGGGCATGGCTGGTGCTGAAAGGCGATAAACTGGCGGTGGTTTCTACTGCTAACCAGGATTCTCCGCTGATGGGTGAAGCTATTTCTGGCGCTTCCGGCTTCCCGATTATGGGCCTGGATGTGTGGGAACATGCTTACTACCTGAAATTCCAGAACCG  >  NZ_CP009273/4091011‑4091279

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