Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,128 C→T 8.8% intergenic (+69/‑262) ribAB‑II → / → PP_3814 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,1280CT 8.8% 125.0 / 2.1 45intergenic (+69/‑262)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base T (2/2);  total (25/20)
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.83e-01

GAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCG  >  NC_002947/4339984‑4340274
                                                                                                                                                |                                                                                                                                                    
gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCg                                                                                                                                                 <  1:1123259/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  2:1486152/150‑1 (MQ=255)
            gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  2:1646090/150‑1 (MQ=255)
             cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                    <  2:892363/150‑1 (MQ=255)
             cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCggg                                                                                                                                    <  1:451819/150‑1 (MQ=255)
                      tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGc                                                                                                                           <  2:58269/150‑1 (MQ=21)
                      tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                   <  1:1031801/126‑1 (MQ=255)
                      tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                                   >  2:1031801/1‑126 (MQ=255)
                        tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACGCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTc                                                                                                                         >  1:2108824/1‑150 (MQ=21)
                        ggggCTGGAGGTGGGGCAGAGCAAAACGCACGAGCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTc                                                                                                                         <  2:2108824/147‑1 (MQ=17)
                        cgtgCTCGATGTTGTGCAGAGCATCCCGTTCGTCCCGGGGATGGATGGCTGAGAACACACATGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTc                                                                                                                         <  2:1117727/147‑1 (MQ=38)
                           gCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGAAACCACCGACGCATTCGCGGGTTAACCCGCTCCTa                                                                                                                      >  1:496332/1‑150 (MQ=255)
                                        cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGt                                                                                                         <  2:496332/150‑1 (MQ=255)
                                                        aCCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCGGTAGGAGCGGGTGTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCgg                                                                                         >  1:1869416/1‑150 (MQ=255)
                                                                             aaCACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCtttt                                                                                                   <  2:364608/119‑1 (MQ=255)
                                                                             aaCACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCtttt                                                                                                   >  1:364608/1‑119 (MQ=255)
                                                                                    aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     >  2:1586390/1‑78 (MQ=255)
                                                                                    aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                     <  1:1586390/78‑1 (MQ=255)
                                                                                        aaTTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                      >  1:1153006/1‑57 (MQ=255)
                                                                                        aaTTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGacac                                                                                                                                                      <  2:1153006/57‑1 (MQ=255)
                                                                                          ttGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGt                                                                                      >  2:2200361/1‑119 (MQ=255)
                                                                                          ttGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGt                                                                                      <  1:2200361/119‑1 (MQ=255)
                                                                                                                  gCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTgg                                                                                                          >  2:1683458/1‑75 (MQ=255)
                                                                                                                  gCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTgg                                                                                                          <  1:1683458/75‑1 (MQ=255)
                                                                                                                   cGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGa                                                                                        <  2:358662/92‑1 (MQ=255)
                                                                                                                   cGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGa                                                                                        >  1:358662/1‑92 (MQ=255)
                                                                                                                     ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGcc‑‑cc                                                                    <  1:1082122/108‑1 (MQ=255)
                                                                                                                     ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGcc‑‑cc                                                                    >  2:1082122/1‑108 (MQ=255)
                                                                                                                     ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAatat                                 >  2:2061817/1‑143 (MQ=255)
                                                                                                                     ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAatat                                 <  1:2061817/143‑1 (MQ=255)
                                                                                                                      gggtACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTGAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                         >  1:1159073/4‑150 (MQ=255)
                                                                                                                      gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                         >  1:820424/1‑150 (MQ=255)
                                                                                                                      gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATccc                         >  1:621933/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGTCGGTGGTGTCACGGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCAGTAGaa                   >  2:826555/1‑149 (MQ=37)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                   >  2:466673/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                   >  2:524773/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCtt                                                >  2:1353114/1‑121 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCtt                                                <  1:1353114/121‑1 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                   >  1:2193020/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc                   >  2:1559797/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTAGAATACCATAATATGATATCCCGTAGAc                   >  2:1061749/1‑150 (MQ=255)
                                                                                                                            ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCCCTTGCCAAAATCTTGGAATACCATAATATGATATCCCTTacat                   >  1:216400/1‑147 (MQ=255)
                                                                                                                              gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa                               <  2:368545/136‑1 (MQ=255)
                                                                                                                              gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa                               >  1:368545/1‑136 (MQ=255)
                                                                                                                              gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt                 >  1:1788559/1‑150 (MQ=255)
                                                                                                                              gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCACAACGTT‑‑CGGTTCACCTTCAAGATCGGACGCGAAAATCTTGGAATAACATAATATGATATCCCGTAGACCt                 >  1:284715/1‑150 (MQ=16)
                                                                                                                                          gTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACctct                                                            >  1:2066688/1‑95 (MQ=255)
                                                                                                                                          gTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACctct                                                            <  2:2066688/95‑1 (MQ=255)
                                                                                                                                           tGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGc    <  2:1788559/150‑1 (MQ=255)
                                                                                                                                             acacCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCg  <  2:621933/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCC‑‑CCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCG  >  NC_002947/4339984‑4340274

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