Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 318,879 C→G intergenic (+84/+161) mak → / ← araJ manno(fructo)kinase/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306318,8790CG98.2% 134.8 / ‑3.4 56intergenic (+84/+161)mak/araJmanno(fructo)kinase/predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (31/24);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

GAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAAC‑‑CGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAC  >  AP012306/318738‑318934
                                                                                                                                               |                                                       
gAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATa                                                   >  1:1183860/1‑150 (MQ=25)
gAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATa                                                   >  2:2460575/1‑150 (MQ=25)
gAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATa                                                   <  2:3559512/150‑1 (MQ=25)
 aaGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATaa                                                  >  2:4177132/1‑150 (MQ=25)
  aGCACGGTGATTCCATCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAg                                                 <  2:1389858/150‑1 (MQ=21)
  aGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAg                                                 >  1:1068287/1‑150 (MQ=25)
  aGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAg                                                 >  2:3919399/1‑150 (MQ=25)
  aGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAg                                                 <  1:2143977/150‑1 (MQ=25)
    cACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGc                                               >  1:927389/1‑150 (MQ=25)
    cACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGc                                               <  2:2449113/150‑1 (MQ=21)
     aCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCg                                              >  1:1677722/1‑150 (MQ=25)
     aCGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCg                                              >  1:449245/1‑150 (MQ=25)
      cGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGt                                             >  2:3927270/1‑150 (MQ=21)
       ggTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGtt                                            >  1:4033209/1‑150 (MQ=21)
        gTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTc                                           >  2:2849698/1‑150 (MQ=21)
         tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCa                                          >  1:2255391/1‑150 (MQ=21)
         tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCa                                          >  2:4279452/1‑150 (MQ=21)
         tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCa                                          <  1:1007218/150‑1 (MQ=21)
         tGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCa                                          <  1:4218230/150‑1 (MQ=21)
           aTTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcg                                        >  2:1658420/1‑150 (MQ=21)
            ttCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgc                                       <  1:2261836/150‑1 (MQ=21)
             tCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgcc                                      <  1:1238208/150‑1 (MQ=21)
              ccAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccg                                     <  2:76302/150‑1 (MQ=21)
              ccAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccg                                     <  1:590410/150‑1 (MQ=21)
               cAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccgc                                    >  2:1044179/1‑150 (MQ=21)
               cAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccgc                                    >  2:2395884/1‑150 (MQ=21)
               cAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCAcgccgc                                    <  2:2424510/150‑1 (MQ=21)
                aGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCa                                   <  1:65822/150‑1 (MQ=21)
                aGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCa                                   >  2:3052505/1‑150 (MQ=21)
                aGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCa                                   >  2:1751809/1‑150 (MQ=21)
                   ggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATcc                                >  2:836315/1‑150 (MQ=18)
                   ggcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATcc                                >  2:3525472/1‑150 (MQ=18)
                    gcgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCg                               >  1:670255/1‑150 (MQ=18)
                    gcgTACGCGGCCTGGCGTGTTTTTGGCCACAAGAGTAAAAAAGGTAGGAAATGGGCGATTAATGCTTGATGCGACGTTTGCCGCGTCTTACCAGGCCTGCAAAAGGTGCCAGAACCGCGTAGGGCGGATAGGGCGTTCACGCCGCATCCg                               <  2:2896046/150‑1 (MQ=9)
                     cgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCgg                              <  1:2693247/150‑1 (MQ=18)
                     cgTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCgg                              <  2:3725669/150‑1 (MQ=18)
                      gTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGc                             >  2:265484/1‑150 (MQ=18)
                       tACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCa                            <  2:2730831/150‑1 (MQ=18)
                         cgcgGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAAt                          >  1:1718982/1‑150 (MQ=18)
                         cgcgGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAAt                          >  1:2661267/1‑150 (MQ=18)
                          gcggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGACAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATa                         <  1:3050829/150‑1 (MQ=17)
                           cggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                        <  1:3364160/150‑1 (MQ=18)
                           cggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATaa                        <  1:1517299/150‑1 (MQ=17)
                            ggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAg                       <  1:2752871/150‑1 (MQ=17)
                            ggcgCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAg                       <  1:1080375/150‑1 (MQ=17)
                             gcgcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGt                      <  1:2755962/150‑1 (MQ=17)
                              cgcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGta                     >  2:2538109/1‑149 (MQ=17)
                              cgcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGta                     <  2:657878/150‑2 (MQ=17)
                               gcTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaa                    >  1:2214070/1‑148 (MQ=17)
                                  gCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatga                 >  1:2236075/1‑145 (MQ=17)
                                   cGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatgag                >  2:3558028/1‑144 (MQ=14)
                                    gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatgagc               >  2:1658378/1‑143 (MQ=14)
                                    gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatgagc               >  1:1409631/1‑143 (MQ=14)
                                    gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatgagc               <  2:434297/150‑8 (MQ=14)
                                    gTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTGCAAAATGTGCCAGAACCGCGTAGGGCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGtaatgagc               >  2:168586/1‑143 (MQ=14)
                                               cacaAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAAC‑‑CGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAc  <  1:3385242/150‑1 (MQ=255)
                                                                                                                                               |                                                       
GAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATTGGCGCATCATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAAGGTGCCAGAAC‑‑CGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATAAGTGCTCCGATGCCTGATGCGAC  >  AP012306/318738‑318934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: