Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,075 A→T 40.9% 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,0750AT40.9% 12.2 / 21.6 34intergenic (+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 (10/6);  new base T (5/7);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

AAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCT  >  NC_002947/4339941‑4340222
                                                                                                                                        |                                                                                                                                                   
aaaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCGCAAGCGGGGGTGGCaccaccc                                                                                                                                       >  3:171936/1‑148 (MQ=255)
          ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCATTTTATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc                                                                                                                             <  2:20758/149‑1 (MQ=25)
            cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                           <  3:195554/149‑1 (MQ=25)
            cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                           <  4:39724/149‑1 (MQ=25)
            cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTCATTTAGCGGGCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                           <  1:56806/149‑1 (MQ=255)
                        gAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCt                                                                                                                <  4:63660/148‑1 (MQ=255)
                              gAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGac                                                                                                                                              >  1:39567/1‑112 (MQ=255)
                              gAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGac                                                                                                                                              <  2:39567/112‑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            >  7:177074/1‑149 (MQ=25)
                                                                                                                            ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa            >  8:179866/1‑149 (MQ=25)
                                                                                                                            ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa            >  8:341307/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)
                                                                                                                                        |                                                                                                                                                   
AAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCT  >  NC_002947/4339941‑4340222

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: