Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,120 G→C 11.1% intergenic (+61/‑270) 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,1200GC11.1% 107.3 / 5.4 36intergenic (+61/‑270)ribAB‑II/PP_38143,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base C (2/2);  total (17/19)
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.93e-01

AATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGC  >  NC_002947/4339986‑4340268
                                                                                                                                      |                                                                                                                                                    
aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                       <  2:1288391/150‑1 (MQ=255)
     gcgGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                  >  1:658853/1‑150 (MQ=255)
       aaggtctacgggatatcatattatggtattccaagattttggcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                  <  2:345600‑M1/56‑1 (MQ=255)
       aaggtctacgggatatcatattatggtattccaagattttggcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                  >  1:345600‑M1/93‑148 (MQ=255)
        aggtctacgggatatcatattatggtattccaagattttggcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                               >  1:458899‑M1/92‑150 (MQ=255)
         ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgcg                                                                                                                              <  2:658853/150‑1 (MQ=255)
                 ggCTATGAGCTGGAGGGGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATGGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                      <  2:728574/150‑1 (MQ=255)
                    tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                               <  1:1234958/138‑1 (MQ=255)
                    tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                               >  2:1234958/1‑138 (MQ=255)
                                      caagattttggcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATTCGAGGGTGACACCACCGACGCAtt                                                                                                                                  <  1:929121‑M1/56‑1 (MQ=255)
                                      caagattttggcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATTCGAGGGTGACACCACCGACGCAtt                                                                                                                                  >  2:929121‑M1/62‑117 (MQ=255)
                                            aTCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                     >  1:1116476/1‑92 (MQ=255)
                                            aTCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                     <  2:1116476/92‑1 (MQ=255)
                                                gcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                  >  2:1072358‑M1/52‑107 (MQ=255)
                                                gcaagaggtatggggcgaccggccatcactccgccctgtcaaaagcaccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCAtt                                                                                                                                  <  1:1072358‑M1/56‑1 (MQ=255)
                                                                ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                     <  1:882075/72‑1 (MQ=255)
                                                                ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAgg                                                                                                                                                     >  2:882075/1‑72 (MQ=255)
                                                                   tggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGc                                                                                                                   >  2:1278971/1‑103 (MQ=255)
                                                                   tggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGc                                                                                                                   <  1:1278971/103‑1 (MQ=255)
                                                                                acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAg                                                                                                                                                      >  2:1093899/1‑55 (MQ=255)
                                                                                acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAg                                                                                                                                                      <  1:1093899/55‑1 (MQ=255)
                                                                                  aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCCTTTGACAGGGCGGa                                                                                <  2:23409/123‑1 (MQ=255)
                                                                                  aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCCTTTGACAGGGCGGa                                                                                >  1:23409/1‑123 (MQ=255)
                                                                                               ccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGACACCACCGAc                                                                                                                                       >  2:1267636‑M1/5‑55 (MQ=255)
                                                                                               ccaaGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGTGACACCACCGAc                                                                                                                                       <  1:1267636‑M1/51‑1 (MQ=255)
                                                                                                     GCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTc                                                                                                                                 <  2:1158665‑M1/55‑1 (MQ=255)
                                                                                                     GCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTc                                                                                                                                 >  1:1158665‑M1/1‑55 (MQ=255)
                                                                                                         TGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTc                                                                   >  1:899160‑M1/1‑113 (MQ=255)
                                                                                                         TGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTc                                                                   <  2:899160‑M1/113‑1 (MQ=255)
                                                                                                                GCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata                       <  1:962898‑M1/150‑1 (MQ=255)
                                                                                                                  GGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATc                     <  1:728013‑M1/150‑1 (MQ=255)
                                                                                                                    GTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCa                                                             <  1:865524‑M1/108‑1 (MQ=255)
                                                                                                                    GTTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCa                                                             >  2:865524‑M1/1‑108 (MQ=255)
                                                                                                                              GAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGat                        >  2:937832‑M1/1‑135 (MQ=255)
                                                                                                                              GAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGat                        <  1:937832‑M1/135‑1 (MQ=255)
                                                                                                                                    GAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGc                                                                                   >  1:188555‑M1/1‑70 (MQ=255)
                                                                                                                                    GAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGc                                                                                   <  2:188555‑M1/70‑1 (MQ=255)
                                                                                                                                    GAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataa                              >  2:1349256‑M1/1‑123 (MQ=255)
                                                                                                                                    GAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGCAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataa                              <  1:1349256‑M1/123‑1 (MQ=255)
                                                                                                                                     aGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGa                                                                            >  2:235459/1‑76 (MQ=255)
                                                                                                                                     aGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGa                                                                            <  1:235459/76‑1 (MQ=255)
                                                                                                                                     aGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGc  >  2:714489/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                                    
AATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGC  >  NC_002947/4339986‑4340268

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