Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 327,553 T→C 10.0% Y61C (TAC→TGC)  PP_0269 ← glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947327,5530TC10.0% 105.9 / 2.6 40Y61C (TAC→TGC) PP_0269glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTGAAGACAAGGCA  >  NC_002947/327408‑327692
                                                                                                                                                 |                                                                                                                                           
ggCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGtt                                                                                                                                         >  2:306203/1‑150 (MQ=255)
  cgcTTTCGTTCTTGGCGTTCACATCTACCAACGAGTGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGCCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGc                                                                                                                                       <  1:405564/150‑1 (MQ=255)
     tttCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGTATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTc                                                                                                                                    <  2:69300/150‑1 (MQ=255)
     tttCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTc                                                                                                                                    >  1:396244/1‑150 (MQ=255)
      ttCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCg                                                                                                                                   >  1:27414/1‑150 (MQ=255)
                       cGTAGACCAACCAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCttt                                                                                                                  <  1:24165/150‑1 (MQ=255)
                          aGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGgc                                                                                                               <  2:414466/150‑1 (MQ=255)
                                    gcgTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCAtt                                                                                                     <  2:189243/150‑1 (MQ=255)
                                        tGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAga                                                                                                 <  2:260966/150‑1 (MQ=255)
                                             gCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGAATTGGGATAGGTCATTTAgagcgac                                                                                            >  2:315039/1‑150 (MQ=39)
                                                    aaCGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTg                                                                                     <  1:306203/150‑1 (MQ=255)
                                                      cGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGag                                                                                   >  1:185160/1‑150 (MQ=255)
                                                          aGCTTGTCGTCGTCACATTAGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGAACCGCATGGTTGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGagcagc                                                                               >  2:430342/1‑150 (MQ=255)
                                                              tGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCgg                                                                           >  1:527741/1‑150 (MQ=255)
                                                                                    gTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAg                                                     >  2:166420/1‑150 (MQ=255)
                                                                                        tGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAg                                                 >  1:542556/1‑150 (MQ=255)
                                                                                        tGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAg                                                 <  1:464397/150‑1 (MQ=255)
                                                                                        tGGCGGACTTCCGGGCTGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCAGGGTGACTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTAACGGCCACAGCCAg                                                 <  2:26554/150‑1 (MQ=255)
                                                                                               tttCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGtt                                          >  2:1390/1‑150 (MQ=255)
                                                                                                ttCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTg                                         <  2:27414/150‑1 (MQ=255)
                                                                                                         gATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTa                                                                                                                                            <  1:33927/42‑1 (MQ=255)
                                                                                                         gATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTa                                                                                                                                            >  2:33927/1‑42 (MQ=255)
                                                                                                          aTGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAgg                               >  1:542154/1‑150 (MQ=255)
                                                                                                          aTGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAgg                               >  1:395978/1‑150 (MQ=255)
                                                                                                              tttCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCa                           <  1:166420/150‑1 (MQ=255)
                                                                                                                    tCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTgg                     <  2:527741/150‑1 (MQ=255)
                                                                                                                    tCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTgg                     >  1:324105/1‑150 (MQ=255)
                                                                                                                       ggTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTgg                                        <  2:17068/128‑1 (MQ=255)
                                                                                                                       ggTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTgg                                        >  1:17068/1‑128 (MQ=255)
                                                                                                                           aCGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGt              >  1:424995/1‑150 (MQ=255)
                                                                                                                            cGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTg             <  2:396244/150‑1 (MQ=255)
                                                                                                                                     cccAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTGAAGACAAgg    <  1:1390/150‑1 (MQ=255)
                                                                                                                                     cccAGGATCGGGTAGTTGCGTCCCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTGAAGACAAgg    <  1:404421/150‑1 (MQ=255)
                                                                                                                                     cccAGGATCGGGCAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTAcc                                                             >  2:147390/1‑93 (MQ=255)
                                                                                                                                     cccAGGATCGGGCAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTAcc                                                             <  1:147390/93‑1 (MQ=255)
                                                                                                                                     cccAGGATCGGGCAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTAcc                                                             >  2:147145/1‑93 (MQ=255)
                                                                                                                                     cccAGGATCGGGCAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTAcc                                                             <  1:147145/93‑1 (MQ=255)
                                                                                                                                       cAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTGAAGACAAGGCa  <  2:424995/150‑1 (MQ=255)
                                                                                                                                         ggATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGaaa                <  1:177502/134‑1 (MQ=255)
                                                                                                                                         ggATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGaaa                >  2:177502/1‑134 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACCTTCGATCAGGTACTGGCGGATTTCCGGGCGGATGGTTTCGATCAGGTAACGGATATTGCCCAGGATCGGGTAGTTGCGTCGCACCGCATGGTGGCTTTGGCGCAGGTCATTCAGACCGACCAGGCTGAGCAGCGCGGTGATCAGCGTTACCGGCCACAGCCAGGCATGTTGGGCCAGCAGGGGCAGGCTGGCGAAAGTGAAGACAAGGCA  >  NC_002947/327408‑327692

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