Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,089 G→C 30.5% intergenic (+74/‑257) 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,0890GC30.5% 38.0 / 24.4 36intergenic (+74/‑257)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base C (7/4);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

GGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCA  >  NC_002947/4339951‑4340233
                                                                                                                                            |                                                                                                                                                
ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCATTTTATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                                        <  2:20758/149‑1 (MQ=25)
  cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  4:39724/149‑1 (MQ=25)
  cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                                      <  3:195554/149‑1 (MQ=25)
  cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTCATTTAGCGGGCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                                      <  1:56806/149‑1 (MQ=255)
              gAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCt                                                                                                                           <  4:63660/148‑1 (MQ=255)
                             cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGc                                                                                                                            >  3:146754/1‑132 (MQ=255)
                                               ccGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCCCGAGTGCGAGGGTGACACCCCCCACGCAATAGCGGGTAAACCCGCTCCGACACGGCCCATGGTGCTTTTGACAGGGCg                                                                                           >  4:268835/1‑147 (MQ=37)
                                                             cTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCGCGAATGCGAGGGGGGCACCCCCGACGCATTAGCGGGTAAACCCCCACCTAAAGGGCCCGTGGTGGGTTTTTCAGGGCGGGGTGATTGGCGGtgg                                                                           >  8:201385/1‑147 (MQ=25)
                                                                       acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTTTACCCCCAAGGCCGGGGGGACCCCCCCCCCCCATTGCCGGGTAAACCCCCCCCAAAACGCCCCGGGGTGGGTTTTAAAGGGGGGGGGGATGGCCGGGGGCCCCAAAc                                                                   >  2:311367/1‑148 (MQ=255)
                                                                       acaACCAATTGATTTAGCGGGGCCCTGTAGGCGCGGGTTTACCCGCGAATCCTTTGGGTGACACCACCGACGCATTCGGGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATa                                                                    <  1:311367/147‑1 (MQ=37)
                                                                          aGCAATTGATTTCGCGGGGCCCGGTCGGCGCGGGTTTCCCCGCGAATGCGCGGGTGACACCACCGACGATTGCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCtt                                                              <  7:201385/149‑1 (MQ=38)
                                                                              aTTGATTTAGCGGGGCCCTGTCGGAGCGGGTTTACCCGCGAATGCCCGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCa                                                          <  5:296359/149‑1 (MQ=255)
                                                                                       gCGGGGCCCTGTCGTCGCGGGTTACCCCGCTAATGCGCGGGTTACACCACCGACCCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGCGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTgg                                                 <  1:266726/149‑1 (MQ=38)
                                                                                       gCGGGGCCCTGTAGGCGCGGGTTTCCCCGCGAGTGCGAGGGTGACACCACCGACGCATGCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                  <  3:268835/148‑1 (MQ=38)
                                                                                            gCCCTTTCGGCGGGGTTTCACCCGCTATTGCGTGGGTCACCCCACCGCCTCTTTGGGGTTAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                                            <  4:171936/149‑1 (MQ=255)
                                                                                            ccccTGTAGGAGCGGGTTAACCCGCTAAGGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                                            <  7:117239/148‑1 (MQ=255)
                                                                                                     gAGCGGGTTTACCCGCGAATGCGAGGGTGACACCCCCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGCTGGCCGGGCGCCCCATACCTCTTGCCAAAATCTTGGAATCCCATAatat                                    >  5:193837/1‑148 (MQ=255)
                                                                                                           ggTTTCCCCGCGAAGGCGCGGGGGACCCCCCCGCCGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGGGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                              <  6:193837/149‑1 (MQ=255)
                                                                                                           ggTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                              >  4:322898/1‑149 (MQ=25)
                                                                                                           ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                  >  4:127011/1‑129 (MQ=255)
                                                                                                           ggTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                                  <  3:127011/129‑1 (MQ=255)
                                                                                                            gTTTACCCGCGAATGCGTGGGTCACCCCAGGTCCGCATGCGGGTTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc                             <  3:73667/149‑1 (MQ=255)
                                                                                                                  ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  7:136664/1‑149 (MQ=25)
                                                                                                                  ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  8:341307/1‑149 (MQ=25)
                                                                                                                  ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  7:177074/1‑149 (MQ=25)
                                                                                                                  ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa                       >  8:179866/1‑149 (MQ=25)
                                                                                                                  ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTGTTGACAGGGCGGAGTGATGGCCGGGCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATACCGTAGa                       >  3:274741/1‑149 (MQ=255)
                                                                                                                    gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                     <  7:341307/149‑1 (MQ=25)
                                                                                                                    gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc                     >  4:131265/1‑149 (MQ=255)
                                                                                                                          gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAg                >  3:372628/1‑148 (MQ=255)
                                                                                                                             gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCt             >  5:177102/1‑149 (MQ=255)
                                                                                                                                        caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCg     <  2:89620/146‑1 (MQ=255)
                                                                                                                                        caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCa  >  2:29185/1‑149 (MQ=255)
                                                                                                                                        caccGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTATTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCg     >  1:89620/1‑146 (MQ=255)
                                                                                                                                        caccGACGCATACGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAg         >  2:54868/1‑142 (MQ=255)
                                                                                                                                        caccGACGCATACGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAg         <  1:54868/142‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCA  >  NC_002947/4339951‑4340233

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