Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,075 A→T 23.4% 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,0750AT23.4% 34.0 / 8.3 28intergenic (+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/8);  new base T (4/2);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT  >  NC_002947/4339929‑4340216
                                                                                                                                                    |                                                                                                                                              
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGa                                                                                                                                                <  5:170573/147‑1 (MQ=255)
gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGa                                                                                                                                                <  7:263811/147‑1 (MQ=255)
  gccgccGCTGACATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGGTCGACCCGGTCATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGt                                                                                                                                            <  3:306110/149‑1 (MQ=255)
      ccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGCTGCTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATCCGAGGGAGaca                                                                                                                                        >  1:337004/1‑149 (MQ=255)
         cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGAcacca                                                                                                                                     <  8:93747/149‑1 (MQ=255)
           gAAATAGGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTGGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGAcaccac                                                                                                                                    <  2:337004/148‑1 (MQ=255)
           gAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTTACCCGCCAGGGCGGGGGTGaccacc                                                                                                                                     >  1:156739/1‑145 (MQ=255)
              aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATTGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                <  2:227451/149‑1 (MQ=255)
                      ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                        <  4:265414/149‑1 (MQ=255)
                                                  gCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTCCCGCGAATGCGAGGGTGACACCACCGCCGCCTTCGCGGGGAAACCCGCTCCTACAGGGCCCTTg                                                                                           >  8:294491/1‑149 (MQ=255)
                                                        cATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTACCGCGAATGCGAGGGTGTCACCACCGACGCATTAGCGGGTAAAACCCCGCCGAAAGGGCCCATGGTGcgt                                                                                     >  7:170997/1‑147 (MQ=38)
                                                         aTCCCGTTAGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTAACCCCGAAAGCGAGGGTGAAAACACCGAAGCATTAGCGGGTAAAACCCCTCCTACACGGCCCATGGTGCttt                                                                                    >  1:35346/1‑149 (MQ=37)
                                                                     ccGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCCCGAATGCGAGGGTGACACCACCGACGCAATCGCGGGTAAACCCCCTCCTACACGGCCCTTGGTGCTTTTGACAGGGCg                                                                          >  5:170980/1‑147 (MQ=255)
                                                                             ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTTACCGCGAAATCGAGGGTGACACCACCCACGCATTAGCGGGTAAAACCGCTACGTAAAGGGCCGTGGTGCTTTTTAAAAGGGAGAGTGATTGc                                                                >  2:247151/1‑149 (MQ=255)
                                                                                                 gaatttGATTTCGCGGGGCCGTGTTGGCGGGGG‑‑TTTCCCCGCTATGGCGTGGGTGACCCCCCCGACCCATGCGGGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGGGTGATGGCCGGTCGCCCCATACCTCTTg                                            <  8:170997/145‑1 (MQ=255)
                                                                                                            ccgctcccccTGTAGGCGCGGT‑‑TTACCCCGCTAATGCGCGGGTGACACCCCCGACGCATGCGCGGTAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                 <  5:2849/142‑1 (MQ=255)
                                                                                                             gCGGGGCCGTTTTGGTGGGGG‑‑TTTACCCGCTAATGCTAGGGTTACACCACCTACGAATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                 <  2:35346/148‑1 (MQ=255)
                                                                                                                     cTGTAGGAGCGGG‑‑TTTAACCGCGAATTCGAGGGTGACACCACCGACGCATTAGCGGGTTAACCCGCGACTAAACGGGCCTTGGTGGTTTTTACAGGGCGGTGTGATGGCCGGTCGGCCCAAACCTCTTTCCAAAATCTTGGAATACCa                         >  2:60686/1‑148 (MQ=255)
                                                                                                                      tGTAGGAGCGGG‑‑TTTAAACGCGAAAGCGAGGGTGGCACCACCCACGCATTAGCGGGTATACCCGCTCCTACACGGCCCGTGGTGCGTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCCGACCTCATGGCAAAATCTTGGAATAACaa                        >  6:289579/1‑147 (MQ=255)
                                                                                                                         attcgCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCGAAAATCTTGGAATACCataata                    >  5:231067/5‑149 (MQ=255)
                                                                                                                          ggAGCGGT‑‑TTTACCCGCGATTTCGAGGGTGACACAACCGACGCATTCGGGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGAAAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata                    <  1:13363/148‑1 (MQ=39)
                                                                                                                                 g‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc            >  2:24714/1‑149 (MQ=255)
                                                                                                                                        cccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATAACCCGTa        <  7:294491/148‑1 (MQ=255)
                                                                                                                                         ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa      >  3:76886/1‑149 (MQ=25)
                                                                                                                                         ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa      >  3:73137/1‑149 (MQ=25)
                                                                                                                                         ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAAACttc                                 >  4:25888/1‑121 (MQ=18)
                                                                                                                                          cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa        >  2:124726/1‑146 (MQ=25)
                                                                                                                                          cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTa        <  1:124726/146‑1 (MQ=25)
                                                                                                                                            cgAATGCTAGGGTGACACCAACGACGCATTCGCGGGTAAACCCGATCATACAGGGCCGTTGGTGCTTATGAAAGGGCTGACTGATGGCCGCTCGCCCCATACCTCTTTACAAAATCATGGAAGACCATAATATGATATCCAGTCGAcat   >  8:138576/1‑147 (MQ=255)
                                                                                                                                             gAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCtt  >  6:166442/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTT  >  NC_002947/4339929‑4340216

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: