Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,359,171 G→C 7.7% V1483L (GTG→CTG)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,359,1710GC 7.7% 176.8 / 5.7 52V1483L (GTG→CTG) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base C (2/2);  total (22/30)
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.23e-01

GCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAG  >  NC_000913/3359023‑3359314
                                                                                                                                                    |                                                                                                                                               
gCATGTGCAGCATACCGGCTCTCAGCGCGGTGCAGAGCTTCTGGCGAACGGGCCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGAGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGCTGcccg                                                                                                                                                 >  3:453665/1‑149 (MQ=255)
 cATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGt                                                                                                                                                <  2:1039/149‑1 (MQ=255)
   tgtgCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGc                                                                                                                                              <  3:289572/149‑1 (MQ=255)
     tgCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAg                                                                                                                                            >  2:263563/1‑149 (MQ=255)
         gcaTACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTAGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGgcgc                                                                                                                                        >  2:407619/1‑149 (MQ=255)
               cGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAg                                                                                                                                  >  3:25527/1‑149 (MQ=255)
                 gCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAggg                                                                                                                                <  3:238532/149‑1 (MQ=255)
                 gCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAggg                                                                                                                                <  1:47237/149‑1 (MQ=255)
                 gCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAggg                                                                                                                                <  3:142919/149‑1 (MQ=255)
                    ctcAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTa                                                                                                                             <  3:392986/149‑1 (MQ=255)
                    ctcAGCGCGGTGAAGAGATTCTGGCGAACTGGCCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCGCTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGAAGGCGCAGTAAGGGGTa                                                                                                                             <  4:453665/149‑1 (MQ=255)
                        gcgcGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGgcg                                                                                                                                         <  1:15364/133‑1 (MQ=255)
                        gcgcGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGgcg                                                                                                                                         >  2:15364/1‑133 (MQ=255)
                            ggTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAAt                                                                                                                     <  3:15397/149‑1 (MQ=255)
                                  gagaTTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCa                                                                                                               >  4:205312/1‑149 (MQ=255)
                                         tGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGtt                                                                                                        <  3:96113/149‑1 (MQ=255)
                                          ggCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGttt                                                                                                       <  2:1227/149‑1 (MQ=255)
                                             gAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATc                                                                                                    >  3:246786/1‑149 (MQ=255)
                                                     cAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATc                                                                                            >  4:52246/1‑149 (MQ=255)
                                                      aaCCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCg                                                                                           <  1:407619/149‑1 (MQ=255)
                                                        ccTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGAc                                                                                         <  3:205312/149‑1 (MQ=255)
                                                              ccACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAg                                                                                   >  2:255277/1‑149 (MQ=255)
                                                                aCTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGc                                                                                                                           <  3:470945/107‑1 (MQ=255)
                                                                aCTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGc                                                                                                                           >  4:470945/1‑107 (MQ=255)
                                                                      tttGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGa                                                                           >  2:244999/1‑149 (MQ=255)
                                                                               gTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCaaa                                                                  <  3:52246/149‑1 (MQ=255)
                                                                               gTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCaaa                                                                  >  1:277680/1‑149 (MQ=255)
                                                                                    aCCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGcaa                                                                                                                         <  3:415872/89‑1 (MQ=255)
                                                                                    aCCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGcaa                                                                                                                         >  4:415872/1‑89 (MQ=255)
                                                                                       gAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGc                                                          <  1:297065/149‑1 (MQ=255)
                                                                                              aGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCCTGc                                                                                                                                              >  2:256748/1‑58 (MQ=255)
                                                                                              aGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCCTGc                                                                                                                                              <  1:256748/58‑1 (MQ=255)
                                                                                               gTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGACAAAGAAACCGCTGAAGATc                                                  <  4:25527/149‑1 (MQ=255)
                                                                                                     tAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATCAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCaaa                                            >  4:79602/1‑149 (MQ=255)
                                                                                                              tgctgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGttt                                   <  1:263563/149‑1 (MQ=255)
                                                                                                              tgctgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGttt                                   >  2:64014/1‑149 (MQ=255)
                                                                                                                ctgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAAt                                                                                                 >  2:67714/1‑85 (MQ=255)
                                                                                                                ctgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAAt                                                                                                 <  1:67714/85‑1 (MQ=255)
                                                                                                                 tgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTttgtt                                <  2:51061/149‑1 (MQ=255)
                                                                                                                 tgGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTttgtt                                <  2:20270/149‑1 (MQ=255)
                                                                                                                     tCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGaaa                            <  1:255277/149‑1 (MQ=255)
                                                                                                                         cGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTa                        >  2:331453/1‑149 (MQ=255)
                                                                                                                         cGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTa                        >  4:273359/1‑149 (MQ=255)
                                                                                                                            aGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGa                     >  1:172428/1‑149 (MQ=255)
                                                                                                                                     gcAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTc            >  2:397783/1‑149 (MQ=255)
                                                                                                                                     gcAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTc            <  1:397783/149‑1 (MQ=255)
                                                                                                                                         cTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGaa        <  1:64014/149‑1 (MQ=255)
                                                                                                                                         cTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGaa        <  1:244999/149‑1 (MQ=255)
                                                                                                                                          tcagctGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAag       >  4:236857/6‑149 (MQ=255)
                                                                                                                                          tGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAag       <  2:277680/149‑1 (MQ=255)
                                                                                                                                           ttactGCGCCTGCACGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAAc                                                             >  2:57446/5‑94 (MQ=255)
                                                                                                                                           ttactGCGCCTGCACGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAAc                                                             <  1:57446/90‑1 (MQ=255)
                                                                                                                                               tGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAaggccag  <  2:317298/149‑1 (MQ=255)
                                                                                                                                                    gTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTAc                       <  3:218791/123‑1 (MQ=255)
                                                                                                                                                    gTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTAc                       >  4:218791/1‑123 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAACTGGTCAACCTTCGCCACTAAATTTGCGCTGGTTAAACCGAAGTCCAGTGATGTAAAAGCACTGCTGGGTCACCGTAGTCGTAGCGCAGCTGAGTTGCGCGTGCAGGCGCAGTAAGGGGTAGCAACAATGAGTCAGAATGTTTATCAATTTATCGACCTGCAGCGCGTTGATCCGCCAAAGAAACCGCTGAAGATCCGCAAAATTGAGTTTGTTGAAATTTACGAGCCGTTTTCCGAAGGCCAG  >  NC_000913/3359023‑3359314

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