Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,787,816 A→C 17.1% F914V (TTT→GTT)  gltB ← Glutamate synthase [NADPH] large chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,787,8160AC17.1% 115.3 / 2.2 53F914V (TTT→GTT) gltBGlutamate synthase [NADPH] large chain
Reads supporting (aligned to +/- strand):  ref base A (37/5);  new base C (2/9);  total (39/14)
Fisher's exact test for biased strand distribution p-value = 2.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCTTC  >  ACB122/5787686‑5787948
                                                                                                                                   |                                                                                                                                    
tAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGc                                                                                                                                 >  3:195406/1‑136 (MQ=255)
    gCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGGCCCCCCAACACGGCCgg                                                                                                                              >  6:254552/1‑135 (MQ=255)
      aGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCCGCTTGGCGCCCTGGGCCACCTAGGTCTGCAGCACATCGGCGTTG‑CCCAGAGACTCCGGGGGCCCACCAAAACGGCCCGTgg                                                                                                                           >  3:251647/1‑136 (MQ=255)
           ggCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAACCGGCCGGTGGcccc                                                                                                                       >  2:350119/1‑133 (MQ=255)
           ggCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCAcc                                                                                                                      >  7:275085/1‑136 (MQ=255)
           ggCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCAc                                                                                                                       >  5:7354/1‑135 (MQ=255)
                tCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCt                                                                                                                  >  2:337513/1‑135 (MQ=255)
                     ccGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACCCCAAAACGGCCGGTGGCCCCCTGCTTGATc                                                                                                             >  1:75092/1‑135 (MQ=255)
                      ggttgccgttgccGCCCGCCACCGGGCCGCCCCCCCCGGGGGGGGGGCCGGGGGCAACGTTTAGCGGCAGCACTGCGGCGTGC‑CCCAGATACGCGGGGGCCACACAAACACGGCCGGTGGCCACCTGCTTGATCt                                                                                                            <  2:75092/134‑1 (MQ=255)
                      cgtggcccTTCCCGCCCGGCACGGGGCGGCCCCCCCCCGGGTGGGCGCCCGGGGAAACGTTGCGCGGCAGCA‑GTGGGGGTGGCACCAGACCCTCGGGGGTCACCCCAACACGGCCGGTGGCCACCTGCTTGATCt                                                                                                            <  2:123695/129‑1 (MQ=255)
                         tgcccTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGCTCTGCCGCACTTCGGCGTTG‑CCCAGATCCTCCGGGGTCACACCAACACGGCCGGGGGCCACCGGCTTGCTCGTc                                                                                                          >  8:149926/4‑134 (MQ=255)
                            ccTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGc                                                                                                      >  2:69475/1‑135 (MQ=255)
                               ttCCGCCCGGCAACTGGCCGCCCTCACCCGGGTTGGCGCCCGGGGAAACCTTGATCTCCCGCAGTGGGGCGTTG‑ACCTGATACTGGGGGGGCCCACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGctct                                                                                                   <  3:164336/135‑1 (MQ=255)
                                 ccgcccgcCAACGGGCCGCCCTCCCCCGGCTGGGCGCCCTGGGCAACCTTGATCTGCAGCACTGCGGCGTTC‑ACCAGATACTCCGGGGACACACCAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGctct                                                                                                   <  6:116065/133‑1 (MQ=255)
                                 ccGCCCGGCACCTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTGCGGGGTTG‑ACCAGATACTCCGGGGACACACAAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTg                                                                                                 <  8:177840/135‑1 (MQ=255)
                                 ccGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTg                                                                                                 >  6:312611/1‑135 (MQ=255)
                                 ccGCCCGGAATATGGCAGCCCTCACGTGGCTTGCCTCCCTGGGAAAGCTTGATCTGTATCAATGCGGTTTTT‑ACCATATTCTCGGGGGTCACAAAAAAAGGGACGGTGGGCACCTGCTTGATCTTCGAGctgt                                                                                                   <  6:197530/133‑3 (MQ=255)
                                          aaCTGGCCGCCCTCACCCGCCTTGGCGCGCTGGGCAACCTTGATCTGCTGGGCTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCAGCTGCTTGATCTTCGAGCTCTTGATGGTGccc                                                                                        >  1:308897/1‑134 (MQ=255)
                                               cccgccCTCACCCGGCTGGGCGCCCTGGCAAACCTGGATCTGCAGCACTTCGGCGTTC‑ACCACATACTCCGGGGCCACACCAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAgcg                                                                                   <  5:103646/134‑1 (MQ=255)
                                                        aCCCGGGTTGGCGCCCGGGCCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCGGGGGCAACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCt                                                                         <  8:1181/136‑1 (MQ=255)
                                                                 ggCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCgccgcc                                                                 >  6:158920/1‑135 (MQ=255)
                                                                 ggCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCt                                                                >  8:71614/1‑136 (MQ=255)
                                                                 ggCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTAACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACTGGTCTTCGCCGCCt                                                                >  3:54623/1‑136 (MQ=255)
                                                                 ggCGCCCTGGGCAACCTTGATCTGCAGCAATTCGGCGTTG‑ACAAGATACTCCGGGGTCACACCAAAACGGACGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCTTAGAGCGACGGGTCTTcgccgac                                                                 >  6:287671/1‑133 (MQ=255)
                                                                  gcgcCCTGGCCAACCTTGATCTCCAGCCCTTCGGCGTTC‑ACCAGATCCTCGGGGGTCACACCAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCtt                                                               <  3:365176/136‑1 (MQ=255)
                                                                  gcgcCCGGGGCAACGTTGAGCTGCAGCACTTCGGCGTTC‑ACCAGATACTCGGGGGTCACACCAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCtt                                                               <  8:275085/136‑1 (MQ=255)
                                                                     ccgtGGGCAACGTTCATGTGCAGCACTTCGGCGTTG‑ACATTAACCTCCGGGGCCACACAAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCa                                                             <  2:62567/132‑1 (MQ=255)
                                                                     cccTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCa                                                             >  6:63779/1‑135 (MQ=255)
                                                                     cccTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCa                                                             >  1:43307/1‑135 (MQ=255)
                                                                          ggCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACAAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCAATAGAGCGACTGCTCTTCGCCGCCTTCACCCGa                                                        >  1:222498/1‑135 (MQ=255)
                                                                           gCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAg                                                       >  5:243787/1‑135 (MQ=255)
                                                                           gCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAg                                                       >  1:281378/1‑135 (MQ=255)
                                                                           gCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAg                                                       >  1:308131/1‑135 (MQ=255)
                                                                           gCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGt                                                      >  7:24732/1‑136 (MQ=255)
                                                                           gCAACCTTGATCTGCAGCACTTCGGCGTTG‑AACAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGt                                                      >  7:290646/1‑136 (MQ=255)
                                                                                         caccaCTGCGGCGTGC‑ACCAGATACCCCGGGGTCACACAAACACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTTCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGccc                                        <  7:19037/133‑1 (MQ=255)
                                                                                            caCTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAgg                                     >  4:164210/1‑136 (MQ=255)
                                                                                            caCTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAgg                                     >  7:133802/1‑136 (MQ=255)
                                                                                                               aTACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGc                    <  2:12055/135‑1 (MQ=255)
                                                                                                                  cTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAg                 >  6:91213/1‑135 (MQ=255)
                                                                                                                       gggTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTc           >  4:28168/1‑136 (MQ=255)
                                                                                                                         gTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCg          >  2:289006/1‑135 (MQ=255)
                                                                                                                         gTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCATGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCTGCCCGCGCTTCGt         >  8:368225/1‑136 (MQ=255)
                                                                                                                          cccccccAACACGGCCGGGGGCCACCTGCTTGAGCTTCGAGCTCTTGATGGGGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGt         <  1:192/130‑1 (MQ=255)
                                                                                                                               ccAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGAGGGCt    >  6:255661/1‑135 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTGGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCTTc  >  7:262767/1‑136 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCCCCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCTTc  >  7:159173/1‑136 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCtt   >  2:11404/1‑135 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCtt   >  5:248752/1‑135 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCtt   >  6:144283/1‑135 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCTTc  >  8:358018/1‑136 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTGGATGGTGCCATAGCGCGACGCCACTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGTCCAGCGCTGCCTGGGCTTc  >  8:330579/1‑136 (MQ=255)
                                                                                                                                cAAAACGGCCGGTGGCCACATGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGTCCTCGGCCAGCGCTTCGTGGGCtt   >  1:293701/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                    
TAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTG‑ACCAGATACTCCGGGGTCACACCAAAACGGCCGGTGGCCACCTGCTTGATCTTCGAGCTCTTGATGGTGCCATAGCGCGACGGGTCTTCGCCGCCTTCACCCGAGTTGGAGCGTGCGCCCAGGCGGTTCATGGCCTCGGCCAGCGCTTCGTGGGCTTC  >  ACB122/5787686‑5787948

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