Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 550,552 C→G 13.8% intergenic (+30/‑112) tufB → / → rpsJ elongation factor Tu‑B/30S ribosomal protein S10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947550,5520CG13.8% 132.8 / 14.4 51intergenic (+30/‑112)tufB/rpsJelongation factor Tu‑B/30S ribosomal protein S10
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base G (3/4);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAG  >  NC_002947/550403‑550696
                                                                                                                                                     |                                                                                                                                                 
aacaTTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGCAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTc                                                                                                                                                   >  1:1184977/1‑150 (MQ=255)
    tnCAGATGACTGTTACCCTGATCAAGACCATCGCANTGGAAGACNGNCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCg                                                                                                                                               >  1:140025/1‑150 (MQ=255)
       aGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCgggg                                                                                                                                            >  1:153004/1‑150 (MQ=255)
           gACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGG‑ttt                                                                                                                                       <  1:1108767/150‑1 (MQ=255)
              tGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑tttttt                                                                                                                                    >  1:469501/1‑150 (MQ=255)
                   cccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTAtt                                                                                                                               <  2:167893/150‑1 (MQ=255)
                    ccTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTg                                                                                                                              >  2:1414453/1‑150 (MQ=255)
                      tGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTggg                                                                                                                            <  1:515378/150‑1 (MQ=255)
                       gATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGt                                                                                                                           <  1:204469/150‑1 (MQ=255)
                           aggACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGG‑ttttttcaat                                                                                                                                >  2:892810/3‑137 (MQ=255)
                           aggACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTGAGCGGGGG‑ttttttcaat                                                                                                                                <  1:892810/139‑5 (MQ=255)
                           aaGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGac                                                                                                                       >  1:1035507/1‑150 (MQ=255)
                              aCCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGa                                                                                                                        >  1:1104157/1‑146 (MQ=255)
                              aCCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGa                                                                                                                        <  2:1104157/146‑1 (MQ=255)
                                    ggagtGGGAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTAGGGGCGGCCCCGGGGTTGTGACCAAAATCATCGAGTAAGCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGT‑TTGTTTTTATTGGGTTGACACTAATTgg                                                                                                              <  2:39615/146‑1 (MQ=37)
                                    gCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTgg                                                                                                              <  2:153004/150‑1 (MQ=255)
                                     cAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTggg                                                                                                             <  2:687063/150‑1 (MQ=255)
                                        tGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGc                                                                                                          >  1:871245/1‑150 (MQ=255)
                                            aGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCt                                                                                                      >  1:1393049/1‑150 (MQ=255)
                                             gACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCta                                                                                                     <  2:871245/150‑1 (MQ=255)
                                                           cGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctc                                                                                       <  2:1393049/150‑1 (MQ=255)
                                                            gCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGCGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGcctcc                                                                                      >  1:1354718/1‑150 (MQ=255)
                                                                            gTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                      <  2:140025/150‑1 (MQ=255)
                                                                             tCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGGTTTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGt                                                                      <  2:280778/150‑1 (MQ=255)
                                                                                            cGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACa                                                      >  1:476534/1‑150 (MQ=255)
                                                                                                gTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCt                                                  >  2:578780/1‑150 (MQ=255)
                                                                                                 ttGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTg                                                 >  1:657067/1‑150 (MQ=255)
                                                                                                     ggCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGt                                             >  2:638754/1‑150 (MQ=255)
                                                                                                      gCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGt                                                              >  1:1256417/1‑132 (MQ=255)
                                                                                                      gCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGt                                                              <  2:1256417/132‑1 (MQ=255)
                                                                                                        cAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTg                                          <  2:1354718/150‑1 (MQ=255)
                                                                                                           aaTCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                       <  2:1035507/150‑1 (MQ=255)
                                                                                                           aaTCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGaaa                                       <  2:657067/150‑1 (MQ=255)
                                                                                                                tcGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCaa                                  <  1:1414453/150‑1 (MQ=255)
                                                                                                                 cGAGTAATCACTCGACCGATTGAAAATACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAAt                                 <  1:233914/150‑1 (MQ=255)
                                                                                                                 cGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAAt                                 <  1:638754/150‑1 (MQ=255)
                                                                                                                  gAGTAATCACTCGACCGATTGAAAAAGCCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTggg                                                          <  1:16912/124‑1 (MQ=255)
                                                                                                                  gAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTggg                                                          >  2:16912/1‑124 (MQ=255)
                                                                                                                          aCTCCACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGACTATATAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCa                        <  1:799288/150‑1 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTggg                                                                                                                            >  1:315429/1‑48 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTggg                                                                                                                            <  2:315429/48‑1 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTgg                                                           <  2:16003/113‑1 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTgg                                                           >  1:16003/1‑113 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                               >  2:425419/1‑125 (MQ=255)
                                                                                                                            tCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGa                                               <  1:425419/125‑1 (MQ=255)
                                                                                                                                         aaaaCCTCCGCTGAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAg         >  2:709174/1‑150 (MQ=255)
                                                                                                                                          aaaCCCCCGCTGAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAgg        >  2:163886/1‑150 (MQ=255)
                                                                                                                                            aCCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGtt      <  2:476534/150‑1 (MQ=255)
                                                                                                                                               cccGATCAGTGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGTGTCGTCTATAGAATCACCCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGCGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGaa   <  2:469501/150‑1 (MQ=255)
                                                                                                                                                ccGCTGAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg  <  1:163886/150‑1 (MQ=255)
                                                                                                                                                ccGCTGAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAg  <  1:709174/150‑1 (MQ=255)
                                                                                                                                                ccGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGTGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGAGTAGTGCGTCCGTTGGGAACAGCCTAGAGTCTGAAATCCAATGCAAAATCATCAAATCCGTATCAGGTTCAAg  >  2:1038974/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AACATTCAGATGACTGTTACCCTGATCAAGACCATCGCAATGGAAGACGGTCTGCGCTTCGCTATCCGTGAAGGCGGTCGTACCGTCGGCGCCGGCGTTGTGGCCAAAATCATCGAGTAATCACTCGACCGATTGAAAAAACCCCCGCTCAGCGGGGG‑TTTTTTTTATTGGGTTGACACTAATTGGGGGCGTCTATAGAATCACGCCTCCTTTTAGCGGGCGTAGTGCGTCCGTTGGGAACAGCCTGGAGTCTGAAATCCAATGCAAAATCAGCAAATCCGTATCAGGTTGAAG  >  NC_002947/550403‑550696

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