Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,075 A→T 17.8% intergenic (+60/‑271) 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,0750AT17.8% 46.6 / 3.6 32intergenic (+60/‑271)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base T (3/3);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA  >  NC_002947/4339929‑4340218
                                                                                                                                                  |                                                                                                                                                 
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGa                                                                                                                                                   <  4:236786/147‑1 (MQ=255)
     gccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGa                                                                                                                                                   >  8:55889/1‑142 (MQ=255)
     gccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGa                                                                                                                                                   <  7:55889/142‑1 (MQ=255)
          tGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTggtgtcac                                                                                                                                       >  3:86417/1‑142 (MQ=255)
                      ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                           <  8:63414/149‑1 (MQ=255)
                                cTATGAGCTGGAGGTGGTGCAGACCATCCCGTCCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCGCGCAAGAGAGGGTGGCGCCACCCCCGCAATAGCGCGGTTAAcc                                                                                                                >  3:257882/1‑149 (MQ=255)
                                 tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCGCAATCGCAGGGGTCACCACCGCCCCCATAGCGGGTGAAcccc                                                                                                              >  5:203654/1‑148 (MQ=255)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                        >  4:303111/1‑138 (MQ=255)
                                   tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                        <  3:303111/138‑1 (MQ=255)
                                      gCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCGCGCATGCGAGGGTGACACCACCGCCGCATTCGCCGGTAAACCCGCAcc                                                                                                          >  8:266918/1‑148 (MQ=39)
                                       cTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCt                                                                                                         >  2:48802/1‑148 (MQ=255)
                                                   cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                          >  3:209543/1‑121 (MQ=255)
                                                   cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                          <  4:209543/121‑1 (MQ=255)
                                                             cGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTAACCCGAATTGGAGGGTGAACCCAACGCCCCAATAGCGCGTAAACAGCCTCCAAAAGGGCCCgtggtgggttttg                                                                                   >  5:322386/1‑141 (MQ=255)
                                                                  gACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCGCGAATGCGAGGGTGACACCACCCCCCCATTCACGGGGAAACCCGCACCGACACGGCCCTTGGTGCATTTTACAAGGc                                                                             >  4:272966/1‑149 (MQ=37)
                                                                               atggCTGAGAACACACAAGCAATTGATTTCGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCg                                                                <  1:90779/149‑1 (MQ=255)
                                                                               atggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTAACCCCGAAAGCCACGGTGACAACACCCACCCAATAGCGGGTAAAACCCCAACAACAAGGCCCATGGTGCTTTTTACAGGGGGGAGTTATGGCCg                                                                >  2:90779/1‑149 (MQ=255)
                                                                                       caaCACACATGAAATTTTTTTCGGGGGGCGGTGTCGTCGGGGTTTTACCCGCTAAGGGGTGGGTTACACCACGTACGTTTTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc                                                         <  6:322386/147‑1 (MQ=255)
                                                                                         acacacAAGCAATTGATTTAGCGGGGCCCTGTAGGCGGGGGTTCACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCAt                                                      <  1:167960/149‑1 (MQ=255)
                                                                                             acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTAGCGGGTTAACCCGCTCCTACAGGGGCCTTGGTGCGTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCCTACct                                                  >  5:356476/1‑149 (MQ=38)
                                                                                                        aTTTCGCGGGGCCCTGTCGGAGCGGGTTTACCCGCTAATGCGAGGGTGACACCACCGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGc                                                                  <  5:132622/122‑1 (MQ=255)
                                                                                                             gCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAAGGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                   <  1:48802/148‑1 (MQ=255)
                                                                                                             gCGGGGCCCTGTAGCAGCGGGATTACCCGCGAATGCGAGGGGGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                   <  7:266918/148‑1 (MQ=255)
                                                                                                              ggggcCCCTGTAGGCGCGGGTTACCCCGCTACGGCGAGGGGCACCCCCCCCCCCCTTGCGGGGTAAACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGa                                 <  3:128196/148‑1 (MQ=255)
                                                                                                                           gAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGGTTTTGCCAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCGTGCCAAAATCTTGGAATACCATAatat                     >  5:208254/1‑148 (MQ=255)
                                                                                                                              cgtgTTAACCCGCGAATGCCTGGGTGACACCACCTACCAATAGGGTTTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                 <  5:313151/146‑1 (MQ=255)
                                                                                                                                 gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAgg                                                                                                    >  1:146772/1‑63 (MQ=255)
                                                                                                                                 gTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAgg                                                                                                    <  2:146772/63‑1 (MQ=255)
                                                                                                                                       ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa        >  2:284952/1‑149 (MQ=25)
                                                                                                                                         gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc      >  6:234048/1‑149 (MQ=25)
                                                                                                                                          ctaaggCGAGGGTCACACCCCCGACGCATGCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt     <  6:208254/144‑1 (MQ=255)
                                                                                                                                             ttGCGTGGGTGACACCACCGACGAATTGGGGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa  <  3:351758/148‑1 (MQ=255)
                                                                                                                                             aTGGGAGGGTGACACCACCGACGAATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGa  <  3:272966/149‑1 (MQ=255)
                                                                                                                                                 gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGatat                <  4:289219/131‑1 (MQ=255)
                                                                                                                                                 gAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGatat                >  3:289219/1‑131 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGA  >  NC_002947/4339929‑4340218

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: